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Abstract 



This article reviews our view on how chiral symmetry, its pattern of breaking and 
(3 \ restoration under extreme conditions manifest themselves in the nucleon, nuclei, nuclear 

matter and dense hadronic matter. We discuss how first-principle (QCD) calculations 
of the properties of finite nuclei can be effectuated by embedding the "standard nuclear 
p • physics approach (SNPA)" into the framework of effective field theories of nuclei that 

^ , incorporate chiral dynamics and then exploit the predictive power of the theory to accu- 

D \ rately compute such solar neutrino processes as the proton-proton fusion and the "hep" 

^ ' process and such cosmological nucleosynthesis process as thermal neutron-proton capture 

etc. The Brown-Rho (BR) scaling that implements chiral symmetry property of baryon- 
^ ■ rich medium is re-interpreted in terms of "vector manifestation" of hidden local symmetry 

a la Harada-Yamawaki. We present a clear direct evidence and a variety of indirect ev- 
idences for BR scaling in nuclear processes at normal nuclear matter density probed by 
weak and electromagnetic fields and at higher density probed by heavy-ion collisions and 
compact-star observables. We develop the notion of "broadband equilibration" in heavy- 
ion processes and sharpen the role of strangeness in the formation of compact stars and 
their collapse into black-holes. We revisit the Cheshire-Cat phenomenon first discovered 
in the skyrmion structure of baryons and more recently revived in the form of "quark- 
hadron continuity" in mapping low-density structure of hadrons to high-density structure 
of quarks and gluons and argue once more for the usefulness and power of effective field 
theories based on chiral symmetry under extreme conditions. It is shown how color-flavor 
locking in terms of QCD variables and hidden local symmetry in terms of hadronic vari- 
ables can be connected and how BR scaling fits into this "continuity" scheme exhibiting 
a novel aspect of the Cheshire Cat phenomenon. 
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1 INTRODUCTION 



In a recent paper, Morgenstern and Meziani reported a remarkable result of 
their extensive analysis of the longitudinal and transverse response functions in medium- 
weight and heavy nuclei and concluded that the available world data showed unequivocally 
the quenching of ~ 20% in the longitudinal response function and the only viable way 
to explain this quenching is to invoke the Brown-Rho (BR) scaling [0] formulated to 
account for the change of the strong-interaction vacuum induced by matter density. Partly 
motivated by this development and other indirectly related developments including the 
most recent one on color-flavor-locking and color superconductivity, we shall make in this 
review our overview of evidences that indicate both directly and indirectly that the BR 
scaling is indeed operative in finite nuclei and dense (and superdense) nuclear matter. We 
shall develop our arguments starting with the basic structure of the nucleon, then go to 
that of nuclei and of nuclear matter and finally to that of dense matter that one expects to 
create in relativistic heavy-ion collisions and find in the interior of highly compact stars. 
The principal theme will be that the chiral symmetry of quantum chromodynamics - its 
spontaneous breaking in free space and partial or full restoration in medium - underlies 
the common feature from elementary hadrons to complex dense systems. We admit that 
our views are not necessarily shared by others in all details and that some points may not 
be fully correct. We are however confident that the general theme that we have developed 
and we shall review here will survive the experimental test that is to come. 

It was recognized a long time ago that chiral symmetry, now identified as an es- 
sential ingredient of QCD associated with the light-mass up (u), down (d) and strange 
(s) quarks, with its "hidden" realization in Nambu-Goldstone mode plays an important 
role in nuclear physics. In fact, the emphasis on its role in a variety of nuclear processes 
predates even the advent of QCD proper and subsequent acceptance as the correct theory 
of strong interactions [Q. Among the time-tested and established is the role of chiral 
symmetry in exchange currents in electro-weak processes in nuclei [|, |^ and its subse- 
quent importance in probing nuclear structure at various electron machines, now out of 
operation, throughout the world Q. Since our early review on the subject 0, there has 
been continuous and significant evolution in the field specifically associated with the role 
of QCD in nuclear physics, more or less unnoticed by workers in other areas of physics. 
The evolution touched on the structure of the nucleon, the proton and the neutron, start- 
ing with the notion of a "little bag" the recognition of the preeminent importance of 
Goldstone pion clouds in both the nucleon structure and nuclear forces [0, §], followed by 
the resurrection of the Skyrme soliton model that contains baryons - which are fermions 
- in a bosonic field theory ||^, |10[ and then the emergence of the notion of "Cheshire Cat 
phenomenon." At present, the issue of both "spontaneous" breaking and restoration of 
chiral symmetry occupies one of the central positions in current activities in both nuclear 
and hadron physics communities. The development up to 1995 has been summarized in 
a recent monograph by Nowak, Rho and Zahed [pJ| . 

It is now becoming clearer in which direction the chiral symmetry of QCD will steer 
the next generation of nuclear and hadronic physicists, with the advent of new accelerators 
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such as the Jefferson Laboratory (JL) and the RHIC at Brookhaven - both of which are 
already operating - and ALICE/LHC at CERN - which is to come in a few years. On one 
hand, the electron machine at the JL will probe deeper into shorter distance properties 
of the nucleon-nucleon interactions, exploring how the chiral structure of the nucleon 
in medium changes from hadronic over to quark/gluon picture as shorter distances are 
probed. As will be repeatedly stressed, one does not expect any abrupt phase transitions 
along the way; as the change will most likely be continuous (for a multitude of reasons 
we will discuss), it will then be a matter of economy which language will be more efficient 
for given processes with given conditions, implying that there will be a continuous map 
between the various descriptions. This continuity will play a central role in our discussion 
and will later be referred, quite broadly, to as "Cheshire Cat Phenomenon." 

On the other hand, the heavy- ion machines at RHIC and ALICE/LHC will create 
systems at high temperature and density, mimicking the early Universe. At high temper- 
atures, QCD predicts a phase transition from chiral symmetry in the Goldstone mode to 
that in the Wigner mode with a consequent change in confinement /deconfinement. This 
is more or less confirmed on lattice in QCD. Indications are coming out already from 
available experimental data that such a transition has been seen. Matter at high density 
is a completely different matter. Up to date, it has not been feasible to put density on 
lattice, so there is no model-independent theoretical indication for a similar phase change 
as density is increased. Theoretical models however do predict that there can be a series 
of phase changes including the chiral symmetry restoration. As will be discussed in this 
article, there are a variety of processes involving nuclei and nuclear matter that provide, 
albeit indirect, evidence for such phase changes, most intriguingly of which are possible 
signals from compact (neutron) stars. 

Going to the physics of matter under extreme conditions — high temperature, high 
density or both - encompasses a multitude of scales. At low density p < po (where 
Po is nuclear matter density), one can rely on the rich phenomenology available in the 
guise of standard nuclear physics to work with an effective Lagrangian field theory that 
is anchored on the premise of quantum chromodynamics (QCD), e.g., chiral perturbation 
theory. We may call this "phenomenological perturbation approach (PPA)." However as 
one approaches the density po, an effective field theory constructed in the matter-free 
vacuum encounters difficulty and becomes more or less unpredictive: Presently there is 
no way to derive a nucleus starting from a first-principle Lagrangian, not to mention 
denser hadronic matter. The reason why this is so will be explained in the course of 
this review. The attitude we will adopt here is then to develop a generic description of 
dense matter based on global characteristics of hadrons in the environment of surrounding 
dense medium. For this, we exploit the chiral structure of the QCD vacuum and the 
excitations thereon (i.e., hadrons) following the notion of "vector manifestation" recently 
proposed and developed by Harada and Yamawaki [|T2|. In the Harada-Yamawaki scenario, 
chiral restoration from the spontaneously broken mode to an unbroken mode takes place 
with the massless pions (in the chiral limit) coming together with the would-be scalar 
Goldstone bosons that are the longitudinal components of massive vector mesons, with 
the massless vector mesons consequently decoupling. This scheme requires that near the 
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chiral phase transition — independently of whatever form it may take, the vector {p and 
uo) meson masses drop as the chiral transition point is approached bottom-up. This, we 
will see, provides a compelling theoretical support to BR scaling. We will also see that 
most remarkably, this scenario receives further support from the QCD structure of quark 
condensates in medium associated with color-flavor locking. 

BR scaling as originally formulated in p|| corresponded to a mean-field approxima- 
tion in a large N^. effective chiral Lagrangian theory (where A^^^ is the number of colors in 
QCD) with the strong interaction vacuum Q "sliding" with density or temperature of the 
system. For a given density, say, quantum fluctuations or loop corrections are to be calcu- 
lated for correlation functions etc with a chiral Lagrangian whose parameters are suitably 
defined at the sliding vacuum appropriate to physical processes of given kinematical con- 
ditions. How this works out has been explained in the literature but will be summarized 
in the review at pertinent places. We propose that the way the BR-scaling mean-field 
quantities vary as a function of the medium condition (density or temperature) should 
roughly match with the change of the vacuum suggested by Harada-Yamawaki's vector 
manifestation. As mentioned above, at small external conditions (e.g., at low density), the 
PPA using phenomenological Lagrangians fixed in free space will be in no disagreement 
with the notion of BR scaling. However as density increases, fiuctuations are to be built 
on a vacuum modified from the free-space one. At some point, this effect can no longer 
be accessed adequately by the PPA. This can be seen in what we call "sobar" description 
as explained in a later section. Eventually a correct theory will have to approach Landau 
Fermi-liquid fixed point at normal nuclear matter density. We will show how this comes 
about with BR scaling chiral Lagrangians: BR scaling and a Landau parameter can be 
identified. 

One of the major themes of this review is that this picture combining the vector 
manifestation of chiral symmetry, color-fiavor locking and BR scaling passes several tests 
in nuclear and dense hadronic systems. This, we believe, points to the validity of the 
notion that understanding what goes on in many-body systems of nucleons at high density 
and understanding the chiral structure of the nucleon at the QCD level relies on the same 
principle, a long standing issue in nuclear and hadronic physics. We shall summarize the 
evidences accumulated since our earher publication 0, . 

2 FROM QUARK TO NUCLEON 

2.1 Modeling the Nucleon in QCD 

For concrete applications of QCD to physical hadronic processes, the only model 
that purports to encode QCD and at the same time is versatile enough for a wide range 
of strong interaction physics has been, and still is, the bag model. All other models 
of equal versatility can be considered as belonging to the same class. There is a long 

^Throughout this article, the terminology vacuum will be used in the generalized sense of a ground 
state in the presence of matter. 
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history in the development of this model which, as we shall mention later, continues even 
today in the guise of string theory and large N Yang-Mills gauge theory. The first such 
model that incorporates confinement and asymptotic freedom of QCD is the MIT bag 



model IjT^. However it was realized early on that since this model by construction lacks 



chiral symmetry, it could not be directly applied to the description of nuclear interactions, 
the domain of strong interaction physics most thoroughly and accurately studied. A 
simple stability argument showed that the bag for the nucleon must have a bag size 
typically of ~ 1 fm, essentially the confinement scale implicit in the model. And this is 
simply too big if one naively applies the picture to nuclear systems. As a possible way 
to reconcile this simple bag modelling of QCD with the "size crisis" , it was proposed in 
1979 to implement the spontaneously broken chiral symmetry and bring in pion cloud 
to the model. How to introduce pion cloud to the MIT bag producing a "chiral bag" was 



already known and discussed in the literature [|T^. It consisted of adding into the model 
Goldstone pion fields outside of the bag and imposing chiral-invariant boundary conditions 
that assure both chiral invariance and confinement. The initial idea in ^ was to use the 
pressure generated by the external pion cloud to squeeze the bag to a smaller size so as 
to accommodate the standard nuclear physics picture of meson exchange interactions. At 
the time the idea was proposed, we had no idea how to squash the bag into a "little bag" 
in a way consistent with the premise of QCD. We simply dialed the parameters of the 
model such that the size diminished from ~ 1 fm to ~ 0.3 fm, a size known to be a spatial 
cutoff between two nucleons. We now know that this prescription was not natural and in 
fact totally unnecessary. How to reconcile the bag size and nuclei can now be put in the 
form of what is known as "Cheshire Cat Principle" |TB|p|. 



Apart from the concern with the nucleon size in nuclei and apparent incompatibility 
with the successful standard meson-exchange potential picture, the "little bag" was not 
indispensable for understanding the wide variety of nucleon dynamics as has been argued 



by Thomas and his co-workers P, [18|. In fact, this latter approach where the pion field 
enters only as a fluctuating field at the bag boundary with its size comparable to that 
of the original MIT one has been reasonably successful in a wide range of applications 
to those nucleon and nuclear processes that are dominated by pionic perturbative effects 
and in which nonpertubative chiral effects play a subdominant role 0]. 

It is now well understood how both the "big cloudy bag" of and the "little bag" 
of can be accommodated in a unified framework of the chiral bag model. In fact, 
we now know how to shrink the bag even further from the "little bag" and eventually 
obtain a "point bag," the skymion 1^ which can be considered as the limit that 
the bag size is shrunk to zero. What has transpired from the development since late 
1970's/early 1980's up until now is that the bag size that appears in the simple confinement 
model is more like a gauge degree of freedom [0, ^ and depending upon the gauge 
choice, it can be of any size without changing basic physics. That physics should not 
depend upon the confinement (bag) size is known as "Cheshire cat phenomenon" and 
has been amply reviewed in the literature (for an exhaustive list of references, see [ p!T| , 

historical note: Simultaneously and independently of p^ , a similar idea was proposed based on 
phenomenology by Brown et al [ pT[ . 
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pT| , p2|). Unfortunately there is no known exact formulation of this "Cheshire-Cat gauge 
theory." This is mainly because there is no exact bosonization of QCD known in four 
dimensions, the skyrmion limit corresponding to the totally bosonized theory, and the 
chiral bag partially bosonized. Asking what size of the chiral bag is needed for a particular 
phenomenology is like asking what gauge choice is optimal in gauge theories that are solved 
only partially as is the case with QCD. In practice, what bag size is optimal depends upon 
what physical processes are concerned. 

Thus far our argument has been phrased in the framework of a bag picture. The 
notion of Cheshire Cat is however much more general than in this restricted context. It 
is a reflection of a variety of "duality" in the problem which is a statement that there 
are a variety of different ways of describing the same physical process. In other areas of 
physics (such as string theories, condensed matter physics etc), such a notion can even 
be formulated in an exact way. In the strong interaction dynamics, however, the notion 
is at best approximate and hence its predictive power semi-quantitative. Even so, it is 
sufficiently versatile as to be applied to a variety of problems in nuclear physics. As we 
will describe in later sections, a remarkable recent development in a related context is the 
"complementarity" of two descriptions in terms of hadrons and in terms of quark/gluon 
variables for both dilute matter and dense matter as discussed in a later section of this 
review. This suggests that the Cheshire Cat notion can be generalized to a wider class of 
phenomena, going beyond the "confinement size." 

2.2 The Chiral Bag Lagrangian 

The Cheshire cat phenomenon in our view is a generic feature in modelling QCD or 
other gauge theories. One may therefore formulate it in various different ways. The use 
of the chiral bag is not the only way or maybe not even the best way to exploit it. Since 
we understand the chiral-bag structure better than others, however, we write here the 
explicit form of the Lagrangian that the chiral bag model takes and that should capture 
the essential physics of the nucleon. We shall do this for the flavor SU{3) (with up, down 
and strange quarks) although the nucleon structure is mainly dictated by the non-strange 
quarks, with the strange quark playing a minor and as yet obscure role. It can be written 
generically in three terms, 

S = Sv + Sy + Sav. (|.l) 

The first term defined in the volume V is the "bag" action that refiects the explicit QCD 
variables, quarks ip and gluons G^: 

Sv = d^x {^i - ^tr G^^G^"^ +■■■ (|.2) 

where the trace goes over the color. Here is the quark field with the color, fiavor and 
Poincare indices suppressed, Gf^i, the gluon field tensor and the covariant derivative. 
The ellipsis stands for quark mass terms that we will leave unspecified. The action for the 
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outside sector occupying the volume V contains relevant physical hadronic fields of zero 
baryon charge representing color-singlet effective degrees of freedom. It includes the octet 
Goldstone bosons and non-Goldstone excitations as well as the singlet t]' which acquires 
mass through U{1)a anomaly and takes the form 



^y = ^Iv "^'"^ d.U^d^U + ^mj, [Tr(ln[/ - Inf/t)]^) + . . . + Swzw (|.3) 



P2 

/ d^x (Tr d^U^d^'U+^ 

where Np = 3 is the number of fiavors, Tr goes over the fiavor index and 

U = e'"''//oe2'-//^ (2.4) 



/o = ViV^/2/. 

Here the ellipsis stands for mass terms and heavy-mass fields or higher derivative terms etc. 
The Swzw is the Wess-Zumino-Witten term that encodes chiral anomalies present 
in QCD: 

Swzw = -Nc-r^ I e^^^p^TT^L^UL^LPL'^) (|.5) 

with L = {x, s)dg{x, s) defined as g{x,s = 0) = 1 and g{x,s = 1) = U{x). The 
inside QCD sector and the outside hadronic sector must be connected by an action that 
"translates" them. This is the role of the surface term. It has two terms, 

^dv — ^av + ^av 
with the "normal" boundary term 

S^^ = \f dW{n,i,W-^l^) + ■■■ (|.7) 
/ J av 

with 

^75 _ g«'7'75//0g2i7r75// (@ 8) 

and the "anomalous" boundary term 

S^y^ = ^ /^^ dWK,^{TT ln[/t - Tr \nU) + ■■■ (|.9) 

where is the "Chern-Simons current" 

= e^'^^^iGlGlp - '^-gr'^GlGiGl). (|.10) 

Again irrelevant terms are subsumed in the ellipses. The "normal" term (^.7) assures 
chiral invariance of the model as well as color confinement at the classical level. The 
"anomalous" term - which is not gauge invariant - takes care of the quantum anomaly 
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(Casimir) induced inside the bag that would violate color confinement |2^, |TT| if left 
un-eliminated 0. 

Up to date, nobody has been able to extract the full content of this model. There 
are some unresolved technical difficulties that frustrate its solvability. Even so, whatever 
we have been able to learn from it so far has been found to be fully consistent with Nature. 
We believe that physics of any other viable chiral model, soltitonic or non-solitonic (e.g., 
is essentially - though perhaps not in detail - captured in the model 



Since there have been extensive reviews on the matter , we shall simply summa- 
rize the salient features obtained from the model: 

• Due to a Casimir effect in the bag controlled by the surface term (0.7), the baryon 
charge is fractionized as a function of the bag radius into the V and V sectors 
such that the conserved integer baryon charge is recovered. The baryon charge 
fractionization was first noticed in 1980 by Vento et al. with a reasoning based 



on an analogy to the monopole-isodoublet system studied by Jackiw and Rebbi [EO] 



and established later ||3^ for the special chiral angle ("magic angle") of the pionic 



hedgehog field and then generalized in [^| for arbitrary chiral angle. The result is 
that the confinement size or the bag radius has no physical meaning although the 
baryon size does. The "cloudy bag" model [l^ ignores the hedgehog component of 



the pion field, so in that picture, the baryon charge will be entirely lodged by fiat 
inside the bag. This implies, in the framework of the chiral bag, that the cloudy bag 
is forced to be big to be consistent with Nature. The real question here is whether or 
not such a big bag can capture the essential physics involved. This does not address 
the question of consistency. On the other hand, the limit where the bag radius is 
shrunk to zero with the totality of the baryon charge "leaked" into the "pion cloud" 
describes the skyrmion, with the pion cloud picking up the topological charge. The 
only issue here is whether or not the point-bag description captures the essential 
physics involved. There is no inconsistency in this description of the baryon. 

That the physics should not depend upon the bag radius is an aspect of the "Cheshire 
Cat Principle" (CCP) [l^. It is the trading between the topological character of the 
skyrmion configuration and the explicit fermionic charge of the quarks that makes it 
satisfy exactly the CCP. It is possible to formulate this notion as a gauge symmetry. 
Formulated as a gauge symmetry, picking a particular bag radius R is equivalent to 
fixing the gauge as formulated by Damgaard, Nielsen and Sollacher [^. We will 
call this "Cheshire-cat gauge symmetry." 



■^This form of bulk action plus a surface Chern-Simons term for the bag structure has recently been 
seen to arise from non-perturbative string theory [ p5[ . It turns out that in four dimensions, the interior 
of the hadronic bag can be identified as the 3-brane and the boundary of the bag as the world- volume of 
the Chern-Simons 2-brane in which the dynamics is entirely lodged. This theory does not contain quarks 
as in the model we are considering but the topological term that is added carries the information on the 
anomaly present in the theory. 

We should also point out that the chiral bag structure of the sort we are discussing arises from QCD 
proper in the large Nc limit . 
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While the topological nature of the baryon charge implies an exact Cheshire-Cat, the 
flavor quantum numbers such as isospin, strangeness etc. are not directly linked to 
topology. Therefore their conservation is not directly connected with topology. They 
are quantum variables given in terms of collective coordinates. This means that 
while conserved flavor quantum numbers in QCD remain conserved after collective 
quantization (unlike the classically conserved baryon charge in the skyrmion limit), 
various static quantities such as axial charges qa (both singlet and non-singlet), 
magnetic moments as well as energies and form factors may not necessarily manifest 
the perfect CC property. They must depend more or less on the details of the 
parameters and dynamics taken into account. 

A good illustration of this feature is the A^^^ dependence of the theory (where is 
the number of colors). While the skyrmion sector is strictly valid in the large Nc 
limit, the bag sector containing the QCD variables encompasses all orders of l/N^. 
This means that in order to assure the CC for the quantities that depend on the 
skyrmion sector would have to have appropriate l/iV^ corrections that map to the 
quark-gluon sector. In practice this means that other hadronic (massive) degrees of 
freedom than pionic must be introduced into the outside of the bag. What enters 
will depend on the processes involved. 

Most remarkably, however, extensive studies carried out in the field, chiefly by 
Park [^, confirmed that all nucleon properties, be they static or non-static, topo- 



logical or non-topological, can be satisfactorily and consistently described in terms 
of the Cheshire-Cat picture. Even more surprisingly, the so-called "proton-spin 
problem" can also be resolved in a simple way within the model (@.l) |3^. For this. 



the anomalous surface term (0.9) associated with the U{1)a anomaly and the rj' 
mass figures importantly. As in other Cheshire Cat phenomena (e.g., the baryon 
charge), various radius dependent terms conspire to give a more or less radius in- 
dependent result: Here the cancellation is between the "matter" contribution (from 
both the quark and t]') and the "gauge field" (gluon) contribution, each of which 
is individually bag-radius dependent and grows in magnitude with bag radius. The 
resulting flavor- singlet axial charge for the proton shown in Fig. |l] comes out to 
be small as observed in the EMC/SMC experiments. We should stress that this is 
not the only explanation for the small a°. There are several seemingly alternative 
explanations but our point here is that the Cheshire Cat emerges even in a case 
where QCD aspects are highly non-trivial and subtle involving quantum anomalies 
of the chiral symmetry. 

So far we have been dealing with long-distance phenomena only. What about 
shorter-distance physics? Can Cheshire Cat for instance address such short-distance 
problems as nucleon structure functions seen in deep inelastic lepton scattering 
where asymptotic freedom is operative and hence QCD should be "visible"? The 
answer to this is that as one probes shorter distances, the model (|^.l) should be 
modified with more and more higher- dimension operators implied by the ellipsis in 
the action. There is nothing that suggests that the skyrmion picture when suitably 
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Figure 1: Various contributions to the flavor singlet axial current of the proton as a 
function of bag radius and comparison with the experiment: (a) the "matter" (quark plus 
rf) contribution (a^^ + (b) the gauge field: the static gluon contribution due to quark 
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source {aQ ^j^J and the gluon vacuum contribution 
shaded area corresponds to the range admitted by experiments 
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extended should fail if one works hard enough. It is simply that the calculation will 
get harder and harder and it will be more economical to switch over to the "big 
bag" or partonic (quark/gluon) picture. Logically, therefore, a hybrid description 
that exploits both regimes would have a more predictive power. Thus it is easy 
to understand that the chiral quark-soliton model, a variation on the model (0.1) 



formulated in a relativistically invariant form ||2^, has been successfully applied to 
the problem. We will see later that even at asymptotic density where the QCD 
variables are the correct variables, the notion of effective fields is applicable. 



EFFECTIVE FIELD THEORIES IN NUCLEAR 
PHYSICS 



3.1 Chiral Symmetry in Nuclear Processes 

The next question we wish to address is: Given the Cheshire Cat freedom, how does 
one treat nuclei, namely, two- -or more-nucleon systems? That is, how is nuclear physics 
formulated in QCD? 

For this, it is clear that the most efficient and economical way would be to adopt the 
"point-bag" description 0. Once we adopt the zero-radius bag surrounded by meson cloud 

^In accordance with the Cheshire-Cat principle, one could continue using the hybrid bag-meson de- 
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as a nucleoli, we can represent it as a local field interacting with pions. How to do nuclear 
physics in terms of local nucleon and meson fields and still do QCD is encoded in what is 



called "Weinberg theorem" [g5|. In the present context, the theorem states that doing in a 
consistent way an effective field theory with the nucleon, pion and other meson fields is no 
more and no less than doing the gauge field theory QCD in terms of quarks and gluons. f\ 
What is required is that such an approach preserve the necessary symmetries - Lorentz 
invariance, unitarity, cluster decomposition etc. - and is quantum mechanical. We are 
simply to write down the most general such local Lagrangian with all possible terms 
consistent with the symmetries involved. This strategy is powerfully illustrated in the 
description of the baryon (octet and decaplet baryon) structure in terms of an effective 
chiral Lagrangian constructed of baryon and Goldstone boson fields 0. An initial 



attempt to implement Weinberg theorem in nuclear physics was made in 1981 (it was 
incomplete in the chiral counting at the time but was completed in 1991 as mentioned 
later) but the major development came after Weinberg's paper in 1990 P5|. 

3.2 Objectives of EFTs in Nuclear Physics 

There are two major roles of EFT in nuclear physics. One is to establish that nuclei 
as strongly interacting systems - that have been accurately described in the past in what 
one would call "standard nuclear physics approach (SNPA)" or alternatively "potential 
model (PM)" based on phenomenological potentials - can also be understood from the 
point of view of a fundamental theory, QCD, and the other is to be able to make nontrivial 
and precise predictions that are important not only for nuclear physics per se but also 
for other areas of physics such as astrophysics and condensed matter physics. The former 
renders the impressively rich variety of nuclear processes a respectable domain of research. 
The latter is to provide an invaluable tool for progress in related areas of physics. In this 
review some recent developments that are relevant to the general theme of this review 



scription by incorporating relevant (hadron) degrees of freedom |34 . One advantage of this approach 
might be that the quark-bag introduction can implement more readily short-range physics that enters in 
many-body systems, albeit at a sub-leading level. 

^To quote Weinberg ||3^ : "M^e have come to understand that particles may be described at sufficiently 
low energies by fields appearing in so-called effective quantum field theories, whether or not these particles 
are truly elementary. For instance, even though nucleon and pion fields do not appear in the Standard 
Model, we can calculate the rates for processes involving low-energy pions and nucleons by using an 
effective quantum field theory of pion and nucleon fields rather than of quark and gluon fields. . . . When 
we use a field theory in this way, we are simply invoking the general principles of relativistic quantum 
theories, together with any relevant symmetries; we are not making any assumption about the fundamental 
structure of physics." 

^While the Lagrangian (§.1) is a model of QCD, the baryon chiral Lagrangian used here is an effective 
theory of QCD. It is therefore reasonable to expect that being a theory, the more judiciously organized 
and the more terms computed, the more accurate the calculation will become. There is of course the 
question of convergence but applied within appropriate kinematics, the theory should work out better 
the more one works. This is indeed what is being found in baryon chiral perturbation calculations for 
which the strategy is fairly well formulated. This is not the case with "nuclear chiral perturbation theory" 
where no systematic strategy has yet been fully worked out as we will stress below. 
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will be discussed. 

The examples of remarkable success are the n + p^d + 'ya.t threshold 



classic nuclear physics process - and the inverse j + d —>■ n + p EDI - a process relevant 



to cosmological baryosynthesis. Going beyond the simplest nuclear systems involving two 
nucleons, the method can be applied to processes that involve n nucleons where n > 2, 
such as for instance the solar "hep" process H3] ^He + p — >^ He + e"*" + u^. which 



figures in the solar neutrino problem and the axial charge transitions in heavy nuclei (to 
be discussed later) which provide evidence for BR scaling. 

An exhaustive and beautiful review on the subject with a somewhat different em- 
phasis has recently been given by Beane et al where a comprehensive list of references 
is found. 

In addressing nuclear EFT, there are currently two complementary - and not dif- 
ferent as one might be led to believe - ways of organizing the expansion, one based on 
Weinberg's strategy [03 and the other based on Seattle- Pasadena strategy pS. We will 



loosely refer to the former as "Weinberg scheme" and to the latter as "KSW (Kaplan- 
Savage- Wise) scheme." The Weinberg scheme contains the SNPA/PM (standard nuclear 
approach/potential model) as a legitimate and essential component of the theory whereas 
the KSW scheme by-passes it, although the latter can be effectively exploited to justify 
the former in certain kinematic domains that we are interested in. We will focus on the 
former as we find it more natural and straightforward from nuclear physicists' point of 
view and for the processes we will discuss, sometimes considerably more predictive. 

The Weinberg scheme |]38| , ^ that we will use here is the version that has been 
formulated by Park et al [^, ^ (this will be referred to hereafter as PKMR) and uses the 
following basic strategy. It starts by recognizing that the SNPA/PM (e.g., 0) based 
on solutions of the Schrodinger equation with the "most realistic" potential available in 
the market fine-tuned to a large number of experimental data on one and two nucleon 
systems, captures most accurately the essential physics of, and hence describes, when 
implemented by many-body forces, the bulk of properties of few-nucleon systems, in some 
cases (e.g., low-energy scattering) within ~ 1 % accuracy [|. We will therefore consider the 
SNPA/PM results as input (e.g., "counter terms") in lieu of attempting to compute them 
from first principles. Our aim is not to recalculate them from a more fundamental theory 
but to correctly interpret and incorporate them as a legitimate leading term in a systematic 
expansion of a more fundamental framework [^1], namely, an effective theory of QCD. 



This, we believe, is completely in line with the notion of effective theory, e.g., the one 



^There is a well-known caveat in carrying out this program in general and in assessing the "accuracy" 
of approximations in particular for n-body systems with n > 2. While there are useful constraints that 
help in organizing, in Weinberg's terminology, the "reducible" and "irreducible" terms for two-body 
systems, this is not usually the case for many-body systems. Together with many-body forces in the 
potential which are difficult to pin down precisely although numerically small, there can be a variety of 
organizational ambiguities in classifying "reducible" and "irreducible" contributions, e.g., off-shell effects, 
that could contaminate the calculations resorting to a variety of approximations one is compelled to make, 
even though a full consistent field theoretic approach should in principle be free from such ambiguities. 
This is an issue that will require considerable effort to sort out. We are grateful to Kuniharu Kubodera 
for stressing this point in the contex of the PKMR approach. 
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a la Wilson. Our key point is that it is in how to do embed, and how to systematically 
calculate the quantities that are missing from, the SNPA/PM results into a fundamental 
theory that EFT has its true power. Since the standard nuclear physics approach/potential 
model (SNPA/PM) results are to be inputs to the more fundamental theory, it is clearly 
more profitable to look at the response functions to external fields rather than at scattering 
amplitudes (which the SNPA/PM can address accurately) to make meaningful predictions. 
Our thesis is that this is where chiral symmetry and its broken mode can play their primary 
role in nuclear physics, in both quantitatively describing and accurately predicting nuclear 
phenomena. 



3.3 Effective Chiral Lagrangians 

For very low-energy processes that we are dealing with, the nucleons can be consid- 
ered "heavy," i.e., nonrelativistic while the pions are relativistic. We shall therefore use 
the heavy-baryon formalism although a relativistic formulation can be made, particularly 
for the case where the nucleon mass drops as in BR scaling that we will discuss below. 
The leading-order chiral Lagrangian that consists only of local nucleon and pion fields - 
with other heavy fields integrated out - takes the form 

Co = N[ivD + 2igAS-^]N-]-Y,CA(NTAN 



-f/,^Tr(zA^^A^) + ^Tr(x- 



(1-1) 



with 



1 



2 

1 r 



2 



il.2) 



where TZ^ = ^ (v^ + Al'j and = ^ (v^ — A'^'j denote, respectively, the left and right 
external gauge fields, x is proportional to the quark mass matrix and if we ignore the 
small isospin-symmetry breaking, becomes x = iii the absence of the external scalar 
and pseudo-scalar gauge fields, and 



exp 



r ■ Tc 



(1-3) 



is the chiral field for the Goldstone bosons. The Cq is the leading order chiral Lagrangian 
in the sense that it is leading order in derivatives, in pion mass and in "1/M" where M 
is the free-space nucleon mass. This is a suitable Lagrangian for (chiral) perturbation 
near the medium-free vacuum described by the so-called "irreducible graphs." However 
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for nuclear processes in general, certain "reducible graphs" involve propagators that are 
infra-red enhanced and require terms higher order in 1/M, namely, nucleon kinetic energy 
term. For the KSW approach, the ab initio account of this "correction" term is essential. 
For the PKMR where the reducible graphs are accounted for via Schrodinger equation or 
Lippman-Schwinger equation, the chiral expansion is all that matters, so this Lagrangian 
plays a key role. 

The second term in d^.!) is the leading four-fermion interaction and contains no 
derivatives. We will specify the explicit form later. For convenience, we will work in a 
reference frame in which the four velocity f and the spin operator S'^ are 

= (1, 0) and = (^0, 1^ . (|.4) 

The next-to-leading order Lagrangian including four-fermion contact terms can be written 

as 

= N f ^"^t^ ~ ^'^ D^D, + ^{S ■D,vD} + ciTrx+ + 4 f - ^] {v ■ ^A) 
-f- 4c3^A ■ zA + (^2c4 + ^) \S^, 5^^] [zA„ ^A,]) N 

- 4:idi NS ■ANNN + 2id2 e"''" e^.^xsv^A'^'^NS^^N NS^r^N + ■■■, (|.5) 

where eoi23 = 1, A^ = ^A^. Terms that are not relevant for the present discussion are 
implied in the ellipses. 



3.4 Nucleon-Nucleon Scattering 

Although the genuine predictivity of the EFT is in response functions, we start 
with nucleon-nucleon scattering at very low energies as this issue has been the focus of 
the recent activity in the field. If the energy scale involved is much less than the pion 
mass, ~ 140 MeV, then we might integrate out the pion field and work with an 
effective theory containing only the nucleon field. We will reinstate the pions later since 
the PKMR strategy relies on the pion degree of freedom explicitly. Since the pionless 
theory has interesting aspects on its own that merits study, we will consider this before 
putting the pion field into the picture. 

When the pions are integrated out, the Lagrangian (^.1) further simplifies to 

- ^[CiN^Nf + C\N^(TNy + ■■■] (|.6) 

where the nucleon kinetic energy term (which is formally higher order from the heavy- 
baryon chiral counting) is reinstated. With the pion fields removed from the theory, this 
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Lagrangian has no remnant of chiral symmetry. The nucleon field is just a matter field 
and does not know anything about the symmetry. This does not mean that the theory is 
inconsistent with the symmetry of QCD. We shall consider the S-wave channels ^5*0 and 
^5*1 (in the notation of '^'^'^'^Sj). The C coefficients for these channels are related to the 
C coefficients in (|.6) by 

(jeso) = c' - 3C\ C^'^'^ = C' + CK (|.7) 



x->o(-:>oo< 



+ 



Figure 2: Graphs contributing to the scattering amplitude 
The scattering amplitude is A given by the sum of the Feynman diagrams given in 

Fig.g: 

c 



A 



1 - C{GG) 



where C is the four-fermion interaction constant in (|^.6) and [GG) is the two-nucleon 
propagator connecting the two vertices, both written symbolically. In the center-of-mass 
(CM) frame, this propagator has a linear divergence in dimension D = 4 as one can see 
from a naive dimensional counting of 



{GG) 



(27r)4 \E^qQ- /2M + ie J \-qo - q^/2M + ie ^ 
d^q ( 1 \ 



(27r)3 \E ~qyM + ie^ 



|9) 



where \k\ = a/ ME is the CM momentum. To make the integral meaningful, let us put a 
momentum cutoff at say at A. Then (^9) is given by 

{GGU = -^{^ + ^\k\). (|.10) 

This linear divergence poses no fundamental problem here since the theory is an effective 
one with the degrees of freedom above a certain cutoff having been integrated out. On 
the contrary, in the Wilsonian sense, this cutoff dependence has a physical meaning: It 
delineates the onset of "new physics". Substituting this into the amplitude formula (^8), 
we have 

— C 

^=l + ^(A + .|fc|)- ^-''^ 
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Next, we trade in the theoretical constant C for a physical quantity. To do this, the 
scattering amplitude (for S wave) is written in terms of the phase shift 5 = 5^ as 

47r 1 ,n N 

A=—— — . |.12 

M \k\cot6 -i\k\ ' 

Now using the effective range formula 

|fc|cot5= -- + -ro|fc|2 + C(|fc|^) (|.13) 
a 2 

where a is the scattering length and tq the effective range, we have, to the lowest order 
in momentum. 

So from (|.ll) and (|.14), we find 

where we have explicited the cutoff dependence of the effective Lagrangian which shows 
how an effective theory carries the "memory" of the degrees of freedom that have been 
integrated out. Obviously physical observables should not depend upon what value one 
takes for the cutoff A, since the cutoff can be chosen arbitrarily. This is a statement 
of renormalization group (RG) invariance of the physical quantities, which reads for the 
amplitude as 

A^^O, (S16) 

One can easily verify that (^15) is the solution of this renormalization group equation 
(RGE). 

In standard field theory calculations in particle physics, one uses dimensional reg- 
ularization (DR) instead of cutoff. This is because the DR has the advantage that it 
preserves symmetries, e.g., chiral symmetry in our case, order by order. In the case we 
are concerned with, with no pions, no symmetry is spoiled by the cutoff regularization 
at least up to the order considered. In general when one does higher order calculations, 
however, it is not a straightforward matter to preserve chiral symmetry with the cutoff 
regularization since it requires judicious inclusion of symmetry breaking counter terms so 
as to restore symmetry. This is why some authors (e.g., the KSW collaboration) prefer 
using DR. Now what does the above linear divergence mean in DR? 

For this we go back to the propagator (GG), Eq.(|[9), and work out the integral 
using DR: 



le 
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^n4~d r q 1 



^2' J {2n)D-^ \E - q'^/M + iej 

= Mi-lkl^Y^^-^y^rC--^)-^^ (|.i7) 

This is regular for D ^ 3 but singular for = 3. In renormalizable field theories 
where power divergences are absent in principle, it is legitimate to simply set D = 4, 
thereby dropping power divergences. But in non-renormalizable theories, one cannot 
sweep under the rug the singularities that lurk behind the physical dimension. This 
is because in the Wilsonian renormalization scheme to which our EFT belongs, there 
are no true singularities and as mentioned, the power divergence structure instead carries 
physical information. As mentioned, the situation here is quite analogous to the quadratic 
divergence that needs to be included in hidden local symmetry (HLS) theory discussed 
above 0. The physics behind may be the same. In fact, this must be a generic feature of all 
effective field theories of the sort we are dealing with here, that is, the theories that have 
fixed points nearby. To do a correct regularization with DR in the present case, the linear 
power divergence - the only power divergence there is in this theory at one loop - has to 
be subtracted. In this sense, the "power divergence subtraction (PDS)" scheme used by 
KSW is a proper implementation of Wilsonian effective field theory in nucleon-nucleon 
scattering, the unnatural length scale associated with the large scattering length signaling 
the presence nearby of an infrared fixed point. 

Subtracting the linear divergence, we have the power-divergence-subtracted (PDS) 
propagator 

(GG)r = (GG), - ^^^^^ = -^(/^ + (0-18) 

Comparing with (^.10), we see that the scale parameter fi in DR is equal to the cutoff 
A/tt. 

As noted by ||5^, the large scattering length for the S'-wave scattering, i.e., = 
— 23.714fm ~ (1/(8 MeV)) reflects that nature is close to a fixed point and conformal 
invariance. This is seen in (0.15) in the limit 1/a — > 0: 



Aj^(AC^)|i/a-.o = 0. (0.19) 



In view of the length scale involved, a — oo is close to nature, signaling a quasi-bound 
state. Note that the linear divergence subtraction plays a crucial role here, a feature 
analogous to the conformal invariance of hidden local symmetry (HLS) theory at chiral 
restoration with the quadratic divergence playing a crucial role |12| . 



The discussion made up to this point is general, clarifying the essential feature 
of effective field theories in two-body nuclear systems. The KSW approach has been 

^Taking into account the power divergence (i.e., quadratic divergence) in HLS is crucial in arriving at 
the the vector mani festation of chiral symmetry of Harada and Yamawaki discussed in Section ^ and 



also in Section 11.2 
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developed further, including pions as perturbation |^^. As we will show below, the 
PKMR strategy, while faithful to the EFT strategy, departs from this point in that the 
rigorous counting rule is sacrificed somewhat in favor of predictivity. 



3.5 The PKMR Approach: More Effective EFT 

The faithful adherence to the counting rule has the short-coming in that while post- 
dictions are feasible provided enough experimental information is available, it is hard if 
not impossible to make accurate predictions. The reason is that as one goes to higher 
orders to achieve accuracy, one encounters increasing number of counter terms that remain 
as unknown parameters and these parameters can at best be fixed by the very process 
one would like to predict. This is just a fitting, and not a prediction. Thus instead of 
adhering strictly to the counting rule that requires by-passing wave functions, we shall 
develop a scheme along the line conceived by Weinberg that would allow us to exploit the 
wave functions that are "accurately" determined from the standard nuclear physics/PM 
procedure. 

The key observation that we will exploit is the following. All nuclear processes in- 
volve two sorts of graphs: "reducible" and "irreducible." In the case of two-body systems, 
they will be two-particle reducible and two-particle irreducible (2PI). Instead of sum- 
ming the two classes of graphs order-by-order in strict accordance with a given counting 
scheme as in KSW, one computes irreducible graphs to a given order in chiral perturba- 
tion theory (xPT) and then account for reducible graphs by solving Schrodinger equation 
or Lippman-Schwinger equation with the irreducible vertex injected into the equation. 
This procedure accords to potentials and wave functions a special role as in the standard 
nuclear physics calculations. What we hope to do is then to combine the accuracy of the 
standard nuclear physics/potential model approach with the power of xPT. We might 
call this more effective effective field theory (MEEFT). The advantage of working with 
wavefunctions is that one can make predictions in nuclear electro-weak response functions 
that are difficult to make in an order-by-order calculation because of un-fixed counter 
terms. The disadvantage is of course that one is forced to sacrifice the strict adherence to 
a counting rule. It turns out however that in the processes so far studied with controlled 
approximations, this sacrifice is minor numerically. A clear explanation as to how this 
comes about in two-body systems was given by Cohen and Phillips |^ . 



3.5.1 NN scattering 

To justify the method that we will apply to response functions, we first discuss 
how the method works out in nucleon-nucleon scattering which is well understood in 
various different ways, field theoretical or non-field theoretical. Compared with the KSW 
approach, the present scheme is a bit less elegant but the point of our exercise is to show 



that it can be at least as accurate as the KSW scheme [53 



For simplicity consider proton-neutron scattering in the ^5*0 channel |^. Since we 



^The proton-proton interaction as in the solar proton fusion process requires the Coulomb interaction 
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are to account for the reducible graphs by Lippman-Schwinger equation, we only need 
to compute the potential as the sum of irreducible graphs to a certain chiral order in 
the Weinberg counting. For this we can work directly with the heavy-baryon chiral La- 
grangian (^1) with the pion exchange put on the same footing as the contact four-fermion 
interaction. The potential calculated to the next-to-leading order (NLO) is of the form 



V(q) 



-Ti ■ T2 



9\ 



oi- qa2- q 



4/2 q' + ml 



with 



a 



An 



CTi cr, 



Co + [C^b'' + D2a'^)q'q^ 



+ 



i.20) 



6'^ \ 



1.211 



where q is the momentum transferred. The C2 and D2 terms containing quadratic deriva- 
tives are NLO contributions that are added to (|^.6) since we will work to that order. 
There are one-loop corrections involving two-pion exchange that appear at the NLO and 



that can be taken into account but as shown by Hyun, Park and Min , in the kinematic 
regime we are considering here, they are not important for the key point we want to dis- 
cuss here. As it stands, the potential (^20) is not regularized, so it has to be regularized 
just as the Feynman propagator (GG), Eq.(0.9), had to be. One can do this in various 
different ways; they should all give the same result. A particularly convenient one for 
solving differential equations in coordinate space is to take the following form when one 
goes to coordinate space 



V{r) 



_fq_ 

(27r)3 



1.22) 



with a smooth regulator S^iq"^) with a cutoff A. For our purpose it is convenient to take 
the Gaussian regulator, 

SAiq') = exp (^-^^ . i.23) 

In the PKMR formalism, pions play a special role due to what is called "chiral 
filter mechanism" (to be defined below) but we can illustrate our essential point without 
the pions, so let us drop the pion exchange term for the moment. Given the Fourier- 
transformed potential properly regulated as mentioned, the Lippman-Schwinger equation 
for the S-wave can be solved in a standard way. The resulting wave function can be 
explicitly written down: 



A2 } 



c, 



1 - TeCe 



ZC2 
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(V^ + ME) 



ZD2Suif) d 1 d 



Ce VS dr r dr 



fA(r) 



1.24) 



in addition to the QCD interaction. Although technicaUy dehcate, this part is known, so we will not go 
into the detail. We shall also restrict to the ^5*0 channel. The "^Si channel is slightly complicated because 
of the coupling to the D-wave but involves no new ingredient. 
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where is the free wave function and 
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i-26) 
(1-27) 
(1-28) 

(3.29) 
(1-30) 



where J„ (n = 2, 4) are defined by 
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The phase shift can be calculated by looking at the large-r behavior of the wave- 
function, 



pcot 5 
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where the rj^E) is the D/S ratio, which vanishes for the ^Sq channel and Ia{E) = F^; + 
i\/ M E S"^ {M E / Is?) = AI{^j§-). The expression for the phase shift contains a lot more 
than what was considered above even at LO since the wavefunction is computed to all 
orders in the potential and we are using a Gaussian cutoff. As such it is not obvious 
that the renormalization group invariance, i.e., the A independence, is preserved. This 
may appear to be a shortcoming of the approach and would constitute a defect if the A 
dependence were non-negligible. 

Let us see how this works out with the present pionless example. As before, we use 
the effective range formula (|[13) to fix the coefficients Co,2- We have 
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At threshold, it can be easily verified that the renormalization-group (RG) invariant form 
for Co, (1^.15), is recovered. What happens away from the threshold is less trivial. For 
illustration we take the CM momentum p = 68.5 MeV corresponding to ~ m^/2. The 
result for the ^Sq channel is given in Fig.^; here the phase shift 6 (in degree) is plotted vs. 
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Figure 3: Cutoff dependence/independence in pionless theory: np ^Sq phase shift (degrees) vs. 
the cut-off A for a fixed CM momentum p = 68.5 MeV. The solid curve represents the NLO 
result, the dotted curve the LO result and the horizontal dashed line represents the experimental 
value. 



A. We learn from this exercise that within the present scheme it is clearly inconsistent 
with the renormalization group invariance if the potential is calculated only up to LO. It 
is at NLO and A > m.„ that the cutoff independence is recovered. We see that A ~ is 
the scale at which a "new degree of freedom" enters into the pionless theory, an obvious 
but important fact. 

Next we incorporate pions. Once pions are included, the cutoff should be put above 
m„. Furthermore the pion presence should reduce the A dependence - if any - for given 
CM momenta. This feature is seen in Fig.^. We see clearly that the pion presence 
improves markedly the ability to describe scattering. The interplay between the probe 
momentum, the pion presence, the cutoff and the chiral order becomes transparent: The 
more refined the potential and the higher the chiral order, the more consistent and more 
accurate becomes the prediction. This will be the key point of our next consideration 
when we look at response functions. 



3.5.2 Electroweak response functions: Accurate post- dictions and predictions 



We now turn to our principal assertion that it is in calculating nuclear response functions 
to the electroweak (EW) external fields that the PKMR approach has its power. We 
will first discuss two-body systems and then make predictions for n-body systems with 
n > 2. In doing this, we will put pions to start with since we will exploit what is known 



as "chiral filter mechanism" first introduced in [51| and elaborated further in |48]. The 



chiral filter mechanism which has been verified in xPT states that whenever allowed by 
symmetry and kinematics, one-soft-pion exchange in the electroweak currents dominates, 
with the corrections to which are in principle controllable by xPT and, conversely, when- 
ever forbidden by symmetry or suppressed by kinematics, the corrections to the leading 
order current are uncontrolled by chiral expansion and require going beyond standard 
low-order xPT. We will refer to these two phenomena as "two sides of the same coin." In 
this formalism, therefore, pions play a prominent role. At a momentum scale p <C m^r. 
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Figure 4: The np ^Sq phase shift (degrees) vs. the cutoff A for fixed CM momenta p = 70 MeV 
(left) and p = 140 MeV (right). The next-to-leading order (NLO) results are given by the solid 
(with pious) and dashed (without pious) curves, and the leading-order (LO) results by the dot- 
dashed (with piou) and dot-dot-dashed (without pion) curves. The horizontal line represents the 
experimental data obtained from the Nijmegen multi-energy analysis. The NLO result without 
pions is drawn only up to A > 300 MeV, above which the theory becomes meaningless, that is, 
unnatural as well as inconsistent. 



one can work with a pionless EFT as described in M in which the role of the chiral filter 
would of course be moot. However in this case, one loses the predictive power since one 
would be left with one or more counter terms that are left undetermined by theory or 
experiments. The PKMR strategy circumvents this difficulty. 

The PKMR strategy is quite simple. We take the most sophisticated wavefunctions 
from the SNP (standard nuclear physics)/PM (potential model) approach with the po- 
tential fit to an ensemble of empirical data on two-body systems and compute the EW 
currents using the chiral Lagrangian (^1) to as high an order as possible in the chiral 
counting. Such realistic potentials - which have been extensively studied and confirmed 
empirically - can then be interpreted as resulting from high-order xPT calculations. For 
two-body systems, we use as the input the most "sophisticated" potential available in 
the literature, i.e., the Argonne vl8 potential - and the resulting wavefunctions |^0|- For 



making predictions in n-hodj systems with n > 2, we rely on Ref.[49|. The power of these 



wavefunctions in confronting few-body nuclei is well documented in the literature. In all 
cases, the scattering and static properties of the nuclei involved are extremely accurately 
given and more or less justified within xPT whenever many-body corrections are unim- 
portant. Whenever many-body corrections are important, on the other hand, they can in 
turn be calculated by xPT as we shall show using the chiral filter argument. In all cases 
studied up to date, it is possible to make error estimates by assessing the deviation from 
the RG invariance, that is, the A dependence. 

Before calculating the full amplitudes to confront nature, we first verify that using 
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the "realistic wavefunction" is indeed consistent with actually calculating the wavefunc- 
tions systematically in the chiral counting starting with (^.1) or (§.6). For this purpose, 
it suffices to look at the single-particle Ml matrix elements for the process 

n + p^d + j (|.35) 

and the single-particle Gamow- Teller matrix element in the solar proton fusion process 

p + p^d + e-^ + Ue (|.36) 

and verify that there is a matching between a bona-fide EFT and the hybrid scheme we 
are using. We are interested Q in £^(M1) and £^(GT) defined by 



£(M1) ^ .r"' , f(GT) s :^^S£-#ST, (337) 

where A^mi -^mi denote, respectively, the single-particle Ml transition matrix el- 
ement of our NLO calculation considered above for the np scattering and that of the 
Argonne potential, and similarly for £^(GT). Here we are taking the Argonne potential as 
the most accurate and closest to "experiment." In Fig.^ are plotted the deviations from 
the Argonne result as a function of the cutoff. We see that as in the case of the phase 
shift for p ~ m^, the inclusion of the pion degree of freedom markedly reduces the cutoff 
dependence, in conformity with the stated requirement for a viable effective theory. 

Leaving the details to the literature ||4^, we merely summarize a few concrete 
results, both post-dictions and predictions, of full calculations to NLO or NNLO including 
pions. Some technical details including precise definitions of the quantities involved can be 
found in the next subsection immediately following this where we address the question of 
"hard-core correlations" invoked extensively in standard calculations dealing with short- 
distance nuclear interactions. 

• POST-DICTION 

The most celebrated case of post-diction is the thermal np capture process (0.35). 
The unpolarized cross section is dominated by the Ml operator - defined below 
- that is protected by the "chiral ffiter" and has been calculated in the PKMR 
approach with the Argonne f 18 to the accuracy of ~ 1% P3, 



(7th = 334 ± 3 mb (|.38) 

to be compared with the experimental value aexp = 334.2 ± 0.5 mb. The error bar 
on the theoretical value represents the range of uncertainty in the short- distance 
physics unaccounted for in the NNLO chiral perturbation series that is reflected 
in the cutoff dependence. Here the correction to the single-particle Ml matrix 
element is dominated by the soft-pion contribution with the remaining corrections 



^"^These are of course not the whole story since there are two-body corrections known as 
exchange currents" that we wiU consider shortly. 
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Figure 5: £{M1) (upper) and £{GT) (lower) vs. the cutoff A. The soHd curves represent the 
NLO results with pious and the dotted curves without pious. When pious are included, the 
physically relevant cutoff must lie between the pion mass m.,^ and the lowest resonance, say, 
rria ~ 450 MeV. 



amounting to less than 10% of the unsuppressed soft-pion term in conformity with 
the chiral filter argument. Since this calculation has no free parameters, it is actually 
a prediction whereas the KSW calculation without pions has one unknown counter 
term as the leading correction term |^. The upshot of the xPT result a la PKMR 
is that it confirms the result of Riska and Brown ^ and in addition supplies the 
corrections to their result with an error estimate, which is the ultimate power of the 
approach as far as post-dictions are concerned. Now having explained accurately 
the threshold up capture, one can then proceed to make predictions for the inverse 
process c? + 7 — > n + p for which cross section data are available as well as for the 
process e+d n+p+e as a function of momentum transfer. A good agreement with 
experiments is obtained for the photodisintegration of the deuteron d+'j — n+p. 



Remarkably the simple theory seems to work also for the latter to a momentum 
transfer much greater than justified by the cutoff involved 0. 

PREDICTIONS 

To bring our point home, we cite a few cases of predictions that precede experi- 
ments. For these predictions, the xPT Lagrangian (0.1) is used for calculating the 



^^In order to explain the experimental value in this calculation, the counter term turns out to require to 
be considerably bigger in magnitude than, and to have a sign opposite to that of, the soft-pion correction 
of ||39[| . This seems to imply that the counter term that appears in this pionless theory describes a different 
physics than the chiral-filter protected mechanism. 
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irreducible graphs entering in the electroweak currents. 



1. While the unpolarized cross section for the np capture process is dominated 
by the isovector Ml operator that is protected by the chiral filter mechanism, 
the process 

n + p^d + -f (g.39) 

with polarized neutron and proton can provide information of the suppressed 
isoscalar Ml and E2 matrix elements. The isoscalar Ml matrix element turns 
out to be particularly interesting because a precise calculation of this matrix 
element requires implementing "hard-core correlations" of standard nuclear 
physics in terms of regularization in field theory. Although the understanding 
is still poor, the issue presents a first indication of how the short- distance 
prescription encoded in the hard core may be understood in a more systematic 
way. 

As noted more precisely below, the matrix elements involved are suppressed 
by three orders of magnitude with respect to the isovector Ml matrix element 
- due to the fact that these are unprotected by the chiral filter mechanism, 
so cannot be "seen" directly in the cross section but can be singled out by 
polarization and anisotropy measurements (see below for definitions). The 
ratios with respect to the isovector Ml are predicted to be 

TZmi = (-0.49 ± 0.01) X 10-^ 7^i^2 = (0.24 ± 0.01) x 10"^ (gAO) 

There are no data available at the moment; in this sense these are true pre- 
dictions. The error bars represent the possible uncertainty in the cutoff de- 
pendence, indicating the uncertainty in the short-distance physics. We should 
mention that the identical results were obtained in the KSW scheme, with 
however, somewhat bigger error estimates M . 



The next prediction is on the proton fusion process (^36). Here the domi- 
nant matrix element involved, i.e., Gamow- Teller, is unprotected by the chiral 
filter but it turns out that the single-particle matrix element which can be 
precisely calculated in the standard nuclear physics/PM approach with accu- 
rate wave two-nucleon wave functions dominates. The remaining correction 
terms are not protected by the chiral-filter mechanism and hence can be sub- 
ject to some uncertainty although numerically small ||5^ . It turns out however 
that to next-to-next-to-leading order (NNLO) it is possible to make an ex- 
tremely accurate totally parameter-free estimate of the correction [Q. The 
quantity that carries information on nuclear dynamics is the astrophysical S 
factor, Spp = a{E)EE^^'' where cr is the cross section for the process (0.36) 



and 1] = Ma/{2p). The prediction for 5* is p8 



S'pp(O) = (3.94 X 10^2^) (1 ± 0.0015 ± 0.001) MeV barn (|.41) 
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where the first error is due to the uncertainty in the one-body matrix element 
and the second due to the uncertainty in the two-body matrix element that is 
truncated at the NNLO. Higher-order corrections are expected to be completely 
negligible. This probably is the most accurate prediction anchored on the 
"fundamental principle" achieved in nuclear physics. 

In the solar neutrino problem, the S factor is an input as a Standard Model 
property. In view of the important role that the pp fusion process will play in 
the future measurements of the solar neutrinos, this accurate Standard Model 
prediction will prove valuable for the issue. 

So far we have been concerned with two-nucleon systems that can be more 
or less equally well treated by both the Kaplan-Savage- Wise (KSW) approach 
and the Park-Kubodera-Min-Rho (PKMR) approach. We now turn to a case 
where the PKMR can make a possible prediction ||^, ^ whereas a strict ad- 
herence to the counting rule would frustrate the feat due to too many unknown 
parameters. This is the solar "hep" problem 

p +3 He He + e+ + u^. (|.42) 

This process produces neutrinos of the maximum energy E,^ ^ 19.8 MeV; for 
this lepton momentum, not only the S wave but also the P wave need to be 
taken into account. Computing the neutrino flux as needed for the solar neu- 
trino problem requires both vector and axial-vector matrix elements. While 
the vector matrix element involves no subtleties - and hence is straightfor- 
wardly calculated (including two-body currents), the principal matrix element 
of the axial-vector current turns out to be extremely intricate. Since the P 
wave enters non-negligibly, the time component of the axial current cannot be 
ignored. However this part of the current is under theoretical control since the 
axial charge operator is chiral-filter protected |47[| : the corrections to the 



known single-particle axial-charge matrix element are dominated by one-soft- 
pion-exchange two-body matrix elements (with three-body operators absent 
to the order considered) and are accurately calculable, given the "accurate" 
wave functions available in the literature. The principal difficulty is in the 
matrix element of the space component of the axial current, the most impor- 
tant of which is the Gamow- Teller operator. (Three-body contributions are 
higher-order than those coming from one-body and two-body terms and can 
be ignored). Normally the matrix element of the single-particle Gamow- Teller 
operator qa J2i T^{i)<y°'{i) is of order 1 if the symmetries of the initial and final 
nuclear states are normal since the operator just fiips spin and isospin. How- 
ever in the "hep" case, the dominant components of the initial state |p^He) and 
the final state |^He) have different spatial symmetries so the overlap is largely 
suppressed. One finds, numerically, that to the chiral order that is free of un- 
known parameters (e.g., next-to-next-to leading order, NNLO, or 0{Q^)), the 
single-particle GT matrix element comes out to be suppressed by a factor of ~3 
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relative to the normal matrix element. Now main corrections to the one-body 
Gamow- Teller are expected to be from two-body terms with three-body terms 
down by a chiral order. However these corrections are chiral-filter-unprotected 
as we have explained and hence cannot in general be calculated accurately 
in chiral perturbation theory unlike the corrections to the single-particle axial 
charge matrix elements. Furthermore, the two-body corrections typically come 
with a sign opposite to the sing-particle one with a size comparable to the main 
term, thus causing a drastic further suppression. This shows why it is essential 
to be able to calculate in a controlled fashion higher chiral order corrections to 
the GT matrix element. 

Such a calculation turns out to be feasible given accurate wave functions for 



three- and four-nucleon systems ||41|. The key point is that the calculation 
of this chiral-filter-unprotected term parallels closely the calculation of the 
suppressed isoscalar Ml and E2 matrix elements in the polarized np capture 
mentioned above, thereby enabling the PKMR approach to pin down the most 
problematic part of the "hep" process within a narrow range of uncertainty 



dictated by the range of the cutoff involved Referring for details to the 



literature we simply summarize what enters in the calculation. The 

first is that as in the case of the pp case, there are no unknown parameters 
once the accurate data on triton beta decay are used. The second is that the 
strategy for controlling the chiral-filter-unprotected two-body corrections to 
the Gamow- Teller matrix element parallels that for the pp fusion process, i.e., 
the renormalization-cutoff interpretation of what is known in nuclear physics 
as "hard-core correlation" discussed below. 

The S factor so obtained |^ 



io obtained |^3| is 
Shep{0) = (8.6 X 10-2°)(1 ± 0.15) MeV barn. 



(0-43) 



The error bar stands for uncertainty due to the cutoff dependence which sig- 
nals the uncertainty inherent in terminating the chiral series at NNLO. The 
uncertainty is bigger here than in the pp case since one is dealing with four- 
nucleon systems. Higher-order corrections are expected to considerably reduce 
the error bar. 

A similar calculation in a standard nuclear physics approach based on the tra- 
ditional treatment of exchange currents |Q has been performed by Marcucci et 
al |^T|. This approach uses the same accurate wave functions as in [Q, without 



however resorting to a systematic chiral counting. That the standard nuclear 



physics approach of agrees with (|^.43) within the error bar suggests that 



the standard approach of [41] may well be consistent - to the order considered 



- with the EFT approach to the "short-range correlation. 
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3.5.3 The chiral filter: Two sides of the same coin 

While the processes protected by the chiral filter, such as the Ml operator and the 
axial charge operator are accurately calculable since one-soft-pion exchange domi- 
nates with correction terms suppressed typically by an order of magnitude relative to the 
dominant soft-pion exchange, the situation is not the same for those processes that have 
no chiral-filter protection. Even so, the PKMR approach can still make a meaningful 
statement on such processes. That is the story of the other side of the same coin. 

If the chiral filter does not apply, then the role of the pion is considerably diminished 
and, compounded with higher order effects, clean and reliable low-order calculations will 
no longer be feasible. A prime example of this situation is the process p + p ^ p + p + 7r^. 
Here one-soft-pion exchange is suppressed by symmetry and a variety of correction terms 
of comparable importance compete in such a way that no simple xPT description can be 
obtained. Both the PKMR and the KSW are unsuccessful for this process. 

The proton fusion process (|^.36) presents a different situation at least in the PKMR 
approach |^. Since the process is essentially dominated by the Gamow- Teller operator 
with no accidental suppression, the single particle matrix element in the standard nuclear 
physics (SNP) /potential model (PM) approach simply dominates. On the other hand, 
since one is dealing with low-energy (or momentum) probes with the relevant scale much 
less than the pion mass, one may simply integrate out the pion and work with a pionless 



effective Lagrangian as in ||5^. However the price to pay here is to face higher-order 
counter terms. We learn from the chiral filter that these corrections cannot be computed 
in a low-order expansion. This is because those terms reflect the physics of a distance 
scale shorter than that that can be captured in a few perturbative expansion. While in 
the PKMR approach, the counter terms of high order that are required are numerically 
small compared with the leading single-particle term given in the SNP/PM approach, 
with no help of the chiral filter, there is nothing that says that subleading corrections 
are necessarily small compared with the leading-order correction that one may be able 
to calculate. For example in high-order KSW calculations in the "pionless EFT" for the 



proton fusion (see |^), the corrections occur at fifth order. Now even if the counter terms 



that appear at that order are - eventually - fixed by the inverse neutrino processes that 
may be measurable in the laboratories, the error incurred in this fifth order calculation 
will not be any smaller than what the PKMR approach can obtain at a lower order. 

3.5.4 An interpretation of hard-core correlation in an effective field theory 

The isoscalar matrix elements of the polarized neutron and proton capture process 
(1-39), 

n + p^ d + -f (0.44) 

are not chiral- filter-protected. Nonetheless it was possible to make the predictions (0.40). 
How this comes about is interesting from the point of view of learning something about 
short-range correlations in nuclei. We elaborate on this point in the rest of this section. 
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At threshold, the initial nuclear state in (0.44) is in either the ^So{T = 1) or 
^5*1 (T = 0) channel, where T is the isospin. The process therefore receives contributions 
from the isovector Ml matrix element (MIV) between the initial ^5*0 (T = 1) and the 
final deuteron (T = 0) state, along with the isoscalar Ml matrix element (MIS) and the 
isoscalar E2 (E2S) matrix element between the initial ^5*1 (T = 0) and the final deuteron 
^Si—^Di states. Since MIV is by far the largest amplitude, the spin-averaged cross section 
<^unpoi{iT-p — > d'-f) is totally dominated by MIV. Meanwhile, since the initial ^5*0 state has 
J = 0, the MIV cannot yield spin-dependent effects, whereas MIS and E2S can. 

The T matrix for the process can be written as 

{^,,^{k,^)mnp)=X^dMCk,X)Xnp (1-45) 

with 



(1-46) 



M{k,X) = VA7T-^^ ^(fcxe*)- (0-1-0-2) MIV 

■ ft -*\ / , MIS - - „ E2S 

-i[k X Ey) ■ (cTi + era) — ^ + (cri • k(T2 • + 0-3 ■ kcTi ■ ^^x)-^ 



where f„ is the velocity of the projectile neutron, Ag is the deuteron normalization factor 
Ag — 0.8850 fm"^''^, and Xd {Xnp) denotes the spin wave function of the final deuteron 
(initial np) state. The emitted photon is characterized by the unit momentum vector k, 
the energy oj and the helicity A, and its polarization vector is denoted by eA = e\(k). The 
amplitudes, MIV, MIS and E2S, represent the isovector Ml, isoscalar Ml and isoscalar 
E2 contributions, respectively, all of which are real at threshold. These quantities are 
defined in such a manner that they all have the dimension of length, and the cross section 
for the unpolarized np system takes the form 

a^npoi = |M1V|' + |M1S|' + |E2S|' . (|.47) 

The isoscalar terms (|M1S|^ and |E2S|^) are strongly suppressed relative to |M1V|^ — 
approximately by a factor of ~ O(10~^) — so the unpolarized cross section is practically 
unaffected by the isoscalar terms. As mentioned above, the isovector Ml amplitude was 
calculated very accurately up to 0{Q^) relative to the single-particle operator. The 
result expressed in term s of M IV is: MIV = 5.78 ± 0.03 fm, which should be compared 
to the empirical value \J(T^npoi — 5.781 ± 0.004 fm. Here, we will focus on the isoscalar 
amplitudes. 

Now the isoscalar matrix elements can be isolated by spin-dependent observables, 

J j-qO \ J /qO \ 

namely by the photon circular polarization = and the anisotrophy 77^ = 

j(90°)+j(o°) where I\{0) is the angular distribution of photons with helicity A = ±1, with 

9 the angle between k (direction of photon emission) and a quantization axis of nucleon 
polarization |^. Measurement of P^ and rj^ can determine the empirical values of the 

""^^Explicit forms are not needed for our discussion but we write them down for completeness: 

, ^ , V2(7^Ml - 7^E2) + |(7^Ml + 7^E2)' 

-L ^ '^Ml ^ '^E2 
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ratios 



MIS 
MIV 



7^E2 = 



E2S 
MIV 



(|.4J 



Both MIS and E2S are suppressed relative to the isovector Ml matrix element 
(MIV) due to mismatch in symmetries of the wave functions etc. We will however take 
the single-particle matrix elements of MIS and E2S as the leading order (LO) in chiral 
counting of the higher order terms in calculating the ratios T^mi and 7^e2- Radiative 
corrections to these leading-order terms come at N^LO and are calculated easily given 
the wave functions 0. Relative to the leading order (LO) terms, corrections will then be 
classified by N'^NO for 0{Q"'). For the chiral-filter protected MIV, the leading two-body 
correction appears at N^LO and the one-loops and the counter-terms appear at N^LO. 
For the isoscalar current, not protected by the chiral filter, the counting rules are quite 
different. While the E2S receives only negligible contributions from higher order (N^LO 
and N^LO) corrections, the situation with MIS is quite different. The corrections are 
quite large even though the leading two-body corrections for the MIS turn out to appear 
in tree order at N^LO while the loop corrections come at N^LO and are finite. 




Figure 6: Generalized tree diagrams for the two-body isoscalar current. The solid circles include 
counter-term insertions and (one-particle irreducible) loop corrections. The wiggly line stands 
for the external field (current) and the dashed line the pion. 



Given that the higher corrections separately are quite large as expected from the 
chiral-filter argument, there is a problem which does not arise in chiral-filter-protected 
processes. It has to do with the zero-range interactions that come as radiative corrections 



Vf = pP 



7^2 



^E2 ^ 67?.Ml'7^E2 



4(l-pP) + (4 + pP)(7^2I^ 



7^|2) + 2pP7^Ml7^E2 



where pP = Pp-Pn and P is the polarization three- vector. 

^•^ Although not in the main line of discussion, we quote as an aside the numbers to give an idea of the 
size we are dealing with. With the realistic Avl8 wavef unctions, the numerical results for the sum of the 
LO and NLO come out to be 

MlSiB(fm) = (-4.192- 0.105) X 10-3 = -4.297 X 10^3^ 
E2SiB(fm) = (1.401 - 0.007) x 10"^ = 1.394 x 10"^ 
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Figure 7: One- loop graphs contributing to the two-body currents. They come at ©(Q^) and 
higher orders relative to the LO one-body term. All possible insertions on the external line are 
understood. 



involving the four- Fermi contact interactions (proportional to Ca) in Eq.(||.l) (see Fig.|^) 
and as four- Fermi counter terms (see Fig.^). It turns out that there is only one term of 



the latter type that contributes, i.e., a term in the Lagrangian of the form 



|.49) 



where is the isoscalar external field and (74 is an undetermined coefficient. When the 
chiral-filter mechanism is operative, these zero-range interactions appear at higher order 
and are expected to be suppressed by "naturalness" conditions and more importantly 
by "hard-core" correlation functions that are incorporated in the wave functions. In the 
present case, the chiral filter is not effective and hence it is not a good approximation 
to let hard-core correlation functions "kill" the contact interactions. One way to resolve 



this problem is to exploit the cutoff regularization that was used above in Section |3.4| for 
NN scattering. Operationally, we can achieve this goal by replacing the delta function in 
coordinate space attached to zero-range terms by the delta-shell with hard core r^. 



47rr^ 



(|.50) 



We then have to assure that the total contribution we compute be independent of 
within a reasonable range. This requirement constitutes a sort of renormalization-group 
invariance. When the chiral filter mechanism is operative, dropping the zero-range terms 
is justified since the independence is not spoiled. If this requirement is not satisfied, 
the effective theory is not reliable and the result cannot be trusted. The present case 
belongs to this class. 

It turns out fortunately that for MIS, there enters only one particular linear com- 
bination of four- Fermi terms involving and Ca's and that it is this combination - call 
it X - that figures in the magnetic moment of the deuteron. Fitting to the deuteron 
magnetic moment allows the constant X to be uniquely determined for any given value 
of Tc. The result for the total MIS for wide-ranging values of Vc is given for illustration in 
Table |l|. The remarkable insensitivity to may be taken as a sign of RG invariance and 
a support for the procedure. This is more remarkable considering that the series in the 



32 



rc (fm) 


0.01 


0.2 


0.4 


0.6 


0.8 


MIS 


-2.849 


-2.850 


-2.852 


-2.856 


-2.861 



Table 1: The total contributions to MIS in unit of 10 ^fmvs. Tc in fm. Recall the single-particle 
contribution: MlSis -4.30 x lO'^ fm. 



chiral order does not appear to be converging. To exhibit this, we quote how each term 
of N^LO for n > 3 contributes to the ratio TZmi in the range of 

r™" = 0.01 fm < re < r^^ = 0.8 fm. (|.51) 

Expressed in the chiral order and in the range (^t-^^x^ ^™n-j^ ^j^^ TZui comes out to be 

7^A/l X 10^ = "LO" + "N^'LO" + "N^LO" + "(N^LO + N^LO)x" 
= -0.74 + (-0.48, -0.74) + (0.23, 0.46) + (0.49, 0.54) 
= (-0.50,-0.49) (|.52) 

where the contribution that depends on the single parameter X fixed by the deuteron 
magnetic moment is indicated by the subscript X. This is how the value given in (0.40) 
has been arrived at. The crucial observation which applies directly also to the hep process 
- and to a somewhat less extent, to the hep process - discussed above is that the final 
result cannot be arrived at by any partial sum of the terms. Also note the importance 
of the X-dependent term. Needless to say, an experimental check of the result obtained 
here would be highly desirable. 

It should be noted that the scale we are talking about in connection with the "hard- 
core regularization" is of the chiral scale ~ 47r/^ that delineates low-energy degrees 
of freedom of hadrons from high-energy, and not that of nuclei that enters in low-energy 
two-body interactions discussed above. 

4 "VECTOR MANIFESTATION" OF CHIRAL 
SYMMETRY 

4.1 Harada-Yamawaki Scenario 

The major problem that we have to face in going to many-body systems is how the 
structure of the ground state, i.e., "vacuum" is to be described and how the properties of 
hadronic degrees of freedom get modified in the changed background. This question will 
be addressed more specifically in the next section in terms of an effective in-medium field 
theory. Here we address the generic aspects based on symmetries involved in the problem. 
We will discuss the problem in terms of density but we believe the general argument applies 
as well to temperature at least up to the chiral restoration point. Possible differences that 
appear after the critical point will be mentioned in Section 11. 
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To discuss how the vacuum changes as matter density increases and how the prop- 
erties of baryonic matter are affected by the vacuum change, we need to first address how 
chiral symmetry is affected by the medium. While it is generally accepted that chiral 
symmetry endowed with iVp-flavor massless quarks breaks spontaneously in the (zero- 
density) vacuum as SU{Nf)l x SU{Nf)r SU{NF)L+ji, it is not fully known how the 
broken sjTumetry gets restored in medium under extreme conditions, that is, at high 
temperature and/or at high density. There are several different scenarios with which the 
broken symmetry can get restored as temperature and/or density are/is dialed. As we 
will argue, what happens in nuclei and nuclear matter where matter density is non-zero 
depends rather crucially on what scenario the chiral symmetry restoration adopts. 

For concreteness and simplicity, take the two-flavor case with up and down quarks 
{Np = 2) 0. The "standard scenario" is the a model one where the broken Nambu- 
Goldstone mode gets restored to the Wigner-Weyl mode wherein the triplet of pions vr* 
with i = ±,0 and the scalar a (in the chiral limit) merge to an 0(4) massless degenerate 
multiplet at the phase transition. At the phase transition, other hadrons such as nucleons, 
vector mesons etc. could be either massless or massive as long as they appear in chiral 
multiplet structure. Now in zero density, the scalar a need not be the fourth component 
of 0(4) as in linear a model. It could for instance be a chiral singlet but what is required 
is that at the phase transition it joins the pions into the 0(4) as described by Weinberg 



in his "mended symmetry" scenario ||6y, |61|. One should note that the cr-model scenario 



is moot on the spin-1 mesons, so, e.g., the vector meson p (which will figure importantly 
in later discussions) could become either lighter or heavier as one approaches the chiral 
restoration. This model belongs to the 0(4) universality class and has been extensively 
discussed in the literature At present, it is not inconsistent with lattice measurements 



with temperature but it is not the only viable scheme either. 

There is an alternative scenario to the above standard one recently proposed by 
Harada and Yamawaki |jl2| which we believe provides support to the BR scaling to be 
elaborated below as well as connects with the color-flavor locking scenario of QCD dis- 
cussed in Section |11.2| . The distinguishing feature of Harada-Yamawaki's "vector man- 
ifestation" is that chiral symmetry with, say, two flavors (u and d) is restored by the 
triplet of pions merging with a triplet of the longitudinal components of vector mesons to 
the representation (3, 1) © (1, 3). At the chiral point, the vector coupling gy flows to the 
fixed point gy = 0,a = 1 corresponding to what Georgi called vector symmetry limit 
together with the pion decay constant /,r going to zero. In this case, the vector p will 
become massless together with other matter fields, i.e., constituent quarks and the vector 
coupling will vanish, a feature which is an important factor in our discussions to follow. 

Limiting ourselves to the two-flavor case (the three-flavor case can be similarly 
treated), we consider specifically the hidden local symmetry (HLS) theory of Bando et 
al with the symmetry group [U{2)l x U{2)ii]giobai x [U{2)v]iocai consisting of a triplet 
of pions, a triplet of p-mesons and an cj-meson. Motivated by the observation that in the 

-^"^The Harada- Yamawaki argument was originally developed with three flavors u, d and s. Although 
it has not been verified explicitly, we see no reason why the same argument should not apply to the 
two-flavor case with some minor changes. 
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vacuum, the p and u mesons are nearly degenerate and the quartet symmetry is fairly 
good phenomenologically, we put them into a U{2) multiplet. In this theory, baryons 
(proton and neutron) do not appear explicitly. They can be considered as having been 
integrated out. If needed, they can be recovered as solitons (skyrmions) of the theory. 

The relevant degrees of freedom in the HLS theory are the left and right chiral fields 
denoted by iifi, 

with cr(x) = ij°'T°' and tt = 7r"T" and the hidden local gauge fields denoted by 

V, ^ V^T- = + (|.2) 

with Tr(T-T^) = If we denote the [U{2)l x U [2) R]giom x [U{2)v]iocai unitary 

transformations by {gi,gji,h), then the fields transform ^^ ji h{x)^L^jig\^ and i— >• 
h{x){V^ — id^)h^{x). How chiral symmetry restoration can come about in the HLS frame- 
work at large Np was discussed by Harada and Yamawaki |^ and at high temperature 



by Harada and Sasaki |£6[. Since their arguments are quite general, it is highly plausible 
that they can be applied to high density. This is what we shall do here. 

4.2 Vector Manifestation in Hot Matter 

To motivate our argument to be developed below for density-driven phenomena, 
we briefly summarize, without going into specific details which we will do in the next 
subsection, the work of Harada and Sasaki |^| on how the VM manifests itself as the 



critical temperature Tc is approached. We do this because this work brings out a subtlety 
in the behavior of vector-meson excitations in medium in general and in the neighborhood 
of the critical temperature in particular that has not been duly taken into account in other 
works on the matter available in the literature. 

There are two distinctive features in the HLS approach to the chiral phase transition. 
One is that in order to match with QCD while preserving hidden local symmetry at the 
chiral scale in the T = n = space 0, the p and u mesons must be considered 
as "light" in the sense that the pion mass is light as first pointed out by Georgi for a 
consistent power counting |Q. As noted below, it is in this sense that the HLS theory 



can be mapped to chiral perturbation theory [^|. This feature figures importantly in 



the hot (and/or dense) matter. The second feature is that in hot medium, there are two 
effects that control hadronic properties. One is the "intrinsic" dependence on T through 
the T-dependence of the QCD condensates (both quark and gluon) of hadronic masses 
and coupling constants and the other is the hadronic thermal fluctuations. The former is 
controlled by the matching to QCD and hence so is the latter through the parameters. 

In , Harada and Yamawaki studied the renormalization group structure of HLS 
in terms of the parameters of the HLS Lagrangian that enter at higher orders in the chiral 



15 



In this section, the matter density will be denoted n, reserving p for the p meson. 
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expansion and showed that effective low-energy theories with the given ffavor symmetries 
can ffow in various different directions as scale is varied but there is only one fixed point 
that is consistent with QCD, namely the Georgi vector limit. This implies that the 
intrinsic temperature dependence must be constrained by the Georgi vector limit, which 
means that the mass parameter My and the gauge coupling constant gy will approach 
zero as T — T^. When T is large, i.e.. My <^ T, |^ the pole mass of the vector (p) meson 
mass takes the simple form 

ml{T) = Ml{T)+glJ{T). (|.3) 

Here = af^gy where gy is the (hidden) gauge coupling and a ^ 1 signals deviation 
from the Georgi vector limit and J{T) is a finite temperature-dependent function. As 
T Tc, Mp and gy 0, so the pole mass goes to zero. This is consistent with BR 
scaling. 

We should stress that while the mass goes to zero at near Tc, the theory is nonetheless 
consistent with what one would expect at low T. For T <C Mp, the pole mass is found to 



be of the form ||66 



mJ(T) = M^2(T) + c^ + --- (|.4) 

n 

where c is a numerical constant and F.^ ~ 93 MeV is the pion decay constant in the 
vacuum. The ellipsis stands for terms higher order in T. Note that while this agrees with 
the low-temperature result of |6^, the pole mass need not increase with temperature as 
suggested in A close look indicates that it can even decrease at about the same rate 



as fniT) since Mp decreases faster than fn{T), suggesting the following - BR scaling - 
relation: 

mp{T)/mp ~ UT)/fn. (|.5) 

Note also that an early QCD sum-rule calculation by Adami and Brown of the thermal 
behavior of the p mass contains features that are similar to those in the vector manifes- 
tation. Although the connection to the Georgi vector limit is yet to be established, the 
vanishing of the p mass in the work of at the critical temperature Tc due to the 
approach to zero of the Wilson coefficient that multiplies a certain combination of gluon 
condensates with the vanishing of quark condensate resembles the way the p mass goes 
to zero in the vector manifestation. In the Adami-Brown formulation, the tadpoles which 
give the quark and antiquark masses in NJL are jammed together in going to deep Eu- 
clidean, so that to begin with, it appears that there is {qqqq) condensate. Adami and 
Brown factorized this, since they believed that before deep Euclidean was reached there 
were really two (qq) condensates. Whereas the question of "factorization" remained in 
the QCD sum rule approach, the present work of Harada and Sasaki shows that this was 
exactly right; namely, the scale A changes smoothly as chiral restoration is reached, with 
dynamically generated masses going to zero with it. 

-'^^This condition makes sense since My ^ as T — > Tc. 
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4.3 Vector Manifestation in Dense Matter 



In addressing dense matter, we follow the reasoning of |6^, §7j. Here we supply 
some details left out in the preceding subsection. We take the HLS Lagrangian as an 
effective Lagrangian that results when high-energy degrees of freedom above the chiral 
scale A taken to be higher than the vector meson mass mv are integrated out. Now the 
scale A will in general depend on the number of flavors Np, density n or temperature T 
depending upon what system is being considered. This is to be taken as a bare Lagrangian 
in the Wilsonian sense with the parameters gv{A.), a(A) and /7r(A) which plays the role of 
order parameter for chiral symmetry with = signaling the onset of the Wigner-Weyl 
phase. Reference [^] describes how these parameters can be determined in terms of QCD 
condensates by matching - a la Wilson - the vector and axial-vector correlators of HLS 
theory with the ones of QCD at the chiral scale A and how by following renormalization 
group (RG) flows to low energy scales, one can determine low-energy parameters that 
can be related to those that figure in chiral perturbation theory (for, e.g., vrvr scattering). 
An important observation here is that the (assumed) equality at chiral restoration (where 
(qq) = 0) of the vector and axial- vector correlators, i.e., 

ny|(gg>=0 = nA|(5-g)=0 (|-6) 

where 

iJd'xe"^-^0\Tr^{x)Jt{Om = S^\q,q.-g,.q')Uv{-q'), 

zJd'xe''^-{0\TJ^^{x)JlmO) = S'^\q,qu-9,uq')IlA{-q') (§7) 

implies that at the critical point, independently of how the symmetry change is driven, the 
HLS theory approaches the Georgi vector limit, namely, gy = and a = 1, in addition to 
the vanishing of f-j^. That chiral symmetry is restored with = at the Georgi vector 
limit follows from a proper account of quadratic divergence in the renormalization group 
flow equation for and a. Q In terms of baryon density n, we interpret this to imply 

^^The quadratic divergence is present when the cutoff regularization is used whereas it is absent at four 
dimensions Z? = 4 when the dimensional regularization is used. The quadratic divergence present in the 
cutoff regularization corresponds to a singularity at D = 2 in the dimcnsionally regularized integral from 
one-loop graphs which comes out to be proportional to the T function 

^ r(i - D/2) = - r(2 - D/3). 

In a renormalizable theory, this singularity (or any other power divergence) does not figure aX D — A but 
in an effective (nonrenormalizable) theory like HLS, this singularity present at D = 2 - which is absent 
at D — 4: - has to be subtracted. This procedure is quite similar to the subtraction of the si ngul arity at 



D — 3 that was required in the KSW scheme for two-nucleon scattering discussed in Section p.4| and we 
suspect that a similar physics is at work here. Note that were it not for this quadratic divergence, the 
/tt would not run and hence would not go to zero at the point where the Georgi vector limit is reached. 
Thus the quadratic term plays a crucial role in specifying the symmetry structure of the chiral restoration 
phase transition. 
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that at the critical density n = ric, we must have 



^y(AK);ne) = 0, a(A(nc);ne) = 1 (|.8) 

where we have indicated the density dependence of the cutoff A. The reasoning used here 
is identical to that for the temperature-driven transition. 

In HLS theory, the vector mass is given by the Higgs mechanism. In free space, it 
is of the form 



my = rup = = ^ a{niv)gv{mv)fn{mv) (|.9) 

where the cutoff dependence is understood. Here the parameter a(my) etc means that it 
is the value at the scale my determined by an RG flow from the hare quantity a(A). Note 
that (|[9) is similar, but not identical, to the KSRF relation mp^KSRF) = 
Now in medium with n ^ 0, this mass formula will remain the same except that it will 
depend upon density 



m,> = m-p = m-^ = ^ a[ml.)gv{my)f^{my). (|.10) 



The density dependence is indicated by the star. As in the case of Np discussed in |T^, 
the cutoff A will depend upon density, say. A* implicit in (|^.10). 



""^^However numerically they are very close. We explain as a side remark how this comes about. While 
Eq.(^.9) for mv — Mv{Q = rny) is approximate with numerically small corrections of order (7y/(47r)^ 
ignored there, the corresponding formula for (5 = 0, i.e., My(0) = a(0).g^(0)/^(0) is an exact low-energy 
theorem |7l[. The small corrections arise in going from the off-shell point Q = to the on-shell point 
Q ~ nip. Noting that the KSRF relation is given in terms of the p on-shell constant ^pTnr and the 
pion on-shell constant /7r(0), we first extract the former from HLS Lagrangian given in terms of on-shell 
quantities: 

9p^^ - 2 /2(0) 

where fa is the decay constant of the scalar that goes into the longitudinal component of the p by 
Higgs mechanism (which becomes equal to at the Georgi vector limit). The pion decay constant 
runs with pion loops below the p scale but fa- does not run below rup since the p decouples. Using 
fairUp) = a{m p) f'^{m p) , we can rewrite the KSRF mass formula as 



,2 ..2, a{mp)f,{mp) 

T-p{KSRF) - -^^^pi^p) 2 /2(0) 



By the renormalization group equations given by Harada and Yamawaki, we have a{mp) ~ 1 and ^"j^q^ 
2. Thus, we get 

"^liKSRF) - Mpirup) EE ml. 



Note that the mass in this formula is a parameter in the Lagrangian and not necessarily a pole mass 
except at the phase transition point. The pole mass will generically have additional density-dependent 
terms which will go to zero as the coupling gv reaches the Georgi vector limit. 
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The Harada-Yamawaki argument (or theorem) suggests that at n = nc where {qq) = 
and hence Uy^ric) = n^(nc), the Georgi's vector hmit gy = 0, a = 1 is reached together 
with = 0. This means that both the mass parameter (^.10) and the pole mass vanish, 

m^K) = 0. §.11) 

At the point where the vector meson mass vanishes, the quartet scalars will be 
"de-Higgsed" from the vector mesons and form a degenerate multiplet with the triplet of 
massless pions (and "r/" 0) with the massless vectors decoupled. This assures that the 
vector correlator is equal to the axial-vector correlator in the HLS sector matching with 
the QCD sector. In this scenario, dictated by the renormalization group equations, the 
vector meson masses drop as density (or T) increases. Note that this scenario is distinct 
from the "standard" (as yet to-be-established) picture in which the p and ai come together 
as do the pions and a scalar a. In the standard scenario, there is nothing which forces the 
vector mesons to become massless and decouple. They can even become more massive at 
chiral restoration than in the vacuum as discussed by Pisarksi |^. Thus the vanishing 
of the vector-meson mass is a prima facie signal for the phase transition in the Harada- 
Yamawaki picture. 

Here using the same reasoning as for large Np and high T, we have arrived at the 
conclusion that the same vector manifestation can take place at the critical density 0. 
As in the case of temperature, the "running" in density of the parameters gy, a and 
implies that away from the critical point, the vector meson (pole) masses drop as density 
is increased. In the next section we will interpret this phenomenon in terms of BR scaling. 

In this section as well as in what follows, we are interpreting the vector-meson 
"mass" in the sense of the BR-scaling mass. Of course in the presence of a medium with 
nonzero temperature and/or density, the lack of Lorentz invariance gives rise to different 
components and it will be necessary to specify which component gets "de-Higgsed" at the 
phase transition. The meaning of the vector mass as used here is not the same as the 
standard definition as we will try to clarify in the following sections. 



5 LANDAU FERMI LIQUID FROM CHIRAL 
LAGRANGIANS 



^°We assume that the fourth pseudoscalar "77" gets a mass by U{1)a anomaly, so can be excluded from 
our consideration. 

^^An important point to note here is that Pisarski's argument hinges on the vector dominance picture. 
However Harada and Yamawaki show that at the vector manifestation, the vector dominance is strongly 
violated [Q. This is one of the basic differences in treating vector mesons. 

^^There is a puzzle here. The analysis of for the TV^-driven phase transition is that the top- 
down phase transition from Np > Npc to Np < Npc is neither first-order nor second-order. The 
order parameter is found to be continuous at the phase transition, so it is not first-order. However the 
correlation length does not diverge at the phase transition, so it is not second-order either. If correct, 
this phase transition may appear to be different from that expected in temperature or density. 
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5.1 Fluctuating Around Zero-Density Vacuum 



We have developed the thesis that chiral symmetry can be suitably implemented 
into nuclear problems involving a few nucleons. Furthermore EFT can be formulated in 
such a way that the traditional nuclear physics approach with realistic potentials (i.e., 
PM) can be identified as a legitimate part of a consistent EFT supplying the leading 
term in the expansion. This means that what nuclear physicists have been doing up to 
today - and their highly successful results - can be considered as an essential part of the 
modern structure of strong interactions based on QCD. For the few-nucleon processes we 
have been dealing with, e.g., NN scattering, n + p ^ d + solar p + p^d + e'^ + Ue 
etc., what enters into the theory is the chiral symmetric Lagrangian determined in free 
space, that is, at zero density. The starting point there is an effective chiral Lagrangian 
that describes QCD in the matter-free vacuum, with masses and coupling constants all 
determined in particle experiments in medium-free space. The few-nucleon problems are 
then treated by fluctuating around the matter-free vacuum near a "fixed point" (at infinite 
scattering length) discussed above. The machinery for doing this is the well-established 
chiral perturbation theory (xPT) involving baryons. Since the processes involved occur 
near the fixed point of the theory, one can successfully do the perturbation. 

Now how do we go to heavier nuclei, nuclear matter and denser matter which must 
live near a different fixed point, say, the Fermi-liquid fixed point to be described below 
which is not necessarily near the above fixed point? This is the principal question that 
we wish to address in this review. If there are no phase changes along the way, one 
may attempt to build many-nucleon systems in a perturbative scheme starting with the 
free-space Lagrangian mentioned above. Indeed there have been several such attempts 
in xPT to reach at nuclear matter ||7^, |75[. It is not obvious though that this can be 
done beyond the normal matter density since a strong-coupling expansion is involved and 
large anomalous dimensions arise. This means that it may be simply too hard if not 
impossible to handle various phase changes (kaon condensation, chiral restoration, color 
superconductivity etc.) that are conjectured to occur as one increases density toward chi- 
ral restoration as demanded for relativistic heavy-ion processes and for describing compact 
stars. A different but perhaps more precise way of saying the same thing is in terms of 
the skyrmion picture for dense matter which is closest to QCD if Nc is taken to be large. 
In the skyrmion description, a lump of baryonic matter with a given baryon number B is 
characterized by a "topological vacuum" associated with a conserved winding number B. 
A system with density p'^ ^ ps means that B' ^ B for a given volume. Now the point 
is that the B' vacuum may not be simply connected to the B vacuum by a small pertur- 
bation. In what follows, therefore, we shall take a different route and formulate a chiral 
Lagrangian field theory using a "sliding" vacuum with parameters of the Lagrangian run- 
ning with matter density. This leads to the notion of Brown- Rho scaling and in-medium 
effective field theory for a system with Fermi seas, i.e.. Landau Fermi-liquid theory. 
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5.2 Skyrmion vs. Q-Ball 

Consider bound many-nucleon systems corresponding to nuclei. If the number of col- 
ors can be taken to be very large, then one may approach them starting with a skyrmion- 
type Lagrangian constructed of Goldstone boson fields (and heavy mesons) and looking 
for solutions of the winding number W equal to the baryon number B or the mass number 
y4 in a nucleus. The Skyrme Lagrangian which consists of the current algebra term and 
the quartic term 

Csk = ^Trid^Ud^U^) + ^Tr[Ud,U\ UdM^]' + ■■■ 

which may be viewed as an approximate zero bag size limit of the Lagrangian (^1) has 
been extensively studied mathematically for baryon number up to S = 22 with a fas- 
cinating result [^. (The ellipsis stands for higher derivative and/or heavy-meson and 



chiral-symmetry breaking terms that are to be added as needed.) The Lagrangian of the 
type (^.1) suitably extended with higher derivatives and/or massive boson fields is accu- 
rate only at large- A^c limit, so it is not clear what the classical solution of this Lagrangian 
for B > 2 represents since in nature where iVc <^ oo. Nonetheless this development 
is quite exciting since it is closest to QCD at least for Nc = oo. As stated above, the 
ground states with baryons numbers B and B' with B ^ B' cannot be connected by 
small perturbation because they represent different topological "vacua." 

A direction we will adopt here is not as ambitious as the above but we will attempt 
to account for the different vacuum structure for dense systems (heavy nuclei) than for 
dilute systems (light nuclei). We shall start with a theory for a system with a Fermi sea. 
It is possible that such a system emerges as a sort of nontopological soliton or "Q-ball" 
from a theory based on chiral perturbation theory defined at zero density as suggested 
first by Lynn |]7^ and recently refined by Lutz et al |7^. Although there is no proof. 



we believe that the Q-ball system which is non-topological is essentially equivalent to the 
skyrmion system which is topological. This equivalence is somewhat like in the nucleon 
structure where the two pictures work equally well. The essential point we shall accept 
here is that a non-topological soliton solution exists possessing the liquid structure that 
is supposed to be present in nuclear matter. We will later identify this with what we do 
know but for the moment let us assume that we have a Fermi sea filled with nucleons. It 
is fortunate for our purpose that such a system is accessible to an elegant and powerful 



effective field theory treatment as shown by various authors that leads to Landau's 



Fermi liquid theory applied to nuclear systems WE . 



5.3 Effective Chiral Lagrangian for Many-Nucleon Systems 



What we wish to do in our case is to arrive, following the developments in ||79|, |81|] , 
at the Landau Fermi liquid theory starting from an effective chiral Lagrangian that we 
have been discussing above. Treating the nucleon in terms of a local field denoted by 
and the Goldstone boson field by ^ = VU = e^'^'^l'^^^ ^ we can write a simple Lagrangian 
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of the form 



L = N[t-f^{d^ + tv^ + g\^,a^) - M*]N - ^ QiNT^Nf + ■■■ (|.2) 

i 

where the elhpsis stands for higher dimension nucleon operators and the Fj's Dirac and 
flavor matrices as well as derivatives consistent with chiral symmetry. The star affixed on 
the masses and coupling constants will be defined precisely below. The induced vector and 
axial vector "fields" are given by = —^{^^df,^ + ^d^^^^) and = —^{C.'^d^C, — C.d^C.'^). 
In (§.2) only the pion (vr) and nucleon {N) fields appear explicitly: all other fields have 
been integrated out. The effect of massive degrees of freedom will be lodged in higher- 
dimension and/or higher-derivative interactions. The external electro- weak fields if needed 
are straightforwardly incorporated by suitable gauging. 

When applied to symmetric nuclear matter in the mean field approximation, the 
Lagrangian (^.2) is known to be equivalent ||8^, |8^ to the Lagrangian that contains just 



the degrees of freedom that figure in a linear model of the Walecka-type ^ 

C = Nii-f^id^" + ig*uj^') - + h''a)N 

--/I. + \{d,af + - + ■ • • (§3) 

where the ellipsis denotes higher- dimension operators. We should stress that (|.3) is 
consistent with chiral symmetry since here both the u and a fields are chiral singlets. 
In fact the a here has nothing to do with the chiral fourth-component scalar field of 
the linear sigma model except near the chiral phase transition density; it is a "dilaton" 
connected with the trace anomaly of QCD and is supposed to approach the chiral fourth 



component in the "mended symmetry" way a la Weinberg [|60| , |61[| near chiral restoration. 
In practice, depending upon the problem, one form is more convenient than the other. In 
what follows, both (|^.2) and (^3) will be used interchangeably. 

In (||.2) and (|^.3), the mass and coupling parameters affixed with stars depend upon 
the "sliding" vacuum defined at a given "density." (Putting density into the parameters of 
a Lagrangian is subtle because of chiral symmetry and thermodynamic consistency. We 
will specify more precisely how the "density" is to be defined.) We are thinking of the 
situation where an external pressure (e.g., gravity for compact stars) is exerted and hence 
we can think of the parameters of the Lagrangian varying as a function of density. They 
need not be associated with quantities defined at a minimum of an effective potential of 
the theory without the external pressure. 



The next important observation we shall exploit is that by Matsui ||8^ that a La- 
grangian of the type (|^.2) can give in the mean field the result of Landau Fermi liquid. 
This means that in the same approximation, the chiral Lagrangian (^2) can be mapped 

^•^The ^ field appearing here is the "gauge-fixed" one, £,l = £,li — £,■ Here we are dealing with the 
matter far away from the chiral phase transition point, so we are taking the unitary gauge. 

^^For asymmetric nuclear matter, isovector fields (e.g., p, ai etc.) must be included in a chirally 
symmetric way. 
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to Landau Fermi liquid. From the works in which tell us that Landau theory is 
a fixed-point theory with a Fermi-liquid fixed point which is distinct from that of the 
effective Lagrangian used for two-nucleon systems discussed above, we learn that the the- 
ory (1^.2) in the mean field has two fixed-point quantities, viz, the effective quasiparticle 
(or Landau) mass m% and the Landau parameters J-' for quasiparticle interactions. The 
full interaction between two quasiparticles Pi and P2 at the Fermi surface of symmetric 
nuclear matter written in terms of spin and isospin invariants is defined by 



^ 1 



iV(0) 



F(cos^i2) + F'(cos6'i2)ti ■ T2 + G{cos6i2)(Ti ■ a-2 



q2 

-f G'(cos6'i2)cri ■ CT2T1 ■ T2 + -r2--f^(cos6'i2)S'i2(g) 

+ TjH'{cOSei2)Si2{q)Ti ■ T2 
Kp 



(1-4) 



where 6*12 is the angle between p^ and P2 and iV(0) = (^)f ^'^^ density of 
states at the Fermi surface. Also, q = Pi — P2 and 5*12(9) = 3cri • qcr2 ■ q — <ti ■ 0-2, 
where q = q/\q\. The tensor interactions H and H' are important for the axial charge 
to be considered later. The functions F, F', . . . are expanded in Legendre polynomi- 
als, F(cos6'i2) = J2e FiPi{cos 612) , with analogous expansion for the spin- and isospin- 
dependent interactions. 

The fixed point is associated with a given Fermi momentum kp or density p of 
the given system, so that there will be a set of fixed points JF for each (given) kp or p. 
Thus we will have the Landau parameters that depend upon density once we accept that 
the Landau mass depends on density as is obvious from the Landau mass formula 

©.5) 

where Fi = (M/m^)Fi and Fi is the / = 1 component of the interaction F in (^.4). 



5.4 Brown-Rho Scaling 

How the fixed point quantities vary with density cannot be determined from the 
effective Lagrangian alone. It must be tied to the vacuum property characterized by the 
chiral condensate (ipip) where ip is the quark field. The dependence must therefore be 
derived from the fundamental theory, QCD. In the past QCD sum rules have been used 
to extract information on this matter. For our purpose, it is more convenient to take the 
Skyrme Lagrangian (^.1) as a starting point and try to deduce the necessary information. 
We expect it to give the leading behavior corresponding to a mean field theory in the 
sense to be defined more precisely later. This is the argument presented in the original 
BR paper . Implementing the scalar dilaton field x associated with the trace anomaly 
of QCD to the Lagrangian extended to take explicitly into account the vector mesons 
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= Pf^, 00^, it has been found that in the mean field the following scaling holds, e.g., in 
the Lagrangian (^.S) 0: 

^'^'^ U ^ ^ my ^ M ■ ' 

Here fj^ is the pion decay constant, the /* being the in-medium one |^. We caution that 
the M* is the nucleon in-medium mass (or BR scaling mass) which is not the same as 
the Landau mass as we will see shortly. The quantity related to the quark 

condensate $ ~ {{MY I {'ll)T where a is a constant which we will find in Section | equals 
~ 1, is the scaling factor that needs to be determined from (fundamental) theory or 
experiments. 

Let us now examine more closely what the "density dependence" in general (not 
necessarily in the context of BR scaling) means in the Lagrangian (^2) or (^.3). For this, 
let us take (^.3) and introduce the chirally invariant operator 

pu^ = N-f^'N (|.7) 



where u'^ = (vT^^i^)^^(l, v) = (y — J^)^^(p, j) is the fiuid 4-velocity. Here j = 
{NjN) is the baryon current density and p = {N'^N) = J2i the baryon number density. 
The expectation value of p yields the baryon density in the rest-frame of the fiuid. Using 
p it is easy to construct a Lorentz invariant, chirally invariant Lagrangian with density 
dependent parameters. 

Now a density dependent mass parameter in the Lagrangian should be interpreted 
as m* = m*(p). This means that the model (|§.3) is no longer linear. It is highly non-linear 
even at the mean field level. 

The Euler-Lagrange equations of motion for the bosonic fields are the usual ones 
but the nucleon equation of motion is not. This is because of the functional dependence 
of the masses and coupling constants on the nucleon field: 



with 



[i-f^{d^ + igluo^ - lu^t) - + /iV]iV = (|.8) 



^^The scalar tr as well as the vector V are absent in (||.2) since they have been integrated out. Their 
scaling will appear therefore in the coefScients representing those degrees of freedom, namely, C*'s. The 
connections between the two are given in |7^, |8^ . 

^^Since Lorentz invariance is lost in medium, the constants /*, g\ etc. contain two components, space 
part and time part, that are in general different. The parameters that appear in BR scaling are mean- field 
quantities and are defined more precisely in the sense described above for a generalized Walecka model 
(^.3). The same remark applies to the in- medium vector coupling we will discuss in later sections. 
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Here we are assuming that {d/dp)h* ^ 0. This can be justified using NJL model [|120| 



but we will simply assume it here 0. The additional term E, which may be related to 
what is referred to in many-body theory as "rearrangement terms" , is essential in making 
the theory consistent. This point has been overlooked in the literature. 

The crucial observation is that when one computes the energy-momentum tensor 
with (^3), one finds, in addition to the canonical term (which is obtained when the 
parameters are treated as constants), a new term proportional to S 

T'-' = + mu ■ ^N)g>^r (|.10) 

The pressure is then given by \{Tii)v=Q. The additional term in (|^.10) matches precisely 
the terms that arise when the derivative with respect to p acts on the density-dependent 
masses and coupling constants in the formula derived from Too: 

v = -^ = P^ = w-^ i-ii) 

where 

E = (T°°). (|.12) 

This matching assures energy-momentum conservation and thermodynamic consistency. 
This verifies that interpreting "density-dependent" parameters in our Lagrangian as the 
dependence on p with (p) being the density is consistent with both chiral symmetry and 
thermodynamics. 



5.5 Landau Mass and BR Scaling 

We end this section by writing down the relations between the parameters of the 
Lagrangian (|.2) or equivalently (|^.3) and Landau parameters. We will use the former 
and take (including the isovector channel) 

~J:c:m^Nr ^ -^{n^^ny - ^{N^.rNy + .... (g.i3) 

This may be considered as the leading terms that result when the vectors u and p are 
integrated out. But it is more than just that. In fact, the subscripts can be taken to 
represent not only the vector mesons uj and p that nuclear physicists are familiar with but 
also all vector mesons of the same quantum numbers (which account for the appearance 
of the tilde on the vector mesons), so the two "counter terms" on the right-hand side 
subsume the full short-distance physics of the given chiral order. Since the Lagrangian 
(^.2) in tree order must correspond to the mean field of (^.3), one can simply compute the 

^^In fact apart from the phenomenological requirement, there is no reason to assume this. Ah our 
arguments that follow would go through without difficulty even if the coupling constant were to scale 
with density. 
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relevant Feynman graphs with (^2) and match them to the Landau Fermi hquid formula. 
One obtains |^ 



where -Fi(vr) is the pion contribution to Fi which is completely fixed by chiral symmetry 
so the only parameter that enters in (^.14) is the BR scaling $ M*/M. 
Now using the Landau mass formula (^-5) and 

F, = F,{u) + F,{7r) (|.15) 

we find 

Fi(cD) = 3(l-<l>-i) (1.16) 

where -Fi(u)) includes contributions from a// excitations of the u meson quantum numbers. 
This is the main result of the connection between the chiral Lagrangian with BR scaling 
and Landau Fermi liquid theory. This is a highly non-trivial relation which says that BR 
scaling - which in the context of the QCD vacuum refiects modification due to matter 
density - is tied to short-range quasiparticle interactions. This means that what may be 
considered to be a basic feature of QCD and nuclear interactions via nuclear forces are 
inter-related, indicating another level of the Cheshire Catism. 

• Determining $(p) 

If one accepts the BR scaling (^6), there are several sources to use for determining 
the scaling factor For instance, one can determine it from the QCD sum rule ^for 

the mass ratio m*/mp or from the ratio f*/f-n using the Gell-Mann-Oakes-Renner 
(GMOR) formula for pion mass in medium or a fit to nuclear matter properties ||79| . 
The result comes out to be about the same at nuclear matter density 

$(p = Po) ~ 0.78. (1.17) 

In terms of the parameterization of the form 

*(p) = (1 + yp/p^r\ §18) 

^^There is a great deal of controversy in the application of QCD sum rules in dense and/or hot medium. 
As it stands, the situation is totally confused and it is perhaps too dangerous to take seriously any result 
so far obtained. For a recent discussion on this matter, see Our proposal is that the value (||.17) 

which will result from the fit to a variety of nuclear processes discussed below be taken as the value to 
be obtained from the QCD sum rule approach. 

^^Using the in-medium GMOR formula to arrive at the value (|^.17), it is assumed that the pion mass 
remains unaffected by density up to nuclear matter density. This assumption may have to be re-examined 
in view of the "pion mass" measured in deeply bound pionic states which indicate a mild increase as a 
function of density as discussed in Section g|. This issue is not so simple to address since it is not clear 
that the mass that appears in the GMOR relation is the same quantity that has been "measured." 
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the result corresponds to y{p = po) ~ 0.28. This gives the Landau mass m)^/M ^ 0.70 
which is to be compared - with a due consideration of the caveat mentioned above on 
the difficulty with QCD sum rules in nuclear physics - with the presently available QCD 
sum-rule value 0.69to;o7 [|^. This value is also consistent with the properties of heavy 
nuclei. 



6 INDIRECT EVIDENCES FOR BR SCALING 

6.1 Indications in Finite Nuclei 

Before the publication on BR scaling in 0, Brown and Rho proposed |89, 90| that 



some of the missing strength in the longitudinal response function in nuclei could be 
explained by the dropping masses in medium m*/mp ^ m^j^/rriN- ^ 0.75. The same 
dropping masses (of both the nucleon and the p) could also explain that within certain 
range of kinematics relevant to experiments, the transverse response function does not 
scale in medium. This is essentially the picture that is conffimed by the Morgenstern- 
Meziani analysis mentioned in the introduction and described in detail in Section |^. 

In addition to the electromagnetic response functions, there were further predictions 
in that followed from the dropping masses. One was the anomalous gyromagnetic 
ratio 6gi in nuclei which gets enhanced by the scaling factor. This will be discussed below 
in terms of a more modern formulation with chiral Lagrangians and Landau Fermi liquid 
theory although the physics is essentially the same as in this work. The second is that the 
dropping mass implies that as density increases, the tensor force in nuclei gets weakened. 
The reason is that the tensor force is given by two agencies, one-pion exchange and one-p 
exchange coming with opposite signs. At low density, the pion exchange dominates, so 
the net effect is more or less controlled by the pion tensor. However as density increases, 
the p tensor gets enhanced by the factor {mN /to^nY apart from the increased range of the 
p-exchange force. Therefore the p tensor tends to cancel more and more the pion tensor 
as density increases, thereby suppressing the net tensor force in the outer region where 
pion exchange dominates. There are some indications that the suppressed tensor force 
is consistent with nuclear spectra as well as with Gamow- Teller transitions in nuclei to 
which we shall return shortly. 

There is, however, a possible caveat to this as we will see later and that has to do with 
the possibility that the vector couphng g*^j^ will also drop as density increases, eventually 
decoupling as dictated by the "vector manifestation" discussed above. It is not clear at 
what density this decoupling sets in but there are empirical indications that the tensor 
force weakening is operative at least up to nuclear matter density. The phenomenology 
in nuclei discussed above also indicates that at least at the low densities encountered in 
nuclei, the ratio gpNN/m*p does increase. It should be noted that even when the vector 
decoupling takes place, with the flavor symmetry in the hadronic sector yielding to the 
color gauge symmetry in the quark-gluon sector |^, some vector repulsion remains coming 
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The way this change-over takes place is discussed in Section |ll|. 
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from (Fierzed) gluon exchange 



6.2 The Anomalous Gyromagnetic Ratio in Nuclei 

The development given above in Section]^ of the effective chiral Lagrangians with BR 
scahng and Fermi hquid fixed point theory can provide a more sophisticated and rigorous 
explanation for the enhanced gyromagnetic ratio gi in nuclei mentioned above. We discuss 
this as a strong, though indirect, evidence for BR scaling. The nuclear gyromagnetic 
moment gi in the convection current for a particle sitting on top of the Fermi sea is 
defined as 

We have defined ^f; = 1 + 6gi such that the bare nucleon mass M, not the effective 
(Landau) mass m%, appears in the formula. This is because we want to preserve charge 
conservation associated with the factor 1 in gi, 6gi being purely isovector. 

The 6gi has been determined from experiments. The measurement that is most 
relevant to our theory is the one on giant dipole resonances on heavy nuclei, in particular 



in the lead region. The result of ||9T|] on ''^Bi gives the proton value 



Sgf = 0.23 ±0.03. (|.2) 

This can also be determined from magnetic moment measurements but the analysis is 
somewhat more complex and hence less precise. We shall refer to (^.2). 

Given the effective chiral Lagrangian (|^.2), with the connections we have established 
so far, it is easy to compute in the mean field order the nuclear gyromagnetic ratio gi. It 
suffices to calculate all terms of the same order to assure that charge is conserved or what 



is equivalent, "Kohn theorem" is satisfied. The result is the Migdal formula fTg] given in 
terms of the parameters of the Lagrangian, that is, the BR scaling 

91 = ^ + Sg, (§.3) 



with 



5gi = ^{Fi-F,)Ts=^ 



- 1 - l^in) 



As mentioned above, -Fi('/r) is completely fixed by chiral symmetry for any density 
- possibly up to chiral restoration - so the only quantity that appears here is the BR 
scaling factor. At nuclear matter density, we have |-Fi(7r)|p=pg = —0.153. Now taking 
$(po) ~ 0.78 from (|.17), we predict from (|.4) 

Sgj^ = 0.227r3. (|.5) 

This agrees with the experiment (P.2) providing a quantitative check of the BR scaling. 



48 



6.3 Axial Charge Transitions in Heavy Nuclei 

Another process where the Lagrangian (^.2) with BR scahng together with the 
PKMR's EFT can make a prediction is in the axial charge beta transition of the type 

A{J^)-^A'{J^)+e-{e+) + u{u) AT = 1. (|.6) 

This process was studied by several experimentalists with the objective of exhibiting 
strong medium-enhancement of the matrix element of the axial charge operator that 
governs the process (§.6) in heavy nuclei. Such an enhancement was first observed by 
Warburton in the lead region ||^ and since then several authors confirmed Warburton's 
finding in various other heavy nuclei. Warburton focused on the quantity (that we shall 
refer to as "Warburton ratio") 

eMEC = M,^p/M,p (|.7) 

where M^^p is the measured matrix element for the axial charge transition and M^p is 
the theoretical single-particle matrix element for a nucleon without BR scaling. There are 
several theoretical uncertainties involved in extracting the Warburton ratio. First, one 
has to extract from given beta decay data what corresponds to the axial charge matrix 
element. This involves estimating accurately other terms than the axial charge that 
contribute. The second is what one means precisely by Mgp which is a theoretical entity. 
For these reasons, an unambiguous conclusion is difficult to arrive at. However what is 
significant is Warburton's observation that in heavy nuclei, eMEC can be substantially 
larger than the possible uncertainties involved: 

^HeavJNucle^ ^ _ 2.0. (|.8) 

Furthermore more recent measurements and their detailed analyses in different nuclei [93, 
94| , p5| as shown in Table quantitatively confirm this result of Warburton. 



The prediction from the Lagrangian (^.2) with BR scaling in the same approxima- 
tion as in the gyromagnetic ratio case is [^, |7^ 



e^^i%c = + ^) (1-9) 

where A is, in accordance with the chiral filter mechanism, dominated by the pionic 
contribution with small controllable corrections from vector degrees of freedom the 
magnitude and the density dependence of which are again totally controlled by chiral 
symmetry, with the BR scaling $ appearing as an overall factor. For p = po/2 and 



p = po, we find 



e^fEC7lp=Po/2 ^ 1-63, e^f^^|,=,, ^ 2.06. (|.10) 
This is in an overall agreement with all the available Warburton ratios (see Table P). 
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In terms of the Landau-Migdal fixed-point parameters, A — G'l/i — 10Hq/3 + 4iJ(/3 — 2H2/I5 with 



the dominant contribution coming from the soft pions |51 
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Table 2: "Empirical values" for cmec 



Mass number A eMEC Reference 



12 


1.64 ±0.05 






50 


1.60 ±0.05 






205 


1.95 ±0.05 


m 




208 


2.01 ±0.10 


d2] 



6.4 Axial- Vector Coupling Constant gj^ in Dense Matter 

One can make a rather simple statement on how the axial- vector coupling constant 
g\ scales in dense matter on the basis of BR scaling and the Fermi liquid theory. 

The story of g\ in nuclei is a long story and rather involved with nuclear struc- 
ture effects compounded with strong nuclear tensor correlations, excitations of nucleon 
resonances and possible effects of the vacuum change ^ ^ |100|| . Gamow- Teller 



transitions in complex nuclei, particularly giant resonances thereof, play an important 
role in presupernova stage of the collapsing stars and the effective coupling constant g\ 
figures importantly in their description [|101|| . 

What we will present here is probably not the unique explanation of what is going 
on but it is a version that is immediate from, and consistent, with the general theme of 
this review 

We have arrived above at the Landau mass formula (§.14) given in terms of the 
BR scaling parameter $. The Landau mass differs from the BR scaling mass M* due to 
the pion contribution. From the Skyrme-type Lagrangian from which we have deduced 
the scaling behavior, this means that there must be contributions from other than the 
property of /^r. The only other factor that enters in the skyrmion description is the 
coefficient of the Skyrme quartic term e or in terms of a physical quantity, qa- When this 

^^It is perhaps appropriate at this point to clarify how the absence of Lorentz symmetry manifests 
itself in our treatment of g\, an issue we briefly alluded to earlier. As the two (preceding and this) 
subsections clearly show, the effective axial coupling constant gA in matter is different for the space and 
time components of the axial current. We see that the time component of the axial charge (measured 
through axial-charge transitions) is enhanced in nuclear matter whereas the space component of the axial 
charge (measured through Gamow- Teller transitions as here) is quenched. Although we did not make the 
distinction in discussing BR scaling, it is evident that the 5^ that appears in the relation between the 
BR scaling mass M* and the Landau mass involves the space component measured by Gamow- Teller 
transitions. A similar separation can be made for instance for the pion decay constant /* which will be 
discussed in Section ^. 
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is taken into account, the resulting expression is found to be Q 



(-AT 



M 

from which we get 




jii) 



At nuclear matter density, this predicts 

9\{Po)^l. (|.13) 

What does this mean physically? If we accept that the Landau mass is a fixed point 
of the Fermi liquid theory, then the g\ must be a constant defined at the fixed point. It 
therefore must correspond to decimating all the way down to the Fermi surface, that is 
to say, sending the infrared cutoff kc in the renormalization- group equation to zero. In 
this limit, g\ must be proportional to the Landau parameter G'q. This connection in 
infinite (translationally invariant) matter is being investigated in ||102|| . What this means 



in practice in finite nuclei in which the transition is measured is not yet fully understood. 
But it appears reasonable to suppose that applied to finite nuclei, what we obtained here 
is an effective constant that should be used when the Gamow- Teller matrix element is 
computed within the valence shell with the single-particle states undergoing the Gamow- 
Teller transition restricted to the lowest-lying excitations on top of the Fermi sea. In other 
words, the g\ so obtained is a fixed-point quantity. 

An apt illustration of how this reasoning seems to be working is given by the com- 
plete 0/lp-shell Monte Carlo calculation of Gamow- Teller response functions by Radha 
et al [p.03|| and also big shell-model diagonalization by Caurier et al ||104|| . Remarkably, it 



is found that (7^ ~ 1 is universally needed in these calculations [|103| , |105| , |101|| which may 
be taken as a support for the fixed-point notion of the g\ . 



6.5 Evidences from "On-Shell" Vector Mesons 

The indirect evidences we have discussed above dealt with the properties of hadrons 
that are way off-shell, with ^ 0. The masses and coupling constants that are probed in 
such processes correspond to the parameters of the effective Lagrangian and they would 
be near their physical parameters only if the tree approximations were valid. One would 
of course like to confront with measurements of the "physical hardrons" propagating 
in dense matter. In order to do this quantitatively with a theory that implements BR 
scaling, one would have to compute to higher orders, say, in chiral perturbation suitably 

^•^Since the middle expression of (^.12) involves the pionic contribution _F'i(7r), one might think that g\ 
differs from qa merely due to the presence of the pion. This interpretation is not correct. The Landau 
parameter Fi{-k) contains the density of states A^(0) = 2kF'm^ /tt'^ which carries the Landau mass 
and hence does involve BR scaling. This is clearly seen in the last equality. 
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formulated to incorporate the "sliding" vacuum structure associated with BR scaling. 
This would of course give rise to widths of the hadrons involved and to a further shift in 
mass from that of BR scaling which is mean field. Such a program has not yet been clearly 
formulated. Hence there are no predictions in our approach. For low density, there are 
several calculations that start from the zero density vacuum ||106|| that predict that while 
the p meson may or may not be shifted in mass, the uj will. The p will broaden in width in 
matter but the uj will remain narrow and may even be bound in nuclear matter ||107|| . The 
experimental finding of this property of uj will confirm unambiguously the observations 
made in the preceding and following sections. 

Experiments are presently being performed to check all this at various laboratories, 
notably at GSI and Jefferson Lab. Specifically, the signals from heavy-ion collisions to 
be discussed below address this issue although with results that are compounded with 
temperature effects. 
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Figure 8: Distribution with invariant mass of e'^e'^ pairs [109]: a) from a carbon target; b) 
from a copper target. The difference between histogram and thin solid line represents the effect 
of medium (density) on vector-meson mass, the e~*"e~ invariant mass being determined by the 
in-medium vector-meson mass. 



Here we briefly comment on some recent experiments that seem to indicate that 
the vector-meson masses do decrease with density in nuclear matter or that at least the 
dropping masses are compatible with the observations. One is the TAGX collaboration 
result ||108|| on the process ^He(7,p°)ppra with the tagged photons in the range 800-1120 
MeV. The mass shift of the p^ meson reported was Snip = 160 ± 35 MeV. This is nearly 
the amount predicted for nuclear matter, i.e., mp(l — <I>) ^ 170 MeV for $ = 0.78 used 
above. This seems to be a bit too much a shift for such a small nucleus. It seems to involve 



other than just BR scaling. The other experiment is a more recent measurement ||109|| 



at KEK on the invariant mass spectra of the e^e pairs in the target rapidity region 
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of 12 GeV p + A reactions with the nucleus A = C and Cu. The experimental result 
is given in Fig.^]. This collaboration "sees" the mass shift of the vectors in the heavier 
nucleus indicative of a BR scaling. Although the precise value of the mass shift is not 
determined in this work, it is consistent with the expectation. We learned from private 
communication ||110|| that with their kinematic distribution, about 60% of the p mesons 
and about 10% of the u mesons are estimated to decay inside Cu nucleus, so most of the 
excess on the low-mass side of the uj peak is from the p mesons. The authors of [|109|] 
considered the possibility of an in-medium increase in u mass, but this did not give a 
statistically significant effect. More data and their analyses are forthcoming and will give 
a clearer insight into the properties of the hadrons on the "mass-shell," thereby checking 
the BR scaled effective Lagrangian in more details. 



7 DIRECT EVIDENCE FROM THE (e, e'p) 
RESPONSE FUNCTIONS IN NUCLEI 

The evidences for BR scaling given so far are more or less indirect. They are based on 
consistency with what we already know of nuclear interactions. In this section, we discuss 
what we consider to be the most direct evidence for BR scaling, namely, the longitudinal 
and transverse electromagnetic response functions in medium-heavy and heavy nuclei. 
Although the process involves off-shell properties of the vector meson, we believe that this 
evidence is cleaner and more direct than what we can extract from presently available 
heavy-ion data - including the celebrated CERES data - described below. 

Electron scattering which knocks nucleons out of nuclei have been used for many 
years to tell us about the binding energy of protons in nuclei and also to pin down the 
effective mass that the proton has in the nucleus. To the extent that the electron couples 
through vector mesons to the nucleons, these experiments can also tell us about the in- 
medium properties of the vector mesons. This was worked out first in a paper by Soyeur, 



Brown and Rho ||111|| following the initial suggestion in [q^, We follow here a more 
schematic treatment [ p.l2|| which brings out the main ideas more clearly, and which is in 
semi-quantitative agreement with the more detailed calculation. 

Before going further, let us remark that the experimental analysis has been thor- 
oughly confused during the 1990's and only now is being clarified, chiefly due to efforts of 
Joseph Morgenstern. A global analysis by Jourdan | |113[ which relied on three-dimensional 
numerical calculation of Onley et al. ||114|| indicated that medium effects appeared to be 
absent, or at least the situation was unclear. The analyses involved therein were un- 



necessarily complicated in view of the fact that Yennie et al ||115|| had shown that the 
chief distortion by the Coulomb field could be reliably taken into account by shifting the 
momentum q to an effective momentum 



q,ff = {l-Vj\q\) 



(0-1) 



where Vc is an effective Coulomb potential of the electron. This is a semiclassical correction 
and had to be accompanied by higher order corrections for a quantitative description of 
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Figure 9: The longitudinal response function Si, as a function of q^ff using (a) Saclay data only 
and (b) Saclay data plus SLAC NE3 and Bates data. The SLAC NE9's ^^Fe result is shown 
by a cross and that of the Jourdan analysis by a star in (b). The theoretical curves are: the 
solid black curve for nuclear matter without BR scaling and the dashed curve for ^He; the long 
dashed curve is the same as the nuclear matter one integrated within the experimental limits of 
cj. The dot-dashed curve is with BR scaling. See Section |^ for details . 



the elastic electron scattering off nuclei |115| since the wavefunctions oscillate rapidly. 



Even there the q^.^^ of eq.(|^.l) gave a semiquantitative description of the distortion. 
Quasielastic scattering is very smooth as function of position in the integrand, and one 
would expect the introduction of q^^^ should be adequate there. For a more recent 



discussion on this matter, see Triani 116 



Before the linear accelerator was shut down in Saclay, quasielastic experiments were 
performed on ^^C and ^°^Pb with a positron beam. In this case the sign of Vc changes, and 
the correction (0.1) to the momentum is of opposite sign. Using the effective momentum 
approximation of eq.(^l), it has been shown [117] that both electron and positron scat- 



tering of ^^C and ^osp^ could be well described in a consistent way. Introduction of q^^^ 
handles what is called the lowest order focusing effect in the theory. Whereas Onley et 



al ||114|| found higher-order focusing effects to be important, introduction of these would 
ruin the consistency between the electron and positron scattering. Of course without 
substantial numerical work, it is not possible to check the Onley et al's results. 

• Experimental Situation 

Although as theorists we cannot go into the experimental situation in detail, we can 
see that the early disagreement between the MIT-Bates data and Saclay data changed 
when MIT replaced their old scattering chamber, where background was corrected for by 
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simulation, by a new one. The updated MIT-Bates data and the Saclay data are now in 
substantial agreement. The analysis of the positron scattering introduced a new element. 
It could be straightforwardly explained by the effective momentum approximation (EMA). 



Kim et al ||118|| found that inclusion of higher order screening corrections in an approximate 
treatment of electron Coulomb distortion in quasi-elastic (e, e') reaction was important. 
Indeed for forward electron angles the low-energy side of the DWBA peak looks similar 
to the plane- wave result. The authors suggest that the phase factors effectively cancel the 
effect of the effective momentum on the low uj side of the quasielastic peak. However, if 
so, the next correction beyond the EMA would involve the square of the projectile charge 
which would be the same for electron and positron, negating the simple relation (^1). 
This would also mean that the difference between plane wave and DWBA for positrons 
would be large. In any case, these matters can be only conjectured because the Ohio 
University collaborations have not been extended to positrons although it seems obvious 
that this should be done. 

It is our belief that the whole situation, especially the theoretical analysis, has been 
a real shambles, for which we take partial responsibility. The Yennie, Boos and Ravenhall 
work was so convincing that theorists should have insisted that experiments were analyzed 
in the effective momentum approximation. 

• Theory with BR Scaling 

The mechanism for the change in quasielastic scattering can be quite simply un- 
derstood by considering the vector dominance model in which the electron couples to the 
nucleon through the vector meson. The vector mass enters into this coupling through 

Usually the numerator enters into the electron-p meson coupling and the denominator 
into the propagator, but the terms should be put together as in eq.(|^.2) because X>(my) 
must go to 1 as g in order to conserve the electron charge. If now my my ~ 0.8my 
then V{my) < V{my). 

The longitudinal response measures the charge density for which the operator is 

o' = \+ \n. m) 

The cu-meson couples to the | and the p-meson to the |r3. The scattering via the cu-meson 
exchange will be chiefly longitudinal in nature, since the isoscalar magnetic moment, which 
gives the transverse response, is small. Thus, the suppression in the longitudinal response 
is essentially 

Here what figures is the BR scaling in the uj meson mass. 
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For the transverse response the main operator is 

this response couphng chiefly to the isovector magnetic moment. Now in medium mj^ ^ 
m*^ as well as the propagator changing, so the relevant factor is 

/ 1 + q'^/m'l \ niN 



1 + /m*^ I m 



m 



N 



What enters here is the BR scaling in the nucleon mass. Now taking m*/mp ~ m*^/m]\[ 
according to the BR scaling, one finds that the two factors in (|^6) nearly cancel each other 
in the range of momentum transfers q ~ several hundred MeV in which the longitudinal 
and transverse responses have been measured. 

In Fig.^ we show comparison of longitudinal data with the theory of Soyeur, Brown 
and Rho. Also a point from Jourdan's analysis using the Onley et al theory is shown. 
The drop of the dash-dot hne, which included the Soyeur et al medium effects, below the 
Fabrocini and Fantoni solid line, which does not include these effects, is very clear, an 
~ 20 % effect. 

Both predicted and experimental effects on the transverse scattering are small. This 
is an example where BR scaling has a visible effect. Clearly the separation of longitudinal 
and transverse components was essential to expose BR scaling. 

• Comments 

The Soyeur-Brown-Rho treatment in [|1 1 1|| was more complicated than our above 
schematic model because it built in the quark structure of the nucleon at short distances, 
essentially through inclusion of the chiral bag model. Since the latter was not scaled, 
the effect of scaling in our schematic model was somewhat diminished. Nonetheless they 
were quite large as seen from the difference between the dashed line and the stared point 
"Jourdan analysis" and Soyeur et al's dash-dot line. 



It should be noted that Saito, Tsushima and Thomas ||1 19|| obtained within the 
quark-meson (QMC) model results for the quasi-elastic electron scattering similar to those 
of Soyeur, Brown and Rho. 

In Section |^, we will develop the evidence of medium-dependent vector-meson 
masses from the dilepton and photon experiments carried out at the SPS of CERN. 
Whereas the excess dilepton can be as simply produced in the Rapp-Wambach scenario 
as in the Brown-Rho scenario, the photons take us to higher densities and temperature, 
beyond the chiral restoration transition. We will develop a notion of complementarity 
between the BR dropping mass scenario and quark-gluon plasma (QGP) description for 
the photon treatment. As for quasielastic electron scattering we need not go into QGP 
since the hadron language is amply adequate. The dileptons we consider next fit in also 
quite simply with the quasielastic electron scattering. Our treatment in the next section 
will, however, allow us to unite the low-energy sector with the chiral restoration region of 
energies. 
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8 BR SCALING IN CHIRAL RESTORATION 

One can approach the phenomenology of the chiral restoration from a variety of 
different angles. We define ours by introducing the role of BR scaling by a simple con- 
struction of the chiral restoration transition as mean field in the Nambu-Jona-Lasinio 
model. More precisely we consider the transition as mean field up to both and T^, the 
critical density and temperature, although there may well be a density discontinuity at 
the transition, making it first order. Most likely the transition is second order for p = 0, 
with a tri-critical point on the phase boundary going towards finite density. What is 
important for our discussion, however, is that masses of mesons other than the pion go 
smoothly to zero in the chiral limit. The simplest possible model for this is mean field 
NJL. We follow the construction of Brown, Buballa and Rho ||120|| . Relevant parts of this 



paper were reviewed in Brown and Rho [0. The present description in NJL agrees with 
the vector manifestation scenario of Harada-Yamawaki's hidden local symmetry theory 
(see Section and with the quark-hadron complementarity picture given in Section |T^. 

8.1 Bag Constant and Scalar Field Energy 

Brown, Buballa and Rho showed that the Walecka nuclear mean field theory at 
nuclear matter density could be connected with NJL mean field theory at higher densities; 
e.g., the (T-field coupling the quarks in NJL can be taken to be 1/3 of that to nucleons in 
Walecka theory, 

9aQQ = 9aNN/^ - 10/3. (|.l) 

Furthermore, the scalar field energy ^m'^.a'^ plays the role of the bag constant B in the 
chiral restoration transition. On the quark-gluon side of the transition the appearance 
of the bag constant corresponds to the disappearance of the scalar field energy going 
upwards from below through the transition. The magnitude of the effective bag constant 
is only ~ 1/2 of that that would be given by the trace anomaly: 



i? = -^(0|(G^.f|0) (|.2) 



which to one loop (for Np = 3) is 

9a 



B = -^(0|(GM^|0) = (245MeVr 



Stt 

469 MeV/fm^ ~ 25^// (|.3) 



where Befj = |m^cr^. As found by Miller ||121|| , only about 50% of the glue is "uncovered" 
in the chiral restoration transition. The other 50% remains as "covered" glue. In Shuryak 



et al ||122|| , the soft glue disappears with the dynamically generated quark masses whereas 
the hard glue ("epoxy") is made up of instanton molecules which do not break chiral 
symmetry. The amount of trace anomaly which disappears at the transition gives -Be// — 
235MeV/fm=^. 
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8.2 "Nambu Scaling" in Temperature 

From the Gell-Mann-Oakes-Renner relation 



,2 



2m,{qq)/f^ (|.4) 



and taking not to scale with density, Brown and Rho [|T3] found hadron masses scaling 
as /* with /* oc \ Q It was then with some surprise that Koch and Brown [TI^] 

found that the temperature dependence of the entropy in the many-body system as cal- 
culated in lattice gauge calculations [|124|| could be reproduced in the hadron sector best 



with the scaling 

m* {qqY 



- , , (1-5) 

m {qq) 

for hadron masses other than the pion mass. In Fig.|l^ is shown how Nambu scaling 
behaves in the phase transition region. This scaling follows from the mean field NJL 
model and is referred to as "Nambu scaling" although it is a generic feature of the linear 
sigma model. The entropy is particularly suitable to study, since quark mean fields drop 
out so it tells us directly about the number of degrees of freedom. 

The hadron language is suitable for low temperatures. As stressed above and further 
elaborated later, we believe that our description is "dual" to a quark-gluon language as 
temperature moves towards the chiral restoration temperature. In particular, Koch and 
Brown limited their number of hadrons to 24, the number of light quarks. Rather than 
just cutting the number at this value, the same number can be consistently imposed by 



introducing excluded volume effects [|125|| . In any case the restriction that the hadronic 
description with dropping masses go over, in terms of the number of degrees of freedom, 
to the quark picture at higher temperature, must be built in. Otherwise it would not 
be possible to go beyond the Hagedorn temperature. The scale at which the change-over 
from hadron to quark language is made is one of convenience in accordance with the 
Cheshire Cat principle. Presumably it is in the region of the chiral restoration transition. 
This aspect becomes even more evident with the "hadron-quark continuity" discussed in 



Section 11 



Finally note that the Nambu scaling that fits the entropy in the hadronic description 
comes from mean field, i.e., the coupling constant QaQQ is not scaled. As long as we deal 
with the cr, the fluctuation field connected with chiral symmetry restoration, mean field 
seems to work very well, perhaps at least approximately up to the chiral restoration 
transition. In fact it was noticed in |79[ that the mean field description of nuclear matter 
in terms of chiral Lagrangians implemented with BR scaling required that the scalar 
coupling be unsealed at least up to nuclear matter density. This will also be observed in 



the problem of "strangeness equilibration" discussed in Section 10 



^''in [||, using Skyrme's Lagrangian and looking at pions in medium, the scaling /* cx \{qq)*\^^^ was 
obtained. We now know that because of chiral symmetry broken only slightly in nature, the in-medium 
properties of a pion cannot be given correctly in the same mean field argument as used for other non- 
Goldstone bosons. Therefore it is safe to say that there is no accurate description of the pionic property 
based on model considerations. Future experiments will help in this direction. 
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Figure 10: Comparison of entropy derived by Koch and Brown from the lattice data on energy 
density (dahsed hne) and from dropping masses (full lines) calculated in the Nambu scaling and 
in the naive Brown-Rho scaling. In the lower dashed line for (tpip) the bare quark mass has been 
taken out. 

8.3 "Nambu Scaling" in Density 

Given that there is nothing which dictates that the scaling in temperature and that 
in density be the same, the question is whether the Nambu scaling also works in density, 



at variance with those found in P, |T3l which are different from the Nambu one. To answer 
this question let us first address the question: How must the pion mass m* change with 
density in order to have Nambu scaling; i.e., hadron masses other than that of the pion 
to scale as 

m {qq) 
From Gell-Mann-Oakes-Renner equation 

f^=2m,{qq)/ml (|.7) 
and since the quark mass does not change with density, we have 
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Now, to linear order in density we have 

(qq) flm. 



2 



0.63 for p = po. (|.9) 



Setting (|^.6) equal to (^.7) at p = po and using (|.8), we find that in order for the Nambu 
scaling to apply to density, the pion mass should scale up with density by a factor 

ml _( {qq) \ ' _ ^ 26 for p = po. (g.lO) 



It turns out that this increase seems to be supported empirically although it is not clear 
whether the "measured" mass is the same quantity that appears in the in-medium GMOR 
relation. Experiments which pinned down the increase in pion mass from bare mass by 



a factor of 1.20 ~ 1.24 were performed recently at GSI ||126|| in the deeply bound states 
of ^''^Pb. Since the outer part of the pion wave functions is at the densities rather less 
than nuclear matter density, these factors may represent the lower limits and a factor of 
1.26 for p = Po seems quite reasonable. Thus the Nambu scaling may hold at least near 
nuclear matter density. 

We should mention however two caveats to this. One is that two-loop chiral per- 
turbation calculations find smaller values for the mass shift. The calculation by Park et 



al ||127|| predicts the pion mass at nuclear matter density to be at most 6 ~ 7 % higher than 



that at zero density. This is a calculation where the off-shell ambiguity associated with 
the pion in medium being an off-shell quantity is minimized. For a pion field that contains 
substantial off-shell ambiguity, one can easily obtain the mass factor 1.20 ~ 1.24 but this 
seems to be unreasonable. Another two-loop calculation by Kaiser and Weise ||128|| for 
the pion self-energy in asymmetric nuclear matter obtains about 10 % increase in mass 
of the pion which is consistent with the result of Park et al. The other caveat is that 
the scaling (||.6) with the value (|[9) at nuclear matter density would disagree with the 
value expected from the experimental value (§|.2) for Sgi, Eq.(^.4), for which we have 
<l>(po) = 0.78. 

In the above we have established that NJL theory at mean field [|120|| works well 



in describing the scaling of masses. In ref. it was shown that the Walecka mean field 
theory in terms of nucleons can be carried over with smooth change of parameters to a 
mean field theory of constituent quarks at higher densities. The constituent quark goes 



massless in the chiral limit with chiral restoration. In Section |11.2| , we will suggest how 
this can be understood in terms of color-isospin locking and hidden local symmetry. 

In low-energy nuclear physics we are used to introducing form factors with couplings. 
These decrease with increasing scale. In modern language, the form factors are there as a 
sort of regularization that plays an important role in effective field theories as discussed in 
Section ^. Why do we do so well here without form factors ? We believe that the reason 
lies in the special role that the scalar field a has in chiral restoration. A priori it would 
seem reasonable to continue the linear approximation for cr-exchange between quarks 

Qiq'-m^lr'Q (|.ll) 
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where Q is the constituent quark, down to small m* where higher-order or nonlinear 
effects must enter. We believe that these effects must be counterbalanced, to a large 
extent, by form factors. As we develop in the next section, this does not seem to be 
so with vector mesons, where form factors really do not cut down the exchanges. As 
formulated more precisely in Section |11.2| , the vector meson exchange at low densities 
(i.e., flavor symmetry) changes into gluon exchange at higher density (QCD symmetry) 
and asymptotic freedom assures that their effects go to zero with increasing scale. 



9 SIGNALS FROM HEAVY-ION COLLISIONS 

Extreme conditions of high temperature and/or high density are expected to be 
created in relativistic heavy-ion collisions being probed in RHIC at Brookhaven and to 
be probed in CERN's LHC In this section, we discuss how BR scaling manifests itself 
in heavy-ion processes, leaving high-density zero-temperature situation to later sections. 
Before going into the details, we first define precisely the meaning of BR scaling in the 
context that we shall use in this section. 



9.1 Top-Down and Bottom-Up and How They Connect 

There are effectively two ways of approaching BR scaling. The BR scaling as first 
suggested in p| is a "top-down" version starting, using quark mean field, from a high 
temperature/density regime wherein constituent quarks (or quasiquarks) are relevant de- 
grees of freedom, i.e., near the chiral transition point and then extrapolating down to 
a low temperature/density regime wherein hadrons are appropriate degrees of freedom. 
In contrast, Rapp and Wambach and others ||129|| calculated medium dependent hadron 



properties in hadronic variables in a strong-coupling perturbation approach. We shall call 
this a "bottom-up" approach. The advantage of the latter approach is that one can resort 
to standard nuclear physics information using phenomenological Lagrangians at tree order 
with more or less known coupling constants. This approach has been found to be success- 
ful in reproducing the medium- dependent p-meson properties deduced from the CERN 
CERES dilepton experiments. In these experiments, most of the p-mesons come from 
densities less than nuclear matter density. In the low-density regime, the phenomeno- 
logical approach is justified by construction in a well-prescribed way. The disadvantage 
of this approach is that as density increases beyond nuclear matter density, the possible 
vacuum change makes the tree order calculation unreliable because of the strong coupling 
or more precisely a large anomalous dimension which signals that one is fluctuating at a 
wrong ground state and that it is preferable to shift to a different vacuum even though 
no phase transitions are involved. This aspect was associated above with differences in 
topological vacua. 

In ||130| , |131| , |132|| , the top-down BR scahng and the bottom- up Rapp- Wambach 



approach were combined into a unified description. We will follow this description. We 
consider the most important component of the Rapp- Wambach theory, i.e., the introduc- 
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tion of the "p-sobar" which is the excitation of the / = 1, J = 3/2 iV*(1520) coupled 
together with the nucleon-hole to the quantum numbers of the p. There are of course 
other A^*-hole states of the same quantum numbers which would come in but here we shall 
focus on the dominant component only. At zero density this isobar-hole state is at 580 
MeV. The phenomenological Lagrangian coupling the "elementary p" and the isobar-hole 
p that we shall use is 

= ^V'jv* {%s- Pa - pis- q)T,i;r, + h.c. , 

TTlp 

where /^^*^/47r ~ 5.5 from fits to photodisintegration data. Here (g, go) are the p-meson 
momenta, we shall keep only the go term in our dilepton discussion. The two branches 
of the p spectral function can be located by solving self-consistently the real part of the 
p-meson dispersion relation (taking g = 0) 

9o = "^p + ReSpAr*Ar(go). (||.2) 

The natural width of the iV*(1520) is supplemented at p = po by an additional medium- 
dependent width of 250 MeV. 

Including also the backward-going graph, the iV*(1520)iV~^ excitation contributes 
to the self-energy of the p-meson at nuclear matter density no (in this section we denote 
density by n to avoid confusion with the p-meson)) 

where AE ^ 1520 - 940 = 580 MeV and Ttot = To + Tmed is the total width of the 
A^*(1520). If we neglect any in-medium correction to the width of the A^*(1520) we 
find two solutions, located at Qq ~ 540 MeV and go" ~ 895 MeV. As n increases, S 
will become larger and initially the go of the lower branch will go down. However, this 
downward movement will be halted once go has decreased substantially. We note that 
already at = 0, the A^*(1520) decays ~ 15 % of the time to A^p indicating that there is 
already mixing there. The mixing becomes stronger with increasing n. 

As noted, the Rapp-Wambach theory works well for the low densities n <no impor- 
tant for the CERES experiments. We show in Fig.|Tl| fits to the CERES spectra with both 
Rapp-Wambach and Brown-Rho theories. In the BR approach, however, the vacuum is 
"sliding," that is, for each new density, one has a new "vacuum" to fluctuate around. 
Recalling the skyrmion structure for given baryon number, we observe that "vacua" are 
solitonic in nature and cannot be connected by perturbation theory. Thus, the vacuum 
at n 7^ has no continuous connection with that at n = 0. This means that with the La- 
grangian (^1), the Rapp-Wambach theory can never approach the BR theory as density 
increases. One can see from (^.3) that the self-energy will stop bringing down the energy 
of the lower branch of the p as (go/w^p)^ becomes smaller. This situation is somewhat like 



35 



The non-covariant form of the Lagrangian is not unique as pointed out in |133] 
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Figure 11: Dilepton data from central Pb (158 GeV/u)+Au collisions compared with 
Rapp/Wambach (labelled as "in- medium tttt") and Brown/Rho (labelled as "dropping m^") 
for qf^ < 0.5 GeV (upper) and qf^ > 0.5 GeV (lower). 



the usual Walecka mean field theory where the nucleon self-energy has a scalar density 
Ps = V'WoV' ill it- As ps goes to zero, the self-energy which involves ps linearly does 
likewise, so the nucleon mass can go to zero only at n = oo. However with the introduc- 
tion of negative energy states in Nambu-Jona-Lasinio theory, the mass of the constituent 
quark goes smoothly to zero at a density ric ~ Suq ||120|| . As argued in ||131| , |132|| , this 



can happen with (^.3) only if we replace there the factor (go/'^p)^ by ~ e{qQ/m*Y. The 
unknown factor e can be adjusted to give the same critical density as in NJL model. It 
is plausible that this change-over is related to the requirement in the Harada-Yamawaki 
scheme which is not present in few-order perturbation theory that the phase transition co- 
incide with Georgi's vector-limit fixed point. As mentioned, the Harada-Sasaki work [36] 
in temperature indeed suggests that the Wilsonian matching with QCD with hidden gauge 
invariance requires the parametric p mass to have an "intrinsic" T/n dependence. Once 
one accepts this, what follows turns out to be quite reasonable: with e = 1, m* drops to 
zero at the same density as in the NJL model, Uc ~ 3?7,o. 



9.2 Distinguishing Rapp/Wambach and Brown/Rho 

The modified sobar picture obtained by replacing mp by m* as described above still 
differs in detail from the BR scaling picture of 1991 even though they describe more or 
less equally well the m* — > limit, because at any density less than ric there are two 
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branches of the p quantum number. However a,s n ^ Uc the width of the lower branch 
goes to zero and all of the spectral strength goes into it so we do recover the original BR 



scaling. We will encounter this feature once more in Section |11.2| from the consideration 
of hidden local symmetry and color-flavor-locking in QCD. 

At densities n < uq, the lower p branch most likely becomes very wide because of 
the large number of open channels that it couples to, which has led Rapp and Wambach 
to interpret their fit to the dilepton data from the broadening of the p, the lower parts of 
the spectral function contributing importantly to the dilepton production because of the 
greater Boltzmann factors. 

So how do we distinguish between Rapp/Wambach and Brown/Rho; equivalently 
between broadening and movement towards chiral restoration? Formally we do not violate 
gauge invariance or any other invariance by using m* rather than mp. It is really up to 
nature to decide. 

In a recent paper, Alam, Sarkar, Hatsuda, Nayak and Sinha [ p.34j| have suggested 
that nature has decided for Brown/Rho 0. This work is extended in Alam et al ||137 |. 



In these papers, the calculations were based on mesonic processes such as np n^y, 
TTTT p'-f, etc. leading to thermal production of photons. The observation then was 
that when increasing the width of the p meson to ~ 1 GeV in order to take into account 
medium effects, the photon yield decreases as compared to the case using the free p width, 
thus underestimating the spectra. However these calculations neglect the contributions 
to the photons coming from some of the processes that increase the width, e.g., the 



^^In fact Alam et al suggest two scenarios that can describe the photons, one being the hadronic 
description with BR scaling and the other one in which quark-gluon-plasma (QGP) is initially produced 
at a temperature Ti — 196 MeV. One has a mixed phase of QGP and hadronic matter which persists 



down to a transition to hadronic matter at Tc = 160 MeV. In a later paper [135|, Sarkar et al find that for 
relatively central collisions, the dilepton yield can be described by either of the two scenarios, BR scaling 
or one in which the QGP is originally produced as above. These authors find that "as of yet, it has 
not been possible to explain the observed low-mass enhancement of dileptons measured in the Pb+Au 
collisions as well as in the S+Au collisions at the CERN SPS in a scenario which does not incorporate 
in-medium effects in the vector meson mass." 

As noted earlier, we believe that two descriptions, BR scaling and the one in which the QGP is initially 
produced in a mixed phase, to be "dual" as long as the number of degrees of freedom in the former is 
constrained to equal that in the latter. The mixed phase in the scenario with the QGP formation could, 
however, be just the part of a second-order phase transition in the BR scaling scenario. In the latter the 
energy is used up, going down in temperature through the phase transition, in building up the scalar 
field energy as the hadrons gain back their masses. Most of the entropy decrease will come from the 
higher mass hadrons disappearing as their mass is increased, due to decreasing Boltzmann factors, rather 
than through expansion. It seems therefore reasonable that the various freeze-out species exhibit thermal 
equilibration at the chiral restoration temperature, the one at which hadron masses go to zero in the 
chiral limit in the BR scaling scenario, since the system spends a long time there. 

Although the dropping masses with BR scaling and production of QGP in mixed phase at T ^ 200 MeV 
may be dual around the phase transition, it is excluded that the latter scenario describe the dropping 
masses found in ~ 300 MeV inelastic electron scattering, which we described in the previous section, or 
the many indirect evidences we described earlier in our review. In other words, the duality is useful only 
at higher scales, near and above the chiral restoration transition. One should remember the Beg-Shei 
theorem, that the nature of symmetry realization (i.e., Wigner-Weyl vs. Nambu-Goldstone) is irrelevant 



to discussing the short-distance symmetry [136 



64 



p — s> A^*(1520)A^~"^ excitation (the collective "p-sobar") which provides the lower-mass p 
excitation in the Rapp/Wambach calculation. This contributes to the photon through 
the A^*(1520) N'-f decay. Calculations in progress (by C. Gale, M.A. Hadasz and R. 
Rapp) show that these baryonic processes contribute importantly so the conclusions of 
||134| , |137|| may be premature. 

The bottom-up approach of Rapp and Wambach should be reliable at low density 
since it is based on phenomenological Lagrangians constructed in standard nuclear physics 
with known properties. Our principal argument is that this picture should cede to the 
Brown-Rho picture at higher densities replacing rup by m* in scaling the Lagrangian 
required to be consistent with Harada-Yamawaki's vector manifestation. Of course much 
more detailed experiments will be needed to confirm this scenario. 



10 "BROAD-BAND EQUILIBRATION" OF 

STRANGENESS IN HEAVY-ION COLLISIONS 

10.1 Kaons and Chiral Symmetry 

The strange quark is massive compared with the light (u, d) quarks, with nis ~ 150 
MeV. This gives the kaon a mass of ~ 1/2 GeV. This is a manifestation of both explicit and 
spontaneous chiral symmetry breaking. Unlike the previous cases where the decondensing 
of the condensate (qq) - that is, restoring the spontaneously broken chiral symmetry - 
enters directly into the phenomena we were looking at, here it is the "rotating away" of 
the explicit chiral symmetry breaking intricately tied in with the decondensing of (qq) 
that is principally manifested in the behavior of kaons in dense medium. This section 
deals with how this aspect is probed in the kaon production experiments at GSI. Based 
on work done recently | ]138| |, we will develop, in conjunction with BR scaling, the notion 
of "broadband equilibration" in heavy-ion processes and suggest the vector decoupling in 
dense medium. 



10.2 Equilibration vs. Dropping Kaon Mass 

Following work by Hagedorn [|139|| on production of anti-^He, Cleymans et al [|140|] 
have shown that for low temperatures, such as found in systems produced at GSI, strangeness 
production is strongly suppressed. The abundance of mesons, in systems assumed to 



be equilibrated, is given by ||141 



Here the ^^'s are the degeneracies. Because strangeness must be conserved in the inter- 
action volume V, assumed to be that of the equilibrated system for each which is 
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produced, a particle of "negative strangeness" Q containing s, say, K or A, must also be 
produced, bringing in the ^ or A phase space and Boltzmann factors. The production 
is very small at GSI energies because of the low temperatures which give small Boltzmann 
factors for the K and A in addition to the small Boltzmann factor for the K~^. Note the 
linear dependence on interaction volume which follows from the necessity to include K or 
A phase space. 

In an extensive and careful analysis, Cleymans, Oeschler and Redlich [|141|] show that 
measured particle multiplicity ratios n~^/p, vr'/vr"'", d/p, K~^/n~^, and / K~ - but not 
77/7?° - in central Au-Au and Ni-Ni collisions at (0.8-2.0)A GeV are explained in terms of a 
thermal model with a common freeze-out temperature and chemical potential, if collective 
flow is included in the description. In other words, a scenario in which the kaons and anti- 
kaons are equilibrated appears to work well. This result is puzzling in view of a recent 



study by Bratkovskaya et al [142] that shows that the mesons in the energy range 
considered would take a time of ~ 40 fm/c to equilibrate. We remark that this is roughly 
consistent with the estimate for higher energies in the classic paper by Koch, Miiller 
and Rafelski | J143| | that strangeness equilibration should take ~ 80 fm/c. Such estimates 
have been applied at CERN energies and the fact that emergent particle abundances are 



described by Boltzmann factors with a common temperature ~ 165 MeV ||144|| has been 
used as part of an argument that the quark/gluon plasma has been observed. 

We interpret the result of [|141|] as follows. Since free-space masses are used for 



the hadrons involved, Cleymans et al [|141|| are forced to employ a /i^ substantially less 



than the nucleon mass m^v in order to cut down A production as compared with K~ 
production, the sum of the two being equal to production. This brings them to a 
diffuse system with density of only ~ po/4 at chemical freeze-out. But this is much too 
low a density for equilibration. 

We shall first show how this situation can be improved by replacing the K~ mass 
by the i^'~energy at rest a;^_ = uj-{k = {]) < rriK- (The explicit formula for uj± is given 
later, see eq. (^.25).) In doing this, we first have to reproduce the to K~ ratio found 



in the Ni -|- Ni experiments |145 



nx+/^x--30. (|T0|.2) 

Cleymans et al reproduce the earlier smaller ratio of 21 ± 9 with fiB = 750 MeV and 
T = 70 MeV. How this or rather (p!0[2) comes out is easy to see. The ratio of the second 
term on the RHS of eq.(|TD|.l) to the first term is roughly the ratio of the exponential 
factors multiplied by the phase space volume 

9k \PkJ e ^K-/^ yniKj e ""k-/^ ^ 

where we have used g{^^ gj^ — 2, {pY /2m ~ |T, = 1115 MeV and E^^ = 495 MeV. 
Inclusion of the E and S hyperons would roughly increase this number by 50% with the 

^^By "negative strangeness" we are referring to the negatively charged strange quark flavor. The 
positively charged anti-strange quark will be referred to as "positive strangeness." 
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result that the ratio of K to A, E, S production is 

Uk- 



1/32. (|T0|.4) 



Since a must be produced to accompany each particle of one unit of strangeness (to 
conserve strangeness flavor), we then have 

nK+/nx-~33. ([10|.5) 

This is consistent with the empirical ratio ([T0|.2). It should be noted that had we set hb 
equal to m^v, we would have had the to ratio to be ~ 280 because it costs so much 
less energy to make a A (or S) rather than K~ in this case. In other words the chemical 
potential /i^ is forced to lie well below m^v in order to penalize the hyperon production 
relative to that of the K~^s. 



One can see from Fig. 5 in Li and Brown ||146|| that without medium effects in the 
K~ mass, the /K~ ratio is ~ 100, whereas the medium effect decreases the ratio to 
about 23. This suggests how to correctly redo the Cleymans et al's analysis, namely, by 
introducing the dropping K~ mass into it. 

In the next subsection we show that positive strangeness production takes place 
chiefly at densities greater than 2po- As the fireball expands to lower densities the amount 
of positive strangeness remains roughly constant. The number of ii"+'s is such as to be in 
equilibrium ratio /ti^ with the equilibrated number density of pions at T = 70 MeV, 
'^a 0.37Tj^g7 fm~^. Only in this sense do the ii'+'s equilibrate. It will be noted, however, 
that with the T of 70 MeV, /is is chosen so that the empirical number of hyperons are 
produced. Since the number of i^'+is one greater than that of the hyperons (including the 
i^~in the negative strangeness), it will also be the apparent equilibrated ratio for this 
and T. Thus, with these two latter parameters the vr, i^'+and K~caii be put into apparent 
equilibration. It is fairly easy for the nucleons to equilibrate with the pions at the given 
temperature because of the strong interaction. 

It is amusing to note that the "equilibrated ratio" of ~ 30 for the nK+/nK- holds 
over a large range of densities for T = 70 MeV, once density-dependent masses are 
introduced, in that the ratio R of ([l3.3) is insensitive to density. (Remember that because 
of the small number of ii"~'s, the number of i^+'s must be nearly equal to the number 
of hyperons. A, S~ and S, in order to conserve strangeness.) This insensitivity results 
because a;^- decreases with density at roughly the same rate as /i^ increases. We can 
write R of (|l3-3), neglecting possible changes in T and in our lowest approximation, 

^^In order to reproduce this result with fiB = 750 MeV and T — 70 Mev within our approximation, 
we have assumed only the E~ and S° hyperons to equilibrate with the A. This may be correct because 
the E+ and S couple more weakly. Inclusion of the latter could change our result slightly. Probably 
they should be included in analysis of the AGS experiments at higher energies where they would be more 
copiously produced. 
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as 



R=(!I!±) e-^-'^. (0.6) 



As will be further stressed later, the most important point in our arguments is that the 
fiB + ^^K- nearly constant with density. This is because whereas fiB increases from 860 
MeV to 905 MeV as p goes from 1.2po to 2.1po, (^k- decreases from 380 MeV to 332 MeV, 
the sum fiB + decreasing very slightly from 1240 MeV to 1237 MeV. Indeed, even 
at p = po/4, Pb + m*j^- ~ 1218 MeV, not much smaller. 

We believe that the temperature will change only little in the region of dropping 
masses because in a consistent evolution (which we do not carry out here) the scalar field 
energy in a mean field theory plays the role of an effective bag constant. In 

ref . ||120|] (which was reviewed in Section |) this is phrased in terms of a modified Walecka 
theory, 

Beff = \mla^ ^mliM^/g^NNf, (0.7) 

the a going to M^/gaNN as the nucleon effective mass goes to zero with chiral restora- 
tion. Most of the energy with compression to higher densities goes with this effective 
bag constant, estimated [120 to be ~ 240 MeV/fm^, rather than heat, mocking up the 
behavior of a mixed phase with constant temperature. Moreover at p = 2po where 
may be ~ O.Sm^v, only about 25% of the bag constant B may have been achieved, so there 
may be some increase in temperature. We shall, at the same level of accuracy, have to 
replace in the prefactor of eq.(|TD|.6) by cu^-. We adjust the increase in temperature 
so that it exactly compensates for the decrease in prefactor so that the / K'Taiio is 
kept the same, as required by experiment | fL45| | . We then find that the temperature at 



p = 2po must be increased from 70 to 95 MeV. This is roughly the change given by that 
in inverse slopes of i^~and i^'+transverse momentum distributions found in going from 



low multiplicities to high multiplicities |p.45|| . 

In any case, we see that R will depend but little on density. This near cancellation 
of changes in the factors is fortunate because the K~ + p ^ A reaction, operating in the 
negative strangeness sector, is much stronger than the positive strangeness reactions, so 
the former should equilibrate to densities well below 2po and we can see that the "apparent 
equilibration" might extend all the way down to p ~ po/4. 

The near constancy of R with density also explains the fact the K'^ / K^iaXio does 
not vary with centrality 0. Although R is the ratio of A's to K~^s, both of which are in 
the negative strangeness sector, nonetheless, the number of i^'+'s must be equal to the 
sum of the two and since the A's are much more abundant than the K~^s, R essentially 
represents the / K'laiio. 

^^To be fully consistent, we would also have to consider the medium modifications of the A and 
_ftr+properties. These and other improvements are left out for our rough calculation. 
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We are grateful to Helmut Oeschler for pointing this out to us. 



68 



K , 1.8 GeV 






Ni + Ni, central 






produc. 






reseat. 













12 3 



Figure 12: Calculations by Bratkovskaya and Gassing (private communication) which show the 
density of origin and that of the last interaction of the K~ mesons. 



Detailed transport results of Bratkovskaya and Gassing (see Fig.[T^) show the last 
scattering of the detected K~to be spread over all densities from po/2 to 3po, somewhat 
more of the last scatterings to come from the higher density. This seems difficult to 
reconcile with a scenario of the numbers being decided at one definite density and 
temperature, but given our picture of dropping masses, one can see that the K'^ / K~Ta.tio 
depends little on density, that is, on fis at which the i^^'last interacts. In any case 
we understand from our above argument that the apparent density of equilibration can 
be chosen to be very low in a thermal description and still get more or less correct 
K+/K-mtio. 



10.3 The Equilibrium K+/7T+ Ratio 

We show here how our argument that gives a correct K^/K^ ratio can reproduce 
the K^/n^ ratio. Let us leave T = 70 MeV and choose p ~ 2po as educated guesses. 
We are thereby increasing the equilibration density by a factor ~ 8. We then calculate 
the baryon density for this fiB and T = 70 MeV and find p = 2po which checks the 
consistency. 

According to Brown et al | |147| ] the production under these conditions will come 
chiefiy from BB NAK, with excited baryon states giving most of the production. From 
the sohd curve for p/po = 2 in fig. 9 of ||147|| we find {av) ~ 2 x 10^^ mb = 2 x lO"'' fm^. 
The rate equation reads 

6^K = l{cTEVvBB)nl ~ ^(2 X 10'^ fm2)p| = dn/dt = ^ x 10"^ fm"^ (^.8) 

where B, Y and K stand, respectively, for baryon, hyperon and kaon. For this, we have 
taken {(t^^vbb) from fig.9 of ||147|| and p = 2po = | fm~^. Choosing a time t = 10 fm/c 
we obtain 

nK+ ^ 5m K ^ = ^ X 10"^ fm-^ (|10|.9) 
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Now equilibrated pions have a density 



= 0.37(T/197MeV)^ fm"^ = 0.016 fm"^ 
for T = 70 MeV. From (|10[9) and dT^-lO) we get 

nK+/n^+ ^ 0.0069 



(0.10) 



(0.11 



which is shghtly below the "equilibrated value" of 0.0084 of Table 1 of Cleymans et al ||141|| . 
Production of K^hj pions may increase our number by ~ 25% ||148|| . 

Our discussion of i^"*" pro duct ion is in general agreement with earlier works by Ko 
and collaborators [|148| , |149|] . In fact, if applied at the quark level, the vector mean field 
is conserved through the production process in heavy-ion collisions, so it affects only the 
strangeness condensation in which there is time for strangeness non-conservation. These 
earlier works establish that at the GSI energies the i^+content remains roughly constant 
once the fireball has expanded to ~ 2po, so that in this sense one can consider this as a 
chemical freeze-out density. 

It should be noted that in the papers ||150| , |147|, |148| , |l 
potential - scalar plus vector - on the i^^-meson is slightly attractive at p 



and others, the net 
2po even 

though the repulsive vector interaction is not decoupled. A hint for such change-over was 
noted already at nuclear matter density by Friedman, Gal and Mares ||151[ . Since in our 
top-down description the vector interaction can be thought of as applied to the quark 
(matter) field in the K~^, the total vector field on the initial components of a collision is 
then the same as on the final ones, so the vector mean fields have effectively no effect on 
the threshold energy of that process. Our proposal here is that the vector mean field on 
the i^"*" should be below the values used by the workers in ||150| , p.47| , |148| , |149|| due to the 
decoupling. However, in comparison with ||150|| , our total potential on the i^"'"at p = 2po 
is ~ —85 MeV, as compared with their ~ —30 MeV. We have not redone the calculation 



of [147] to take into account this difference. 



10.4 The Top-Down Scenario of K Production 

Brown and Rho discussed fluctuations in the kaon sector in terms of a simple 
Lagrangian 

6Ckn = ^(NioN)KdtK + -^(NN)KK +... = £^ + + • (0.12) 

oj / 

It was suggested there that at high densities, the constituent quark or quasi-quark de- 
scription can be used with the cu-meson coupling to the kaon viewed as a matter field 
(rather than as a Goldstone boson). Such a description suggests that the u coupling to 
the kaon which has one non-strange quark is 1/3 of the uj coupling to the nucleon which 
has three non-strange quarks. The in the Lagrangian was obtained by integrating out 
the cu-meson as in the baryon sector. We may therefore replace it by the interaction 

Vk± ~ ±^Vn. (0.13) 
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Figure 13: Projection onto the a, K plane. The angular variable 9 represents fluctuation toward 
kaon mean field. 



In isospin asymmetric matter, we shall have to include also the p-meson exchange 
with the vector-meson coupling treated in the top-down approach. 

For the top-down scenario, we should replace the chiral Lagrangian ([T0|.12) by one 
in which the "heavy" degrees of freedom figure explicitly. This means that -tj}^ in the first 



2/* 

in the second term by |mx 



term of (|T0|.12) should be replaced by g*'^/m*J and ^ff 
assuming that both uj and a are still massive. We will argue below that while the uj mass 
drops, the ratio g*"^ jm^J' should stay constant or more likely decrease with density and 
that beyond certain density above nuclear matter, the vector fields should decouple. On 
the other hand, is not scaled in the mean field that we are working with; the motivation 
for this is given in Brown, Buballa and Rho | 120 | who construct the chiral restoration 
transition in the mean field in the Nambu-Jona-Lasinio model. Thus 



3 ""mf 



(0.14) 



In this framework, m* is the order parameter for chiral restoration which drops a la BR 
scaling P|: 



(0.15) 



with y ~ 0.28, at least for p < Po H- 

Once the vector is decoupled, a simple way to calculate the in-medium kaon effective 
mass, equivalent to using the is to consider the kaon as fluctuation about the "a"- 



axis in the V-spin formalism [[L52|] as depicted pictorially in Fig.|13[ The Hamiltonian for 
explicit chiral symmetry breaking is 



H 



xSB = ^KN{NN)cos{e) + -mlr\in\e) 
^ S;,^(iViV)(l-^) + im^rV 



(0.16) 



may well be different from this value for p > pQ. In fact the denominator of ^{p) could even be 
significantly modified from this linear form. At present there is no way to calculate this quantity from 
first principles. 
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where the last expression is obtained for small fluctuation 6. Dropping the term indepen- 
dent of 6', we find 



K 



Using eq. (0.14) we obtain 



3 m* rriK 

In accord with Brown and Rho ||83| we are proposing that eq.([l0|.17) should be used for 
low densities, in the Goldstone description of the K^, and that we should switch over to 
eq.(0.18) for higher densities. It is possible that the mx appearing in (0.18) should be 
replaced by for self-consistency but the dropping of m*^ makes the of (p^.l8) 
decrease more rapidly than that of (0.17) so that eq.(0.17) with {NN) set equal to the 
vector density p a much used formula valid to linear order in density 



(0.19) 



K , 



obviously gives too slow a decrease of with density. 

Although the above are our chief points, there are two further points to remark. 
One, even without scaling, our vector interaction on the kaon is still too large. Two, more 
importantly, there is reason to believe in the large Tj^n term, 

T.KN ~ 400 MeV. (0-20) 

This comes from scaling of the pion sigma term 

_ {iTLu + Tns){N\uu + ss\N) 

= 7 \ \IAT\- 1 AA\AT\ ^-^^- (0-21) 

[mu + md){N\uu + dd\N) 

Taking m, ~ 150 MeV, m„ + md ~ 12 MeV, S^^v = 46 MeV and (Ar|ss|A^) ~ \{N\dd\N) 
from lattice calculations ||153|| , one arrives at (p^.20). 



Other authors, in adjusting the E term to fit the kaon-nucleon scattering amplitudes, 
have obtained a somewhat smaller T^kn- This can be understood in the chiral perturbation 



calculation of C.-H. Lee ||154|| where the only significant effect of higher chiral order terms 
can be summarized in the "range term''^; namely T^kn is to be replaced by an effective 



'KNjeff 



0.37uj*j^/ml)J:KN- (0-22 



^^The correction to this approximation which may become important as the nucleon mass drops comes 
as a "1/m" correction in the heavy-fcrmion chiral perturbation theory (as for the "range term" mentioned 
below) and can be taken into account systematically. It can even be treated fully relativistically using a 
special regularization scheme being developed in the field. Our approximation does not warrant the full 
account of such terms, so we will not include this correction here. 

^•^As mentioned, in the language of heavy-baryon chiral perturbation theory, this corresponds to a 
"1/m" correction term. 
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It should be pointed out that although the T^kn is important at low densities, ujk decreases 
with m^, this "range-term" correction becomes less important at higher densities. This 
effect - which is easy to implement - is included in the realistic calculations. 

10.5 Partial Decoupling of the Vector Interaction 

10.5.1 Theoretical and experimental evidences 

In Section ^, a general argument for the vector decoupling in hot/dense matter 
was presented based on Wilsonian matching with QCD of HLS. Here we supply more 
specific theoretical and empirical reasons why we believe that the vector interaction should 
decouple at high density. 

1. We first give the theoretical arguments why the vector coupling g* should drop with 
density. 



The first is the observation by Song et al |[L55|| that describing nuclear matter 
in terms of chiral Lagrangian in the mean field requires the ratio g'^/m*^ to 
at least be roughly constant or even decreasing as a function of density. In 
fact to quantitatively account for non-linear terms in a mean-field effective 
Lagrangian, a dropping ratio is definitely favored 0. For instance, as discussed 
in |155|, the in-medium behavior of the cu-meson field is encoded in the "FTSl" 



version of the non-linear theories of ref. |p,56|| . In fact, because of the attractive 



quartic uj term in the FTSl theory, the authors of [ |156|] have (for the parameter 
?7 = — 1/2 favored by experiments) g^j^ /rn*J^ ~ Q-^gJ^ /rnj^ as modification of 
the quadratic term when rewritten in our notation . In other words, their 
vector mean field contains a partial decoupling already at p po although 
they do not explicitly scale g^ as we do. 

Historically, Walecka-type mean field theories with only quadratic interactions 
(i.e., linear Walecka model) gave compression modulus K ~ 500 MeV, about 
double the empirical value. This is cured in nonlinear effective field theories like 
FTSl by higher- dimension non-renormalizable terms which effectively decrease 



the growth in repulsion in density. As suggested in ||155|| , an effective chiral 



Lagrangian with BR scaling can do the same (by the increase in magnitude of 
the effective scalar field with density) but more economically and efficiently. 

This decrease of (y'^/m* as density increases toward nuclear matter density 
contrasts with the increase of g^p^^ I ""^p within the same density range discussed 
in Section |6TT| . This can be explained as follows. In free space g^i = gu)NN = 
^QpNN whereas near chiral restoration, we expect (as in the case of the quark 
number susceptibility mentioned below) that 5'*jvAr — 9pNN- This implies that 
at low density, the u coupling must fall faster than the p coupling by "shedding 
off' the factor of 3. 



Since the non-linear terms - though treated in the mean field - are fluctuation efi^ects in the eff'ective 
field theory approach, this represents a quantum correction to the BR scaling. 
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Implementing baryons into the HLS Lagrangian, Kim and Lee ||157|| showed 



exphcitly that both gy and my (where V stands for hidden gauge bosons) are 
predicted to fall very rapidly with baryon chemical potential ||157|| The main 
agent for this behavior is found in ||157|| to be the pionic one-loop contribution 



linked to chiral symmetry which is lacking in the mean-field treatment for BR 
scaling 0. The Kim-Lee arguments were made for the p which has a simple 
interpretation in terms of hidden gauge symmetry but it will apply to the u if 
the nonet symmetry is assumed to hold in nuclear medium. 

Finally, close to chiral restoration in temperature, there is clear evidence from 
QCD for an equally rapid drop, specifically, from the quark number susceptibil- 



ity that can be measured on the lattice |13]. The lattice calculation of the quark 



number susceptibility dealt with quarks and the large drop in the (isoscalar) 
vector mean field was found to be due chiefly to the change-over from hadrons 
to quarks as the chiral restoration temperature is approached from below. The 
factor of 9 in the ratio gljNN I Q^qq (where Q is the constituent quark) should 
disappear in the change-over. Now since the electroweak properties of a con- 
stituent quark (quasiquark) are expected to be the same as those of a bare 
Dirac particle with gA = 1 and no anomalous magnetic moment (i.e., the QCD 
quark) ||158|| with possible corrections that are suppressed as l/N^ [ [159|| , there 



will be continuity between before and after the chiral transition. This is very 



much in accordance with the "Cheshire-Cat picture" developed elsewhere |TT 
In fact, it is possible to give a dynamical (hadronic) interpretation of the above 
scenario. For instance in the picture of [ [I32|| sketched in Section ^ this may be 
understood as the "elementary" uj strength moving downwards into the "nu- 
clear" u;, the [N*{l'b2Q)N~^Y^^'^^^ isobar-hole state involving a single-quark 



spin flip |p,33|| . The mechanism being intrinsic in the change-over of the degrees 
of freedom, we expect the same phenomenon to hold in density as well as in 
temperature. The upshot of this line of argument is that the suppression of 
the vector coupling is inevitable as density approaches the critical density for 
chiral transition. In Section |11.2| ,we will have an additional support for this 
from color-isospin locking in QCD both in low density and in high density. 

We believe that the different behavior of vector and scalar mean fields, the latter to 
be discussed below, follows from their different roles in QCD. With the vector this 
is made clearer in the lattice calculations of the quark number susceptibility which 
involves the vector interactions. In Brown and Rho [|13|, it is shown that as the 
description changes from hadronic to quark/gluonic at T ~ Tc, the critical temper- 
ature for chiral restoration, the vector interaction drops by an order of magnitude, 
much faster than the logarithmic decrease due to asymptotic freedom. We expect 
a similar feature in density, somewhat like in the renormalization-group analysis 



^^Kim and Lee found that even the ratio ffy/my feh rapidly. This is at variance with what one would 
expect if the vector manifestation a la Harada-Yamawaki were realized. As noted, this may be due to 



the incompleteness of the renormalization group analysis of |157| 
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for the isovector vector meson p of Kim and Lee ||157||. The scalar interaction, on 



the other hand, brings about chiral restoration and must become more and more 
important with increasing density as the phase transition point is approached. 

2. From the empirical side, the most direct indication of the decoupling of the vector 



interaction is from the baryon flow ||160|| which is particularly sensitive to the vector 



interaction. The authors of [|160|| find a form factor of the form 

MP) = (0-23) 

with Ay =0.9 GeV is required to understand the baryon fiow. Connecting momenta 
with distances, one finds that this represents a cutoff at 

\/6 

i?cnto// ~ -X- ~ 0.5 fm. (0.24) 

Ay 

Furthermore it is well known that the vector mean field of the Walecka model must 
be modified, its increase as E/m^ removed, at a scale of ~ 1 GeV. The reason for 
this is presumably that inside of i? ~ 0.5 fm, the finite size of the solitonic nucleon 
must be taken into account. A repulsion still results, but it is scalar in nature as 



found in [161| and for which there are direct physical indications [162 . 



10.5.2 Kaons at GSI 

The i^'"^and i^'" energies in the top-down scenario are given by 



uj± = ±—Vk + Jk^ + (110.25) 

with Vx given in (^.27) below. Although at high densities it will decouple, the term 
linear in Vk that figures in the range correction in {IlKN)eff will give a slightly different 
effective mass for i^+and /('"before decoupling. Although the large distance vector mean 
field must arise from vector meson exchange, this must be cut off at a reasonably large 
distance, say, ~ 0.5 fm as indicated by the baryon-flow mentioned above. 

For the GSI experiments with temperature ~ 70 MeV, the nucleon and kaon mo- 
menta are ~ 444 MeV and ~ 322 MeV, respectively, and 



Mp) ~ 0.82. (0.26) 



We therefore propose to use 



7t/vVW(p = 0)/mK ~ 0.15 (0.27) 

This is small. Furthermore we assume it to be constant above po- This assumption 
amounts to taking the vector coupling to drop as ~ 1/p. 
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The above arguments could be quantified by a specific model. For example, as 
alluded to above, the low-lying p- and cu-excitations in the bottom-up model can be 
built up as A^*-hole excitations [ |132|| . At higher densities, these provide the low-mass 
strength. One might attempt to calculate the coupling constants to these excitations 
in the constituent-quark (or quasi-quark) model, which as we have suggested would be 



expected to be more applicable at densities near chiral restoration. Riska and Brown [133 



find the quark model couplings to be a factor ~ 2 lower than the hadronic ones 163 



10.6 Schematic Model 

10.6.1 First try with the simplest form 

On the basis of our above considerations, a first try in transport calculation might 



use the vector potential with the Song scaling ||155|| as g^/m*^ = constant and the effective 
mass 



2 _ ™2 1 PK^KN)eff 



- I ^ 



10.28 



K 



with T^KN = 400 MeV and {T,KN)eff given by eq. ([10|.22). While as argued above the 
vector coupling will decouple at very high densities, as ujk drops, the vector potential will 
become less important even at moderate densities since the factor ujk/^k comes into the 
coupling of the vector potential to the kaon. 



Our schematic model (10.28) gives roughly the same mass as used by Li and 
Brown ||146|| to predict kaon and antikaon subthreshold production at GSI. For p ~ 3po 
it gives ~ 230 MeV, less than half the bare mass. We predict somewhat fewer K~- 
mesons because the (attractive) vector interaction is largely reduced if not decoupled. 
Gassing et al ||1 50|| have employed an somewhat lower than that given by eq.(p!0|.28) 
to describe a lot of data. 

The experimental data verify that our description is quite good up to the densities 
probed, i.e., ~ 3 po- In order to go higher in density, we switch to our top-down description 
which through eq. (|l0| .14) involves g'^/m^J'. Although the scalar interaction could have 
roughly the same form factor as the vector, cutting it off at ~ 0.5 fm as mentioned above, 
we believe that this will be countered by the dropping scalar mass m* which must go to 
zero at chiral restoration (viewed as an order parameter). Treating the scalar interaction 
linearly as a fluctuation (as in (|10[14)) cannot be expected to be valid all the way to chiral 
restoration but approaching the latter the a-particle becomes the "dilaton" in the sense 
of Weinberg's "mended symmetry" 0, ^ with mass going to zero (in the chiral limit) 
together with the pion . 

At high densities at which the vector interaction decouples, the and K'^ will 
experience nearly the same very strong attractive interactions. This can be minimally 
expressed through the effective mass m^. At low densities where the vector potential not 
only comes into play but slightly predominates over the scalar potential, the will have 
a small repulsive interaction with nucleons. It is this interaction, extrapolated without 
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medium effects by Bratkovskaya et al |p.42|| which gives the long equihbration time of 40 
fm/c. However, clearly the medium effects will change this by an order of magnitude. 



10.6.2 Implications on kaon condensation and maximum neutron-star mass 

While in heavy-ion processes, we expect that taking m* to zero (or nearly zero in 
the real world) is a relevant limiting process, we do not have to take m* to zero for kaon 
condensation in neutron stars, since the i^~-mass must be brought down only to the 
electron chemical potential fie — Ep{e), the approximate equality holding because the 
electrons are highly degenerate. Arguments based on effective chiral Lagrangians that 
are consistent with low-energy phenomenology in kaon-nucleon and kaon-nuclear systems 
typically give the critical density in neutron-star matter of pc ~ 3po l|164| , |165| , |154|| . We 



should however mention that it has been suggested that the electron chemical potential /le 
could be kept low by replacing electrons plus neutrons by hyperons (or more generally 
by exploiting Pauli exclusion principle with hyperon introduction) in neutron stars ||166| . 
In this case, the fie might never meet m^. 

Hyperon introduction may or may not take place, but even if it does, the scenario 
will be more subtle than considered presently. To see what can happen, let us consider 
what one could expect from a naive extrapolation to the relevant density, i.e., p ~ 3po; 
based on the best available nuclear physics. The replacement of neutron plus electron will 
take place if the vector mean field felt by the neutron is still high at that density. The 
threshold for that would be 

E^ + VN + Pe^ Ms- + "^Vn + Sj:- (0.29) 

where Ep is the Fermi energy of the neutron, M^- the bare mass of the S~ and the 
5's- the scalar potential energy felt by the S~. Here we are simply assuming that the 
two non-strange quarks of the E~ experience 2/3 of the vector mean field felt by the 
neutron. Extrapolating the FTSl theory | 156|| P^ and taking into account in Vn the effect 
of the p-meson using vector dominance , we find Ep + Va^ ~ 1064 MeV at p ~ 3po- 



From the extended BPAL 32 equation of state with compression modulus 240 MeV ||167| 



the electron chemical potential comes out to be /ig — 214 MeV. So the left-hand side of 
(1^.29) is Ep + Vn + Pe ^ 1278 MeV. For the right-hand side, we use the scalar potential 



energy for the E at p ~ 3po estimated by Brown, Lee and Rapp ||168|| to find that 



Ms- + I^AT + 5*2- ~ 1240 MeV. The replacement of neutron plus electron by S~ looks 
favored but only slightly. 

What is the possible scenario on the maximum neutron star mass if we continue 
assuming that the calculation we made here can be trusted? A plausible scenario would 
be as follows, if "condensation supposedly occurs at about the same density and both 
the hyperonic excitation (in the form of SA^~^ - where stands for the nucleon hole 

^^There is nothing that would suggest that the effective Lagrangian vahd up to p ~ po will continue to 
be valid at p ~ Spo without addition of higher mass-dimension operators, particularly if the chiral critical 
point is nearby. So this exercise can be taken only as indicative. 
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- component of the "kaesobar" [|168|| ) and condensation would occur at T ~ 50 MeV 
relevant to the neutron-star matter. Now if as is likely the temperature is greater than 
the difference in energies between the two possible phases, although the hyperons will be 
more important initially than the kaons, all of the phases will enter more or less equally 
in constructing the free energy of the system. In going to higher density the distribution 
between the different phases will change in order to minimize the free energy. Then it is 
clear that dropping from one minimum to another, the derivative of the free energy with 
density - which is just the pressure - will decrease as compared with the pressure from 
any single phase. This would imply that the maximum neutron star mass calculated with 
either hyperonic excitation or kaon condensation alone must be greater than the neutron 
star mass calculated with inclusion of both. 

The story will be quite different if the vector field decouples. We showed in Section 



10.5.1 that the isoscalar vector mean field must drop by a factor > 9 in the change- 



over from nucleons to quasiquarks as variables as one approaches the chiral restoration 
density. Hyperons will disappear during this drop. It is then highly likely that the kaon 
will condense before chiral restoration and that the kaon condensed phase will persist 
through the relevant range of densities which determine the maximum neutron star mass. 



10.7 Discussions 

In this section, we have given arguments to support that not only the kaon mass 
but also its coupling to vector mesons should drop in matter with density. In particular, 
with the introduction of medium effects the apparent equilibration found in strangeness 
production at GSI can be increased from the baryon number density of ~ ^po up to 
the much more reasonable ~ 2po- These properties do not appear at first sight to be 
connected to BR scaling per se. But they must be connected in an albeit indirect way. 

From the baryon flow analysis we have direct indications that the vector interaction 
decouples from the nucleon at a three-momentum of |p| ~ 0.9 GeV/c or at roughly 0.5 
m]s[C. In colliding heavy ions this is reached at a kinetic energy per nucleon of ~ ^m]^(? 
which means a temperature of 78 MeV when equated to |T. This is just the temperature 
for chemical freeze-out at GSI energies. We have given several theoretical arguments why 
the vector coupling should drop rapidly with density. 

Once the vector mean field, which acts with opposite signs on the and K~ 
mesons is decoupled, these mesons will feel the same highly attractive scalar meson field. 



Their masses will fall down sharply; e.g., from eq.([10|.28) with proposed parameters. 



rriK 



0.5 dH-SO) 



at p ~ 3po and possibly further because of the dropping m*. The differing slopes of 
and K~ with kinetic energy will then develop after chemical freeze-out, as suggested by 
Li and Brown ID-461 . 



In this section we have focused on the phenomenon at GSI energies. Here the 
chief role that the dropping i^'-mass played was to keep the combination fiB + ^*k- 
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nearly constant so that low freeze-out density in the thermal equilibration scenario became 
irrelevant for the K'^ / K~ia.tio. We suggest that the same scenario applies to AGS physics, 
where the freeze-out density in the thermal equilibration picture comes out to be ~ 
0.35po ||169|| . In fact, there is no discernible dependence on centrality in the K~ / K^idXios 
measured at AA GeV, QA GeV, 9>A GeV and 10. GeV [ p70[] . From this it follows either 
that the ratio of produced K~io hyperons is nearly independent of density or that the 
negative strangeness equilibrates down to a lower freeze-out density and then disperses. 
Given the relative weakness of the strange interactions, we believe the former to be the 
case. In fact, we suggest that near constancy with multiplicity of the K~ / K^iaXio found 
experimentally be used to determine temperature dependence of tu^- in the region of 
temperatures reached in AGS physics. As was done for GSI energies, the temperature 
can be obtained from the inverse slopes of the kaon and antikaon distributions of and 
then the temperature dependence of tu^. can be added to the density dependence as in 
147| , in such a way that hb + ^k- stays roughly constant as function of density. At least 



this can be done in the low-density regime considered in ||169|| when the approximation 
of a Boltzmann gas is accurate enough to calculate fiB- Our "broad-band equilibration"; 
i.e., the production of the same, apparently equilibrated ratio of /('"-mesons to hyperons 
over a broad band of densities, avoids complications in the way in which the K~ degree 
of freedom is mixed into other degrees of freedom at low density ||171|| . Most of the K~- 



production will take place at the higher densities, as shown in Fig. [I2|, where the degrees 
of freedom other than K~ have been sent up to higher energies by the Pauli principle. 

Unless the electron chemical potential in dense neutron star matter is prevented from 
increasing with density (as might happen if the repulsive E-nuclear interaction turns to 
attraction at p > po), kaon condensation will take place before chiral restoration. This 
has several implications at and beyond chiral restoration. For instance, its presence would 
have influence on the conjectured color superconductivity at high density, in particular 
regarding its possible coexistence with Overhauser effects, skyrmion crystal and other 
phases with interesting effects on neutron star cooling. 

The phenomenon of vector decoupling, if confirmed to be correct, will have several 
important spin-off consequences. The first is that it will provide a refutation of the 



recent claim [172] that in the mean-field theory with a kaon- nuclear potential given by 



the vector-exchange (Weinberg- Tomozawa) term - both argued to be valid at high density 
- kaon condensation would be pushed up to a much higher density than that relevant in 
neutron-star matter. Our chief point against that argument is that the vector decoupling 
and the different role of scalar fields in QCD (e.g., BR scaling) described in this article 
cannot be accessed by the mean field reasoning used in ||172|| or by any other standard 
nuclear physics potential models. The second consequence that could be of a potential 
importance to the interpretation of heavy-ion experiments is that if the vector coupling 
rapidly diminishes with density, the strong-coupling perturbation calculation of the vector 



response functions used in terms of "melting" vector mesons to explain | 129 |, e.g., the 
CERES dilepton data must break down rapidly as density increases. This could provide a 
specific mechanism for the quasiparticle description of BR scaling found to be successful 
in nuclear matter []7S|). It is of course difficult to be quantitative as to precisely where 
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this can happen. 

Finally we should stress that the scenario presented in this section - which is an- 
chored on Brown-Rho scaling- should not be considered as an alternative to a possible 
quark-gluon scenario currently favored by the heavy-ion community. It is once again more 
likely a sort of complementarity of the same physics along the line that the Cheshire-Cat 
Principle embodies that would continue to apply at higher energies. 



11 EFFECTIVE FIELD THEORIES FOR DENSE 
QCD 

Thus far we have been climbing up in density in what may be aptly called "bottom- 
up" approach. In this section, we will start at some asymptotic density and come down 
to the density regime relevant to dense compact stars and possibly in certain regimes 
of heavy-ion collisions. This will be a "top-down" approach exploiting once more how 
chiral symmetry in nature manifests itself through (Nambu-)Goldstone mode. There is 
an intense activity on this matter, comprehensively reviewed in the recent literature ||1 73|| , 



so we will not go into most of the topics. Here we will focus on the most intriguing 
aspect of dense QCD that reflects an intriguing repeat of the Cheshire-Cat phenomenon. 
We wish to highlight the seemingly continuous role of chiral symmetry from low-density 
strong-coupling regime to high density weak-coupling regime in the form of a skyrmion 
replica that maps high-density excitations to low-density excitations. In view of the rapid 
development on this matter, our discussion will necessarily be incomplete and primitive 
but the general feature of the Cheshire Cat phenomenon is likely to survive the test of 
the time. 

11.1 Color-Flavor Locking for A^^;^ = 3 

11.1.1 Instability at high density 

Asymptotic freedom of QCD suggests that at very high density with yU — > cxd, the 
QCD gauge coupling becomes very weak and so interactions disappear. However if there 
is an attraction in a certain channel, now matter how weak, the same renormalizationg 



group consideration discussed above in connection with Landau Fermi liquid state 
tells us that the system will be unstable against a phase change. Indeed in the weak 
coupling limit, one gluon exchange that survives asymptotic freedom has an attraction 
in the channel that involves diquarks sitting on the opposite side of the Fermi sphere, 
namely, the BCS channel and this inevitably leads to the Cooper pair condensation as 
in superconductivity. That this must occur at asymptotic density is without doubt as 
dictated by the renormalization group argument. Whether or not - and in what manner 
- this will occur in a physically relevant regime of density cannot be addressed by the 
asymptotic QCD and it appears most plausible that the scenario is much richer in various 
different ways than just the Cooper pairing type as one approaches top-down to the 
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density regime of relevance. This matter is amply discussed in the literature ||173|| but the 



situation is totally unclear. It is not our aim to dwell on the multitude of possibilities. We 
will focus on one particular aspect of this phenomenon, that is, the situation where color 
and flavor of QCD get locked, i.e., the color- flavor- locked (CFL) phase which renders the 
discussion simple and transparent. 

11.1.2 Symmetry breaking and excitations 

We shall discuss first the case of = Np = 3 which is relevant at asymptotic 
density, returning to the case of A^c = 3, Np = 2 in the next subsection. 

Gluon exchange between two quarks is attractive in the color 3 channel and repulsive 
in the 6 channel. So one expects diquarks would condense in the 3 channel. We shall 
ignore a small mixing to the 6 channel which is allowed by symmetry. The diquark 
condensate would be in the form 

= -ilRal^Rf^) = '^^''^"''^api (0.1) 

where k is a constant, i,j are SL{2,C) indices, a,b are color indices and a, (3 are flavor 
indices. This condensation breaks the symmetry 

G^H ([n].2) 

with 

G = SUi3)cxSUi3)LxSUi3)RxU{l)B, (11.3) 

H = SUISUl+rx Z2. (|11].4) 

At superdensity, instantons are suppressed, but the anomaly never disappears, so the 
U{1)a symmetry is never truly restored. We shall not be concerned with it here. The 
color is completely broken and since the vectorial color is locked to both L and R, chi- 
ral symmetry is also broken so the resulting invariance group is SU{3)c+l+r- Since 
(llalLis^Ra'lRb) ~ ^^"^ ^a/SF^'^^^ IS gaugc invariant, it can be used as an order parameter but 
while it breaks chiral symmetry, it leaves Z2 invariant, thus its appearance in H. 

The spontaneous symmetry breaking induces Goldstone bosons: eight from the 
breaking of the color SU{3)c, eight from chiral symmetry breaking and one from U{1)b —>■ 
We shall forget the last one since it does not concern us here. The striking feature 
we would like to focus on here is the uncanny resemblance of the spectrum in the CFL 
phase to the that of low density in terms of a hidden local symmetry |Q. The scalar 



Goldstones are eaten up by the gluons thereby the vector bosons becoming massive. The 
octet of the massive gluons in the CFL are the analogs to the (octet) massive light-quark 
vector mesons in the low-density vacuum. The octet pseudoscalar Goldstones in the CFL 
are the analog to the octet pions in the low density vacuum, the dynamics of which can be 
described by an analogous chiral Lagrangian. The quarks in the system are gapped - and 
hence massive - and bear the quantum numbers of the octet baryon of the zero-density 
vacuum. We shall show below that this can be understood by looking at the solitons 
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in the chiral Lagrangian [ |174|| . Since the sohton corresponds to a gapped colored quark, 



we call it "qualiton" in analogy to the qualiton considered by Kaplan for the constituent 
quark ||175|| . (The breaking of U{1)b leads to superfluidity which resembles the same 
in nuclear matter, another uncanny correspondence.) The one-to-one matching between 
the low density spectrum in terms of hadronic variables and the high density spectrum 
in terms of quark-gluon variable is referred to as "quark-hadron continuity." |p.76|| . We 



will identify this as belonging to the class of the Cheshire Cat phenomenon in the strong 
interaction physics [11]. 

11.1.3 Chiral Lagrangians and qualitons 

As suggested in ||174| , |177|| , the dynamics of the surviving Goldstone modes can be 
described by an effective chiral field theory as in the zero-density situation. We introduce 
the chiral effective field 

ibaaix) ~ e'^eabc^ap^qLi{-'"F,x)qj]j{VF,x) ([n].5) 

corresponding to the map SU{?>)c x SU{?>)l/ SU{?>)c+l and similarly for i?, 

iRaa{x) ~ e'^eabc(^af3'rqRi{-VF,x)q'flj{VF,x) (|ll|.6) 

corresponding to the map SU{3)c x SU{3)r/ SU{3)c+r. 

Under an SU (3)c x SU (3) l x SU{3)ji transformation by unitary matrices {gc, gL, gn), 
transforms as ^ Qc^igL and C,r transforms as i— >• g*^Rgli- In the ground state of 
the CFL superconductor, and take the same constant value. QCD symmetries with 
(g.l) imply 



i^Laa) = - i^Raa) = f^^aa- ([n].7) 



The Goldstone bosons are the low-lying excitations of the condensate, given as unitary 
matrices ^l(x) = gJ{x)gL{x) and ^r{x) = g^{x)gR{x). For the present decomposition, we 
note the extra invariance under the (hidden) local transformation gc+L+R{x) - which is 



an analog to the hidden gauge transform h{x) of - within the diagonal SU{3)c+l+r 
through 



9l{x) gJ{x)gl+L+R{x) 9l,r{x) 9c+L+R{x)gL,R{x) (|TT].8 



Hence, the spontaneous breaking of SU{3)c x {SU(3)l x SU{3)r) — > SU{3)c+l+r can be 
realized non-linearly through the use of ^l,_r(x) or linearly through the use of gc,L,R{x) 
with the addition of an octet vector gauge field transforming inhomogeneously under local 
gc+L+R{x)- This can be identified with the hidden local symmetry |^4[ Q. Clearly the 

^^One has to be careful with a singularity in the product of two fermion fields at one point. See [ |174| 
for a more proper definition. 

^^We comment as a side remark on a technical detail which is somewhat outside of the scope of this 
review but may be helpful to those interested in subtleties involved in hidden-gauge-symmetry aspects of 
the problem. 
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composites carry color-flavor in the unbroken subgroup, with a mass of the order of the 
superconducting gap. 

The current associated with det^^ (det^/j) is the left (right)-handed U{1)l {U{1)r) 
baryon number current. Because of the U{1) axial anomaly, the field det det C,r is mas- 
sive due to instantons. We decouple the massive field from the low energy effective action 
by imposing det^^/ det^R = 1. The Goldstone boson associated with spontaneously bro- 
ken U{1)b symmetry is described by det .^L'det ^r and responsible for baryon superfluidity. 
But since it is not directly relevant for our problem, we further choose det - det = 1 to 
isolate the Goldstone bosons resulting from the spontaneous breaking of chiral symmetry, 
from the massive ones eaten up by the gluons. We now parameterize the unitary matrices 



As mentioned, when the color is completely broken, the octet gluons become massive by Higgs mech- 
anism with the mass proportional to the chemical potential fi. These were found to have the same 
quantum numbers as those of the light-quark vector mesons present at zero density. The hidden gauge 
bosons excited at high density discussed here belong to the unbroken diagonal subgroup SU{3)c+l+r- 
How are these vectors related to the massive gluons? 



In |17S|, the hidden gauge bosons £ SU{3)c+l+r were identified with the spin-1 two-quark bound 
states with the properties that match with the Georgi vector limit at some non-asymptotic density. At 
the leading order in these states seem to have nothing to do with the Higgsed gluons: The overlap 
is zero. But this cannot be correct. In nature, the low-lying vector excitations at large density must be a 
coherent mixture of the two but at the leading order, they must be equivalent. This is analogous to the 
complementary description of the pion as a Goldstone mode and as a zero- mass bound state of a quark and 
antiquark. Indeed th is comple mentarity has also been verified in the pion channel in dense QCD in terms 



of bound diquarks [179, 17S|. To see this aspect in the vector channel, we follow the CCWZ (Callan- 



Coleman-Wess-Zumino) formalism |18C]. Consider the symmetry group G — SU{3)l x SU{3)r x SU{3)c 
broken down to H ~ SU (3) l+r+c- Using a notation slightly different from what is used in this subsection, 
let the G transformation be represented in the block diagonal form g — diag(L, R, G) corresponding 
to the left, right and color transformations. Let the generator of the diagonal subgroup be given by 
h = diag([/, J7, J7). Now we would like to write € G, the rotation matrix that transforms the 

standard vacuum configuration to the local field configuration. One choice consistent with the CCWZ 
prescription is 



E;(x) = diag(CL(a;), S,r{x), 1) 



which transforms as 



gS/i-i =diag(LaC/"\ R^rU-\ GU-'). 
It is clear that we are required to take 

GU-^ ^ 1 

which means that U — C and that S,l,r transform as 

^R R£,rC ^■ 

Now U \s a. non-linear function of L, _R, and ^^j, and so must be G. In HGS, this is the group that is 
gauged. Therefore the "hidden gaug e symmetry" can be equated to the color symmetry of QCD as was 



done by Casalbuoni and Gatto |177|. This can be seen also in Eq.(^.10) given below. 
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^L,R as 

= exp {tntT^/f) , U^) = exp {tn^T^/f) , (0.9) 

where are SU{3) generators, normalized as Tr {T^T^) = 6^^/2. The Goldstone 
bosons TT^ transform nonhnearly under SU{3)c x SU{3)l x SU{3)r but hnearly under the 
unbroken symmetry group SU{3)c+l+r- 

The effective Lagrangian for the CFL phase at asymptotic densities follows by inte- 



grating out the 'hard' quark modes at the edge of the Fermi surface ||181|| . The effective 
Lagrangian for Ul^^ is a standard non-linear sigma model in D=4 dimensions including 
the interaction of Goldstone bosons with colored but "screened" gluons G. Expanding in 
powers of derivative, the effective Lagrangian for the Goldstone bosons is then 

C= + ^ Tr(9oeL5oei) - ^ Tr{d,^Ldr^i) + g^G ■ -h + {L ^ R) + ■ ■ ■ (^.10) 

Here the ellipsis stands for higher derivative terms including the Wess-Zumino-Witten 
term p3| that is needed to account for anomalies. The vectorial coupling to the gluon 



field Gfj, via the color current indicates that the left-right chiral symmetry is also locked 
to each other. In the presence of a chemical potential, the SL{2, C) breaks down to 0(3), 
so the "decay constant" / has two different values for the time and space components. 
The temporal and spatial decay constants /t,5 are fixed by the 'hard' modes at the Fermi 
surface. Their exact values are determined by the dynamics at the Fermi surface. They 
go like /5 ~ /j. ~ /i. 

The effective Lagrangian (|ll|.10) in the CFL phase bears much in common with 



Skyrme-type Lagrangians |Tl| , with the screened color mediated interaction analogous to 
the exchange of massive vector mesons. Indeed in the CFL phase the 'screened' gluons 
and the 'Higgsed' gauge composites are the analog of the massive vector mesons in the 
low density phase. In the CFL phase the WZW term is needed to enforce the correct 
flavor anomaly structure. 

The 'Higgsed' gluons are very massive in the large /i limit, so can be integrated 
out from the Lagrangian (|lT[lO). This is very much like the hidden gauge symmetry 
theory wherein the vector mesons can be integrated out to give rise to a nonlinear sigma 



model. The theory then can have a soliton as suggested by Hong, Rho and Zahed ||174|| if 
there are higher derivative terms that stabilize it. The qualiton will be colored. However 
the pair-condensed ground state is colored also, so combined with the background, the 
excitation will be effectively color-singlet. In other words, the chiral field from which the 
soliton arises can be considered as color-singlet fluctuation. This means that one can 
rewrite the effective field theory in terms of the color-singlet field 

Ul^iLa^' (0.11) 

which transforms U — > giUgli- The Lagrangian (|Tl|.10) can then be rewritten in the form 
familiar from low density 

C = ^Tr{doUdoU^) - ^TT{diUd,U^) + 0{d^) + 0{M^) + ■■■ (^.12) 
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The resulting theory is essentially identical to the familiar nonlinear sigma model giving 
rise to the Skyrme soliton except that the mass terms which we did not specify 



here should be quadratic in the quark mass matrix Ai because of the Z2 invariance. 
The qualiton can be collective-quantized in a way paralleling the low-density skyrmions. 
This was worked out recently in ||182|| in which the authors argue that since a qualiton 



is a quark on top of the condensed state, it would correspond to the quasiparticle of the 
particle-hole complex and the qualiton mass that is obtained as a soliton solution must 
therefore correspond to 2A plus the ground-state (or condensation) energy where A is 
the gap. 

11.2 Complementarity of Hidden Gauge Symmetry and Color 
Gauge Symmetry and BR Scaling 

The Nc = 3, Np = 2 case is a bit more intricate and involved than the A^^^ = Np = 3 
CFL discussed above. It turns out, however, that color-flavor locking in this case |]183|] 



is more relevant in the density regime closer to nature, both in high density and in low 
density. It can also be mapped as we shall argue to hidden local symmetry discussed in 
the previous sections. 

11.2.1 Quark number susceptibility 



It was observed in [T^, |184|| that the quark number susceptibility x± = {d/d^u ± 



d/djj,d){puipd) where pu,d and iJ,u,d are respectively u, (i-quark number density and chem- 



ical potential measured on lattice as a function of temperature ||185|| exhibited a smooth 
change-over from a flavor- gauge symmetry or hidden gauge symmetry to QCD color- gauge 
symmetry at the chiral phase transition critical temperature Tc . It was suggested [|l^ 
that at the phase transition, the flavor gauge symmetry - which is induced and hence not 
fundamental - gets converted directly to the color gauge symmetry - which is fundamen- 
tal, implying that they could be related in an intricate way. In this subsection, following 
|186|| , we suggest how this can be realized in terms of color-flavor-locked (CFL) quark- 



antiquark and diquark condensates and "vector manifestation" (VM) of chiral symmetry 
explained in Section ^. It will be seen how BR scaling can be flt into the general scheme 
that results from these developments: Its proof will be the most convincing case for the 
validity of BR scaling. 

Our argument relies on two recent remarkable developments that come from seem- 



ingly unrelated sectors. One is the suggestion by Harada and Yamawaki ||12| that the 
phase transition from the Nambu-Goldstone phase to the Wigner-Weyl phase involves 
"vector manifestation" of chiral symmetry, that is, at the phase transition, the longi- 
tudinal components of the light-quark vector mesons (i.e., the triplet p in the 2-flavor 
case) and the triplet pions (vr") come together becoming massless in the chiral limit, the 
vectors decoupling a la Georgi's vector limit |Q but with the pion decay constant 
vanishing at that point. The other important development is the proposal by Berges and 



Wetterich ||187|, |183|| that color and flavor get completely locked in the Nambu-Goldstone 
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phase (a) for three flavors {Nf = 3) ||187|| by the quark- ant iquark condensate in the color- 
octet (8) channel 

X = {q:.j:{r^U{X^r%) (0.13) 

i=l 

and (2) for two flavors {Np = 2) ||183|| by (|Tl].13) together with the diquark condensate 
in the color-antitriplet (3) channel 

A={q:{T,U{X,r%). (0.14) 

In (^.13) and (|11].14), the indices a,P denote the flavors and a,b the colors. In the 
Berges-Wetterich scenario, the chiral phase transition and deconfinement occur through 
the melting of the condensates x ^i^nd A. 

11.2.2 Color-isospin locking 

The scenario for color-flavor-locking (CFL) is a bit different depending on the num- 
ber of flavors Np. For our purpose, the case of two flavors is more relevant, so we shall 
focus on this case. In ||183|| , Berges and Wetterich argue that both the x ^ conden- 



sates can be nonzero in the vacuum, thereby completely breaking the color. Should one 
of the condensates turn out to be zero, then color would be only partially broken (see 



1 1 73|| for a simple explanation for this and references). This pattern of color breaking 
and color-flavor locking renders the octet gluons and six quarks massive by the Higgs 
mechanism and generates three Goldstone pions due to the broken chiral symmetry. All 
of the excitations are integer-charged. Among the eight massive gluons, three of them are 
identified with the isotriplet p's with the mass 

rrip = KQcX dni-lS) 

where k is an unknown constant and Qc the color gauge coupling. The fourth vector meson 
is identified with the isosinglet uj with the mass 



k'qcA (0.16) 



where k' is another constant. The remaining four vector mesons turn out to have exotic 
quantum numbers and are presumably heavy. We assume that they decouple from the 
low-energy regime. As for the fermions, there are two baryons with the quantum numbers 
of the proton and neutron with their masses proportional to the scalar condensate 0, 

0=(«)- (0-17) 

The four remaining fermions are also of exotic quantum numbers with zero baryon number 
and heavier, so we assume that they also decouple from the low-energy sector. What 
concerns us here is therefore the three pions, the proton and neutron, the p-mesons and 
the a;-meson. 
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11.2.3 Implications of the lattice measurements of quark number susceptibility 



It was argued in [^] that the "measured" singlet and non-singlet QNS's |p.85|| in- 



dicate that both the p and cu couplings vanish at the transition temperature T^. This 
meant that the ujNN coupling which is ~ 3 times the pNN coupling at zero temperature 
became equal to the latter at near the critical temperature. This also meant that the 
both vector mesons became massless and decoupled. Viewed from the CFL point of view, 
it follows from (|lT].15) and (|Tl|.16) that the condensates x ^i-nd A "melt" at that point. 
This is consistent with the observation by Wetterich ||187|| that for three-flavor QCD, the 



phase transition - which is both chiral and deconfining - occurs at Tc with the melting of 
the color-octet condensate x- The transition is first-order for Np = 3 in agreement with 
lattice calculations, so the vector meson mass does not go to zero smoothly but makes 
a jump from a finite value to zero. We expect however that in the case of Np = 2 the 
transition will be second-order with the vector mass dropping to zero continuously. 

Noting that both the CFL formulas ([TT1.15)-(|TTil6) and the vector manifestation 
result (H.9) are of the Higgsed type, we invoke the lattice results to arrive at 



agvfn ~ KQcX ~ /t'^cA. (M.18 



We admit that this relation follows neither from the group-theoretical considerations of 
Berges and Wetterich |p.83| , |187|| nor from the vector manifestation of HLS [1^. We are 



proposing that this is indicated by the lattice data, at least near the phase transition 
point. It holds empirically at zero temperature and zero density, so we are led to assume 
that it holds at least approximately from T = up to T = Tc 0. It should be noted 
that the vanishing of the hidden gauge coupling gv corresponds to the vanishing of the 
condensates x ^i^nd A with the color gauge coupling remaining non- vanishing, Qc ^ 0. 

Next, we have shown in that above the chiral transition temperature T > 
the QNS's can be well described by perturbative gluon exchange with a gluon coupling 



,2 



constant |^ ~ 0.19 and argued that the flavor gauge symmetry cedes to the fundamental 



QCD gauge symmetry. Now the HLS theory is moot on what it could be beyond the chiral 
restoration point since the theory essentially terminates at Tc. We propose that this is 



where the color-flavor locking of |p.87| , |183|| phrased in the QCD variables takes over by 
supplying a logical language for crossing-over from below Tc to above Tc. Indeed (]TT|.18) 
describes the relay that must take place in terms of the hidden flavor gauge coupling 
gv on one side and the color gauge coupling gc on the other side. Now above Tc, the 



^^This would imply that the nonet symmetry for Np = 3 or the quartet symmetry Np = 2 is a good 
symmetry not only at T = but also for T 0. Why this symmetry should hold at any non-zero 
temperature or density is not obvious either in CFL QC D or in HLS effective theory. In discussing color- 
flavor locking in two-flavor QCD, Berges and Wetterich |l83|] entertain among others the possibility that 
the two condensates x and A could be different at the critical point. As for the HLS theory, at one-loop 
order, the RG flows are expected to be different for the p and w properties, so it is not obvious that the 
p and Lo mesons would reach the chiral restoration with the Georgi vector limit at the same temperature. 
Nonetheless if our interpretation of QNS is correct, it seems most plausible and appealing that the nonet 
or quartet symmetry does hold at the phase transition. Proving this conjecture remains as a theoretical 
challenge. 
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color and flavor must unlock, with the gluons becoming massless and releasing the scalar 
Goldstones. The dynamics of quarks and gluons in this regime will then be given by hot 
QCD in the proper sense. The way the two condensates melt as temperature is increased 
is a dynamical issue which seems to be difficult to address unambiguously within the 
present scheme. It will have to be up to lattice measurements to settle this issue. Our 
chief point here is that their melting is intricately connected. 



11.2.4 Link to BR scaling and Landau parameter Fi 

The situation appears to be quite different in dense medium. Since there is no 
guidance from lattice as it is impossible at present to put density on lattice except for 
unphysical cases of two colors or adjoint quarks, we shall simply assume that the above 
scenario holds in density up to n = nc Certain models indicate that the phase structure 
near chiral restoration could be quite involved and complex. As suggested by Schafer 
and Wilczek [ |176|| , an intriguing possibility is that the three-flavor color-flavor locking 



operative at asymptotic density continues all the way down to the "chiral transition 
density" (nc) in which case there will be no real phase change since there will then be 
a one-to-one mapping between hadrons and quark/gluons, e.g., in the sense of "hadron- 
quark continuity." However the non- negligible strange-quark mass is likely to spoil the 
ideal three-flavor consideration. One possible alternative scenario is that viewed from 
"bottom-up," one gets into the phase where x = o" = and A 7^ corresponding to 



the two-flavor color superconducting (2csc) phase ||173|| . Unless A goes to zero at pc, this 
would mean that the p mesons become massless but the uo meson remains massive. One 
cannot say that this is inconsistent with the vector manifestation since the HLS does not 
require that U{2) symmetry hold at the chiral restoration point or in medium in general. 
But this seems unlikely. On the other hand, it is highly plausible that both x ^"^^ ^ 
approach zero (or near zero if it is first-order) from below ric and then A picks up a 
non- zero value at or above ric in which case we will preserve the mass formula (§.9) as 
one approaches ric- This is the scenario that we favor. 

Among the scaling relations implied by BR scaling 0, the one most often 
discussed in the literature is the dropping of the p-meson mass in medium. This relation 
has been extensively discussed recently in connection with the CERN-CERES data on 
dilepton production in heavy-ion collisions. The simplest explanation for the observed 
dilepton enhancement at an invariant mass ~ 400 MeV is to invoke BR scaling for the 
excitations relevant in the process [ |188| |. It turns out however that this explanation is 
not unique. One could explain it equally well if the p meson "melted" in dense medium 
with a broadened width [[L29|] . Since the process is essentially governed by a Boltzmann 
factor, all that is needed is the shift downward of the p strength function: the expanding 
width simply does the job as needed for the dilepton yield. If one calculates the current- 
current correlation function in low-order perturbation theory with a phenomenological 
Lagrangian, it is clear, because of the strong coupling of the p meson with the medium, 
that the meson will develop a large width in medium and "melt" at higher density. The 
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In this subsection, we again denote density by n reserving p for the vector meson. 
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upshot of the dilepton experiments then is that they cannot distinguish the variety of 
scenarios that probe average properties of hadrons in the baryon density regime - which 
is rather dilute - encountered in the experiments. 

In the vector manifestation scenario, the width should become narrower, decreasing 
like ~ Qy. Then the vector mesons become more a quasiparticle at high density than at 
lower density. This is the underlying picture of BR scaling. 

Thus far, we have made a link between the color-flavor-locked condensates and hid- 
den gauge symmetry. We can go even one step further and via BR scaling, make an 
intriguing connection between QCD "vacuum" properties and many-body nuclear inter- 
actions. This connection was discussed in Section |[ Simply put, it comes about because 



nuclear matter owes its stability to a Fermi-liquid fixed point ^7^. Certain interesting 



nuclear properties were found to be calculable in terms of the Fermi-liquid fixed point 



parameters ||79[. Among others, it was shown that the Landau parameter Fi - which is 
a component of quasiparticle interactions - can be expressed in terms of the BR scaling 
factor $(n) = m*(n)/mp(0) 

Fi = 3(l-$-^) + Fi(7r) (|Tl|.19) 

where Fi (vr) is the contribution from the pion which is completely fixed by chiral dynamics. 
One of the most remarkable prediction was Eq.(P.4) for the anomalous gyromagnetic ratio 
6gi in nuclei 



. 4 
S9i = ^ 



T,. (0.20) 



At nuclear matter density n = no, we have -Fi(vr) = —0.153. Note that ([TT|.20) 
depends on only one parameter, $. This parameter can be extracted from various sources 
and all give about the same value, $(no) ~ 0.78. Given that this is not very accurately 
determined, it is probably a better strategy to determine $ from the data on 6gi. In any 
event, given $ at nuclear matter density, Eq.(|rT[20) makes a simple prediction. 



Sgi = 0.23x3 (11.21 



which should be compared with the measurement in the Pb region ||91[] , 



Sgf = 0.23 ±0.03. (0-22) 

We believe this to be a good agreement within the theoretical uncertainty involved. It 
is intriguing that what is an intrinsic QCD quantity can be related to what appears to 
be a standard many-body nuclear interaction summarized in the Landau quasiparticle 
parameter Fi. 

11.3 Kaon Condensation: Encore Cheshire Cat 

In addition to the "continuity" in excitations and chiral phase transitions between 
the hadronic phase and the quark phase, meson condensations can occur in high density 
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matter with n > Uc and in hadronic matter with n < ric where Uc is the presumed critical 
density for the so-called "chiral restoration" - whatever that may be. Pion condensation is 
unlikely in the density regime that is relevant for laboratory or astrophysical observations, 
but the negatively charged kaon K~ can condense at a relatively low density as well as at 
an asymptotically high density providing yet another support for "continuity" or Cheshire 
Cat. For completeness we briefly describe this phenomenon although the story is by no 
means final. 

• K~ Condensation in the Hadronic Sector 

Since this matter was discussed already in Section 0, we shall simply summarize 
the pertinent feature in a slightly different language. As originally pointed out by Ka- 
plan and Nelson | 189 | and reinterpreted by Brown et al [ 152|| , the S-wave condensation 



of K~^s is driven by "rotating away" of the kaon mass associated with both the explic- 
itly broken chiral symmetry and spontaneously broken chiral symmetry. As is widely 
discussed, hyperons can also participate through P-wave coupling with the kaons but 
near the condensation transition, the P-wave coupling would be "irrelevant" in con- 



trast to the S-wave that has to do with "relevant" terms ||190|| and hence may be ignored 
in this qualitative discussion. The condensation therefore occurs when the kaon mass 
is eaten up by attractive interactions. Negative kaons have attractive interactions with 
nucleons by exchanging the vector mesons p and u as well as scalar mesons, e.g., the a 
meson. The vector exchange is dictated by the vector current conservation and is given 
by the Weinberg- Tomozawa term and the a exchange by the T^kn ^ the KN sigma term. 
When many-body interactions through many-Fermi contact terms in the Lagrangian are 
implemented through BR scaling, one then has a Lagrangian of the form (|TU[12). One 
can generate this effective Lagrangian in chiral perturbation theory as reviewed e.g. by 



Lee | J154| ]. The condensation occurs when these relevant terms drive the system to insta- 



bility ||190|| . In neutron star matter, the kaon energy uk need not go all the way to zero. It 
suffices to drop to the electron chemical potential /ig- Since the fie increases as a function 
of density in nonrelativistic nuclear systems and the uk must fall, the crossing is bound 
to occur at some density. However exactly at what density it will occur will depend upon 
details of the dynamics and this is still a controversial issue ||172|| . In Section we have 
proposed that the critical density is rather low, ~ 3po- 

• Condensation in the CFL Phase 



The effective chiral Lagrangian discussed above, (|TT[12), implemented with mass 
terms predicts ||191| , |179| , |192| , |193| , |194[ at high density where the color-flavor locking sets 



in that the kaon mass is of the form 

m^± = cmrf(m„ + rris) (Ill|-23) 

where c is a constant that can be evaluated in the weak-coupling QCD, c = where 
A is the pairing gap and the subscripts u, d and s stand for up, down and strange 

^-'^This presumably is related to Migdal theorem. See Polchinski in Ref. ifrzl for a related discussion in 
condensed matter physics. 
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quarks respectively. Only the pairing gap A appears in the mass formula despite that an 
antiquark is involved. This is because the Majorana mass of the antiquark is equal to the 



gap |195|. Because the Goldstone bosons in this phase are of q q configuration, the kaon 



is less massive than the pion which has the mass 

m^± = cmsiniu + m^). (|il|-24) 

It is found numerically that for ^ ~ 500 - 1000 MeV, the kaon mass ranges from ~ 5 MeV 
to ~ 1 MeV. The reason for this small mass is easy to understand. First of all, because 
of the Z2 symmetry, the Goldstone mass is quadratic in the quark mass and secondly, the 
mass is proportional to (A//i)^ that goes to zero as — >• 00. 

For large density, one can ignore the attractive kaon-quark interaction analogous to 
the kaon-nuclear since the interaction is suppressed by (l//i)^. Considering the Fermi sea 
filled with non-interacting massless up, down and strange quarks, the electron chemical 
potential /Ze is found in this case to drop as ~ rUs/kp. This contrasts with the increasing 
fie found in the nonrelativistic hadronic system. Schafer [196| | found that for a reasonable 
range of parameters involved for the CFL phase, the kaon mass is always less than the 
electron chemical potential. However this does not lead to a K~ condensation of the type 
seen in the hadronic sector since as argued by Rajagopal and Wilczek ||197|| , the CFL 



phase is charge-neutral more or less independently of the mass of the strange quark and 
hence no electrons need to be present in the phase. 

Kaon condensation can however occur in superdense matter, not just because of the 
small mass of the kaons (which helps) but because of the presence of the chiral symmetry 
breaking mass term. Though perhaps different in character, this is analogous to the case of 
kaon condensation in the hadronic sector where the phase change is essentially effectuated 
by "rotating away" the large S term arising from the chiral symmetry breaking due to 
the strange-quark mass. Indeed in superdense regime, a term of the form 

- MV2/i (0-25) 

where M is the quark-mass matrix plays the role of a flavor chemical potential that 
provides an effective attraction favoring neutral kaon condensation [|198|| . 

Although interesting purely from the theoretical point of view, the physical relevance 
of this observation to the physics of compacts stars is yet to be established. Whereas 
the phenomenon lends itself to a simple and elegant analysis for asymptotic densities, 
the relevant process for the formation and structure of compact stars requires instead 
a "bottom-up" approach starting from low density at the stage of supernovae explosion 
and climbing up to high but non-asymptotic density as the matter is compressed in the 
interior of compact stars. This must involve, along the way, such hadronic phases as kaon 
(or pion) condensation, hyperon presence etc. for which QCD is intractable and it is not 
clear that the simple picture based on QCD at asymptotic density is directly relevant 
to what actually takes place in the stellar matter. Working out the change-over from 
hadronic variables to QCD variables necessary for a realistic description of the process - 
which we believe involves a "Cheshire cat" mechanism - remains a challenge for nuclear 
theorists in the sense discussed in this review. 
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